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CHAPTER 1 
THE MOTIVATION FOR ANALOG COMPUTERS 


Analog Computation, based on the modern electronic 
analog computer, is of fairly recent date. The first 
commercially available general purpose electronic 
analog computers appeared on the market in the 1940's. 
These early machines were an outgrowth of an emerging 
electronics technology anda critical needfor automatic 
computing machines that could solve. complex dynamic 
problems. Slide Rule and manual equation solving could 
no longer be relied upon as a practical approach to 
seeking engineering solutions. Therefore, analog com- 
puters became important tools in the design of air- 
craft, jet engines, atomic reactors, oil refineries, 
chemical plants, etc. 


Many types of analog computers have evolvedover the 
years. The family has included the mechanical differ- 
ential analyzer, electromechanical differential ana- 
lyzer, and most recently, the iterative differential 
analyzer. 


The analog computer has always had several advantages 
compared to a digitalcomputer. These are primarily 
speed, more simulation capability per dollar, an ability 
to integrate, and an excellent man-machine interface. 





The mainfeature of the analog computer is that it cam 
integrate time-varying voltages. There is no easy 
way to differentiate. Consequently a mathematical 
model of aphysical system which is expressed in terms 
of differential equations cannot be solved with the 
machine directly. It is necessary to reformulate the 
mathematical model in terms of integral equations, 
either implicitly or explicitly. 


The analog computer can integrate only with respect 
to time. Thus, a mathematical model which contains 
partial integrals (corresponding to partial derivatives) 
with respectto several variables must be approximated 
by a set of ordinary integral equations with respect to 
time. Of course, computer time need not correspond 
to time inthe physical world, although it usually does. 


Before 1959 the analog computer was a synchronous 
machine; all its integrators operated in unison. In 
1959 the DYSTAC! was introduced. The name was an 
acronym for 'dynamic storage analog computer'. This 
machine was the forerunner of the iterative differential 


1 
DYSTAC is a registered trademark of CSI. 


Two large SD 80 computers are used in this engineering laboratory to stimulate the behavior of 
drone helicopters with different load configurattons and vartous automatte flight control sys- 


tems under a wide range of operating conditions. 


Simulation saved time and money, reduced the 


exposure of personnel and materiel to possible damage when actual flight tests were made, and 
itt also eliminated the necessity of watting for spectal envtronmental (weather) condittons. 


(Photo courtesy, Gyrodyne Company of America, Inc.) 
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Analog computers have long served tn the 
fields of chemistry and process control as 
a conventent, low cost means to observe, 
analyze, control, and predict the effect 

of varying parameters in a dynamic problem. 
Some of the problems may pertain to enzyme 
reaction, chemical kinetics, continuous 
distillation, heat transfer or transport 
delay -- just to mention a few baste applt- 
cations of analog computers. The university 
student shown here is performing a research 
problem involving the effect of potential 
barriers on kinette energy levels. 


(Photo courtesy, University of California) 


analyzer which appeared in1960. The iterative differ- 
ential analyzer is an asynchronous computer; the 
integrators need not be controlled in unison. They 
can operate independently either in groups or singly. 


In 1962 the analog computer was augmented with digital 


logic. This innovation first appeared in a machine 
called the HYDAC! The HYDAC had a very large 
quantity of synchronous digital logic. Since that time 
the inclusionof asmaller complementof asynchronous 
digital logic has become accepted practice. 


Digital logic can be used for the implementation of 
logical decisions. These are based on results obtained 
from the analog portion of the computer during the 
solutionof the problem. Digital logic also can be used 
for mode control of the analog computer. All of the 
Systron-Donner 10/20 and 40/80 series analog com- 
puters can be operated in the iterative mode and can 
be augmented with digital logic. 


Analog computers have found widespread acceptance 
in virtually every area of scientific investigation. This 
growing interest in analog computers has created a 
need for complete software, specially designed for the 
beginner and less experienced user. Also, the recent 
addition of digital logic control has greatly improved 
the problem-solving capability of analog computers. 
How this new feature can be used inanalog computation 
is thoroughly illustrated in this publication. 


It is the purpose of this handbook to provide students 
as well as experienced computer users with compre- 
hensive and up-to-date analog computer software. 
Chapters 2 to 5 develop the basic fundamentals of com- 
puter operation and illustrate the solution of elementary 
problems. Chapter 6 provides a detailed description 
of the operating controls and computer logic of the 
SD 10/20 and 40/80 series computers. This informa- 
tion serves as useful reference material to problem- 
solutions illustrated inthe more advanced discussions 


1 
HYDAC is a registered trademark of EAI. 
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An SD 80 computer mounted inside a Boeing 
experimental jet transport, selected for 
the NASA sponsored Supersonie Transport 
(SST) Program. In thts actual in-flight 
application, the SD 80 ts inserted between 
the pilot's controls and the aerodynamic 
control surfaces of the jet plane. This 
permits the total control system to assume 
the dynamics of any of a wide variety of 
SST types. 


(Photo courtesy, the Boeing Company) 





A student in mechanical engineering ts shown 
how to stmulate a mass-spring-damper system 


on an SD 33800 analog computer. Stating 
equations is unnecessary. Using block 
programming techniques, illustrated tin 
Chapter 4, the student need only understand 
the baste relattonshtp of physical vartables 
and constants. Following the program block 
schematic for a given system, the student can 
easily program the problem on the computer, 
observe results on the oscilloscope, and make 
further parameter adjustments to seek an 
opttmum solutton. 


which relate iterative programming techniques to the 
Systron-Donner computers. Chapters 7 to 22 develop 
more sophisticated programming techniques and appli- 
cations on a progressive basis. 


On the Apollo Program, Douglas Atreraft 
coupled an SD 40 computer to a large centrt- 
fuge. The computer calculated and integrated 
errors in human performance during Apollo 
lifting body reentry stmulatton studtes. The 
Systron-Donner computer ts seen next to the 
programmer, in upper left portton of ptcture. 


(Photo courtesy, Douglas Aircraft Co.) 








Analog computers are now standard computing equipment for classroom teaching and research work 
in colleges and universities. Students in the departments of Electrical Engineering, Mechan- 
teal Engineering, Chemistry and Biosciences receive instruction in the use of desk top analog 
computers as basic electronic model builders of dynamie problems. The ease and swiftness of 
presenting a solution on a readout (oscilloscope, XY recorder), and the ability to vary problem 
parameters and observe immediately their corresponding effects, have made the analog computer 

an important teaching aid. 


(Photo courtesy, University of Santa Clara) 


1-3/4 


CHAPTER 2 
BASIC ANALOG COMPUTING ELEMENTS 


This chapter shows how electronic equipment and 
circuits are used to implement mathematical rela- 
tionships in an analog computer. 


OHM'S LAW 


Ohm's law describes the relationship between the 
current through, and the voltage across a passive 
impedance. A passive impedance is a collectionof 
passive elements such as resistors and capacitors 
connected together in anarbitrary way. Such an impe- 
dance is generally considered (in analog computing) 
to be a two-terminal network which can be denotedby 


where Z is the dynamic impedance of the element. 
Ohm's law states 


E = ZI 


for 


where I is the time-varying current (in the direction 
indicated) through the passive element generated by 
the impressed time-varying voltage, E (with the 
polarity indicated). 


In order to simplify what follows, transform notation 
will be used where 


and 
ts tat 
: 


For a resistor Z = R and Ohm's law is 


E = RI 


R J 
—ww——— 
+ E - 
For a capacitor 
dE _ 
C at I 
so that Ohm's law is 
_ ll 
E = Gs I 
or E = ZI 
where Zo J 
~ Cg 
IT 
ey 
+ E - 


SERIES and PARALLEL IMPEDANCE 


Impedances are additive in series: 


21 =Z,+Zo 


This canbe proved in a simple way. The total voltage 
from point a to point b is Ep = Ey + Eg and by 
Ohm's law 


E, = Z,1 + ZI = (Z 


T 1 2 1 


When impedances are connected in parallel the total 


impedance can be found as the reciprocal of the sum 
of the reciprocals: 
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Jyt 
1 Zo 23 


N 
NM 
it 


N 
N 


T 


This again can be proved by simple application of 
Ohm's law: 


Lp = I + I, + I, 

1 2 3 
ie, ee eee ae a 
Zn E Z, Zo Ze 


The total impedance of two parallel impedances has 
the simple formula: 


oi, pete apt 
to 1 Z,+2Z, 
Z,+2 


iw) 


DIODE 
A diode is a non-linear resistor. Its resistance or 


impedance depends on the direction of the current flow- 
ing through it. It is denoted by 


=== 


+E- 
where 
E = ZI 
Z = Ry E>0 
= Ro, E <0 
and generally 
Ry Say 
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Typically, Ry, ranges between 1 ohm and 100 ohms; 
and Rg between 100, 000 ohms (100K®) and 1, 000, 000 
ohms (1MQ2). Thus the diode is an approximation toa 
switch for which 


R, = 0 


Ryo = # 


OPERATIONAL AMPLIFIER 


An operational amplifier is of the type that is called a 
d-c amplifier; it amplifies not only time-varying volt- 
ages, but also d-c or constant voltages. It is charac- 
terized by its excellent stability and extremely high 
low-frequency gain (amplification). It is denoted by 


where 


These voltages are measured with respect to ground 
(zero reference). The gain, A, is usually frequency 
dependent and will decrease with increasing frequency. 
As will be seen presently, this becomes a limiting 
factor in the use of a computer at high speeds. An 
equivalent circuit for the amplifier is 





where, generally, Rj is greater than 10° ohms and 
Ro less than 1074 ohms (closed-loop) at zero fre- 
quency. However, the amplifier is current limited, 
That is, it will only perform satisfactorily if the out- 
put current, Ih, is less than some value. The cur- 
rent limit for the S-D! amplifier is +25 milliamperes 
,at +100 volt output. The amplifier is also voltage 
limited; it will not function satisfactorily if the output 
voltage, €9, is greater, in absolute value, than some 
upper limit. The limit for theS-D solid-state 100 volt 
amplifier is 105 volts. Since the output of these amp- 
lifiers contains unwanted noise whose magnitude typi- 
cally can be 1072 volts, their effective useful rangeis 
about three and one-half decades (5x1072 to 102 volts). 


GENERATION OF TRANSFER FUNCTIONS 


In what follows, it will be assumed that the input 
impedance of the operational amplifier is infinite, the 
output impedance is zero, and the gain is infinite. 


1 Abbreviation for Systron-Donner. 


These assumptions introduce negligible error at zero 
frequency and are a good approximation at mid-fre- 
quencies. At high frequency, the assumptions cannot 
be made (particularly for gain). 


Consider an amplifier with input and feedback (from 
output to input) impedances and applied input voltages 
as shown. 








(e, -e.) (e, -e.) (e, -e.) 
Ip = a + = 
1 2 3 
©; ar 
I = 
Z 
: f 
tp Se Fay 


Now if the input impedance, R;, of the amplifier is 
assumed to be infinite 


°j 
Fee 
1 
for finite e, . Thus 
ip = de 
and 
(e, -e.) (e, -e.) (e, -e.) e. -e 
1 he IS a eee) 
Z, Zo Zs Ze 


Further, if the gain, A, of the amplifier isassumed 
to be infinite, then 


ea Ae 
oy > ek 0 


This is a reasonable assumption at low frequencies 
since as noted above €,, at most, will be in the neigh- 
borhood of 100 volts in absolute value. Thus, equation 
(1) becomes 


Me oe 
Z Z 


bo 
wo 
Hh 


— 
do 
oo 


The simplest case is 








ee fe scene 
oO Z. i 
1 
which can be expressed 
fo Ze 
= AG ae Samer eS ose s ne (3) 
1 i 


By a suitable choice of impedances many desired trans- 
fer functions can be generated. (Henceforth, the units 
of megohm and microfarad will be used for resistance 
and capacitance respectively. ) 


SUMMER (ADDER) 


Consider the configuration shown below, which is 
called a summer. 


Ry Re 

e| 
Ro 

eo & 
R3 

e3 


If, in equation (2), the substitutions 


1 1 
Ag Rs 
oa 
Ze = Ry 
are made, then 
° 1 Rp 3 








A set of typical values (in megohms) for these resistors 
' in a summer in a computer is 


R, = Rp = R; = 1.0 
Rg = 0.1 

so that 
e, = "ey “ee -10e, . 


In a computer, access is usually provided to the input 
(summing) junction so that additional input resistors 
or feedback components can be added to the summer 
circuit externally. Thus, the summer circuit is 





where the gains (multiplying factors) are indicated at 
each input, and the summing junction terminal by J. 
The program symbol for the summer is 


=> 


where the input gains are omittedif they are unity. If 
it is desired to indicate the junction or summing junc- 
tion (high-gain input), the symbol becomes 


INTEGRATOR 


The high-gain input is labeled only if pertinent. Con- 
sider the configuration shown below, which is called 
an integrator 


ej Te) 
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Substitution in equation (3) of 


Z, = R, Zs = as 
leads to the transfer function 
% 1 
e RCs 
or 
o> RCS “ - Ae ene 


If R=1, C=1, then 
t 


Gane = ff etoet 


O 


so that the output voltage, e., of the amplifier will be 
the integral with respect to‘time of the time-varying 
input voltage, ej 


If a switch is included as shown 


and if 


1) there is an initial charge stored on the capacitor 
which results in the voltage, €y (0) 


2) the switch closes at t=0, 


then 
t 
G1 76, (0) - fede. 


oO 
It is common practice to call this initial voltage the 
‘initial condition’ for the integrator. The derivation 
of this terminology is obvious from mathematics. 
Itis necessary to find a practical way to establish the 


initial condition or initial voltage for the integrator. 
To accomplish this, the circuit below is used. 


“Cre 


eo 


Here 
Le 1 2. RCs +1 
go ee Te R 
f 
R 
Ze RCs +1 
Z = R 


and substituting in equation (3) 


C2 = C2 )Ge ) eines 
oO RCs +1 IC RCs +1 “IC 


or 


de, 
RC rs + ey, = ec 


The solution of this differential equation is 


2 —t 
ey = Aexp (Ra) + eo 


where A is a constant depending on the initial voltage 
stored on the capacitor before Cra was applied. Then 


Lim 
ee £5 = ater 


For practical cases, it is only necessary that t>10RC, 
since e~10<, 0001, whichis compatible with the accu- 
racy of the circuit. The result in simple terms is: 
to guarantee the establishment of the initial condition 
it is necessary to wait 1ORC seconds after the voltage, 
“810 has been applied to the circuit. 


The two previous circuits, combined with appropriate 
switches as shown below, constitute a practical inte- 
grator 


0.1 0.1 
~€Ic 


The operation is as follows: 


With Sg closed and Sj open, eo = ejc after 
1.0 second (or 0.1 second if C = 0.1, etc.). 


With Sg open, eo still equals eyc, since the 
current through the capacitor is zero. (This 
must be true since both Sj andSg are open and the 
input impedance of the amplifier is assumed to be 
infinite. ) 


With Sy closed at, say, t=0, then 
t 
e (t) = €iq 7 f e, (t)dt. 
oO 


Ifatt = T, Sy is again opened, e, will stop chan- 
ging and remain at the last value before Ss, 
opened, namely 


T 
SCD) Se J e, (t)at. 


Oo 


If S2 is closed again while Si is open, the out- 
put of the integrator will return to 


eo = fc . 


Each integrator in a computer has these switches as 
part ofitscircuit. They are open or closed depending 
on what the programmer wants the integrator todo. 
The state of these switches is called the integrator 
mode. If all the integrators are controlled in unison, 
the switch states are determined by the main computer 
mode. 


The modes have simple names with obvious interpre- 
tations. During the R, or reset (IC) mode, initial 
voltages are impressed on the integrator capacitors. 
During the compute or C mode, the integrators inte- 
grate input voltages. During the hold or H mode, the 
integrator outputs remain constant at the lastvalue 
achieved before entering the hold mode. Thus, the 
mode permits the programmer to stop the computation 
at any time, enabling him to evaluate what has happened 
thus far in the calculation. The table below shows the 
states of Sy andSg inthe various modes. Thenumeral 1 
indicates the switch is closed (logical 1) while 0 
indicates the switch is open (logical 0). 


MODE S 1 So 

R(IC) 0 i. 
H 0 0 
C 1 0 


Finally, the general integrator circuit is shown below 
together with a definition of its transfer function for 
various modes. 


“ete 
1 Lo 
e| 
1 Lo 
£2 
io Ol 
es 
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MODE TRANSFER FUNCTION 
R @y = &0 
Cc ett) = €10 
t 
- J {ex + ey (t) + 10¢3(0] dt 
oO 
H e,(T) = & 
T 
- See + eg (t) + 1040) dt 
oO 


Note that the gain (multiplying factor) of the computing 
inputs is noted next to the terminals and. is determined 
by the value of the input resistor. The SJ terminal is 
provided to allow the other external input resistors to 
be added. The IJ terminal provides the ability to gen- 
erate an initial condition which is the sum of several 
voltages. The J terminal permits the connectionof 
additional feedback elements around the amplifier (to 
be discussed in chapter 10). Also, these terminals 
provide an external connection for the use of solid- 
state switches in place of the mechanical switches 
(relays) S; and Sg. In this case S; and So are con- 
strained to be open regardless of integrator mode and 
the necessary switching is done by the external solid- 
state switches. 


The program symbol for the integrator is 


The input gains are usually omitted if they are 
unity. 


POTENTIOMETERS 


It is necessary to have a device for entering constant 
parameters inthe computer program. This is accomp- 
lished with a potentiometer. The circuit for this ele- 
ment is shown below. 


ej 


eo 


If the output voltage, e,, of the potentiometer (pot) is 
applied to an input of another element of the computer, 
then the input impedance, RL, of the other element is 
connected from the pot output to ground. 
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Then 
I, = 1 + ro 
so that 
HY 20h Ohne ot 
Ry Ry Ro 
and 
A = e il + de + gh 
Ry o R, R, Ry 
Since 
R, + Ry = R, 
where R, = total potentiometer resistance, 
Ro = R- Ry 
and 
e = 1 e 
BP A ag es a 
Ry Ry, Ry 
a R, R 
= = 2 
Oo Ri Ry + Ry Ry 
e =ae 
oO i 


Now, obviously, it would be time-consuming to deter- 
mine Rj, with the knowledge of Rt, Ry, for eachnew a. 
Consequently, in practice, the pot is set with the load 
connected, by reading the output voltage, e , with a 
meter, for aknown input voltage, which is ustially 100 
volts. 


Thus the meter reading is equal to 100 aand the pot is 
changed until the desired value for a is obtained. The 
symbol for the potentiometer is 


ae 


where His the notation for the high end, or input. The 
His usually omitted since it is obvious from the com- 
puter program which side is the input to the pot. The 
element discussed above is called a two-terminal pot 
since it has two available terminals (or connections) 
on the computer program board (patchboard). Some- 
times it is desirable to connect the bottom (low) end 
of the pot to some computing-element instead of ground. 
In this case the low end is made available at the patch- 
board. The circuit is 


and the program symbol is 


or 


where notation for the arm, A, is omitted when it is 
obvious from the computer program. 


ARBITRARY FUNCTION GENERATOR 


Many mathematical problems to be solved with a com- 
puter require the generation of an arbitrary function. 
This is accomplished with a device called an arbitrary 
function generator (or frequently, diode function gen- 
erator, since the internal circuitry uses diodes). This 
device allows the programmer to approximate the de- 
sired function with straight linesegments. Anexample 
is shown below 





Here the function f(e,)is approximated by f‘(e,) with 
three line segments. In general, each function gen- 
erator, depending on howitis used, will provide either 
10 or 11 line-segments with maximum slope changes 
of 20or2.5:1. Severalof these devices can be used to- 
gether if more segments are required. The location 
of the slope discontinuity is called the breakpoint and 
is adjustable. Detailed instructions ‘for the setup of 
this element appear in Chapter 6. 


The principle of operation of the function generator 
depends on the use of a non-linear input impedance 
for an operational amplifier. That is, the impedance 
generates an input current proportional to the function 
which in turnconstrains the output voltage to have this 
functional relationship to the input voltage. A simpli- 
fied circuit is shown below. 





“I, = I, 
1 = - f (e,) 
e, = f’(e,) 


The program symbol for an arbitrary function gener- 


ator is 


FIXED FUNCTION GENERATOR 


Fixed function generators are used to generate often- 
used functions such as sine, cosine, log. etc. These 
are similar to arbitrary function generators in opera- 
tion, but do not permit the programmer to change the 
parameters within the device. Generally, a fixed 
function generator is more accurate and has better 
frequency and noise specifications than an arbitrary 
function generator. The program symbol for a sine 
generator, for example, is 


x(t) Sin x(t) 


MULTIPLIER 
The multiplier used with the S-D computer is calleda 
quarter-square multiplier. The name derives from the 
equation 

XY = E: + yy - ax)? | 


The actual equation to be implemented with hardware 
in order to provide multiplication is 


@) &) = 70 {l¢1 + e, |” - | ee e, || (4) 


The multiplier module contains two fixed function gen- 
erators. Ablock diagram for one of the generators is 
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| 

Xx 
eo E 

7 I 
-e F 

Cc 

N 
-e5 1 


(Note that both polarities are required as inputs for both 
ey and Eo: ) The transfer function is 


le, + e,|” 
ao 


where K is aconstant that determines the feedback 
resistor of the output amplifier (to be shown below). 


The block diagram for the other is 


¥ | F 
| 
Xx 
~@o E 
D I 
-e, F 
Cc 
N 
eo 2 
The transfer function is 
ley ~ €g : 
T= - K—~ 790 


The multiplier module block diagram is 
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The transfer function is 
_ «K 2 Ol 23 et 
1 = dalle, * el -|e, - e9|?] = «53 


A simplified block diagram for the module is 


e; M R 
U 
-e, L 
T 
— 
eo M I 
O 
a D 


If the module is connected to an output amplifier 





“ey; &9 





e 
0 100 


The program symbol for the multiplier module (with- 
out the output amplifier) is 


WZ 
a 


-e, 


With the designations 


M: The module is used as a multiplier 
I: The output is a current 
X, Y: Bipolar inputs 


If the bipolar inputs are not naturally available from the Note that the ''minus sign" associated with the ''M" indi- 


program (e.g. other variables) then summers can be cates a negative current is generated for a positive 
used to generatethem. This can be done in two ways: product. 
Case I: The program symbol for the multiplier module together 


with an output amplifier is 


ex Co 


Cx Cy 
100 





As before, summers can be used to generate the bi- 
polar inputs if they are not otherwise available. 


CIRCUIT 


Case I: 


ey 


PROGRAM SYMBOL 





CIRCUIT 


Note that the "plus sign" associated with the ''M" indic- 
ates apositive current is generated for a positive pro- 
duct, e,ey. The sign is generally omitted since this 
is the normal condition. 


Case II: 
ey C9 
eye 
e =4+ —4°¥ 
Cy : 100 


PROGRAM SYMBOL 





(Note that here the sign inversion is due to the output 
amplifier. ) 


CIRCUIT 


Case II: 





PROGRAM SYMBOL CIRCUIT 
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PROGRAM SYMBOL 


A multiplier module canbe used to make a divider with 
the same output amplifier. (This connection can be 
made conveniently at the patchboard. ) 








Here 
ene e, 
I= Kig 7 - 
so that for R = i 
K 
ee 
-~. 27 
x" 100 
or 
100e 
x 
e = - 
oO e 
y 


The programmer's symbol is 


ex €o 
lOO e, 
ey "ey 





From the circuit it can be seen that the condition ey>0 
must hold in order to avoid instability. This is easy 
to show. The multiplier module can be considered to 
be a non-linear feedback resister, R,. The sign of 
R, is the same as the polarity of e y Suppose in the 
circuit 
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R¢<0. Then, If¢ will be in the direction shown, for 
€9 > 0, and I, will produce ane; <0. Since the gain 
of the amplifier is negative and large, the circuit will 
be unstable. An additional restriction, |e,| <|ey|, is 
necessary to ensure that the output amplifier not ex- 
ceed the voltage limit. 


The divider can be used to generate square root in a 
simple way. Ife. =e 


y 0” 
e, Cy 
CIRCUIT 
then 
ex 

oO 

or 
e 2 = -100e 

x 

so that 


= 10 fe e_<0. 
© x? °x% 


The diode insures that the system will not saturate in 
the wrong direction if e, goes positive inadvertently. 


This circuit has the program symbol 


ex =) e@, lOex. 


FUNCTION SWITCH 


A function switch is a manual switch which can be oper- 
ated from the controlconsole. It is used to change the 
computer program during execution. The switch has 
three positions: up, down, andcenter (off). Its circuit 
is shown below. 


,o 
— 


This is the program symbol as well. The terminals 
of the switch are located at the patchboard. 


FUNCTION RELAY 


A function relay is used to make program changes auto- 
matically during execution. It is energized (logical 1) 
for an input equal to or greater than +28 volts and is 
de-energized (logical 0) for an input equal to or less 
than 0 volts. The circuit is shown below 


Ux 

n ~~! 

ee ee 
U U,=X if U=0O 
U,=X if U= 1 


This is also the program symbol where N is the relay 
number. 


ELECTRONIC SWITCH 


An electronic switch is used in much the same way as 
afunction relay; the only restriction being that it must 
be inseries with the junction of an amplifier. Its speed 
of operation is much greater (10-5 sec) than afunction 
relay (10-3 sec). It requires the same input voltages 
(logic levels) as the function relay (i.e. logical 1: 
28<e, <100, logical 0: -100<e,_ <0). 


The program symbol is 


a Ug 


Ug=a if U=| 


whereN is the switch number. 


COMPARATOR 


Frequently, itis required to determine the sign of the 
sum of two variables. 


A circuit which accomplishes this is 





If (ce, + eg) > 0, then eg is limited to approximately 
O volts by Dj because of its low resistance when it con- 
ducts (logical 1). If (ey + eg) < 0, then eg is limited 
to approximately +28 volts by Do (logical 0). The 
latter is true since 


1) Ip mustbe in the direction shown (e, + €,>0) 


and eS > ey: 


2) €9 = 28 volts since ex must be approximately 
0 volts due to the high gain of the amplifier. 


The program symbol for the comparator is 


. U 
b 

U=! if (a+b)<O 
U=0 if (at+b)>0 


If both logical comparator outputs are required, then 
€, is complemented with abiased analog inverter. The 
program symbol is 


and the analog circuit is 


C2 


-|00 
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CHAPTER 3 
ELEMENTARY ANALOG PROGRAMMING 


All dependent variables in an analog computer are vol- 
tages. Consequently, it is necessary to equate these 
voltages to physical variables. Since the elements of 
the computer are voltage limited, scaling will be re- 
quired. The art of scaling is discussed inChapter 5. 
In this chapter scaling is ignored since the necessary 
concepts can be developed without reference to scaling. 


As pointed out in Chapter 1, it is necessary to remem- 
ber that an analog computer cannot differentiate easily. 
It is a machine whose main feature is its ability to 
integrate. Consequently, problems which are defined 
by differential equations must be reformulated in terms 
of integral equations either explicitly or implicitly. 
Consider the problem 


y = y(x) 
dy x 
a By = 0 

y(0) = A. 


This can be rewritten as the integral equation 


For this problem computer time will represent x, and 


y(t) + Bf yat = 0, t = x 
y(0) = A 


An integrator, as shown in Chapter 2, integrates from 
t = 0 onward. It also requires an initial condition: 
namely the value of the integral at t = 0, which in 
this case is A. Thus, the problem is formulated for 
the computer as 


t 

yi) = — f By(t)at, (1) 
0 

y(0) = A 


The right half of equation (1) can be generated by 


A/\00 
- 100 


t 
+y/(t) -f By (t) dt 
0 


where the constant B and the initial condition A are 
inserted by means of pots. Remember that the inte- 
grator inherently has asignchange. Equation (1) states 
that y(t) is equal to the negative integral. All that is 
needed, to complete the computer program for this 
problem, is to connect the output of the integrator to 
the input of the pot set to B. Hence, the computer 
program is 


A/|00 
-l00 


y (t), t=x 


The independent variable, y, can be recorded from the 
output of the integrator. 


In this example the differential equation was reformu- 
lated explicitly in terms of an integral equation. It is 
frequently possible to do this in an implicit way, as 
the next example will show. 


Consider the problem 


y + Ay + By = f(t), y = y (t) (2) 
y(0) = C 
y(0) = D 


This is equivalent to 
y = f(t) - Ay - By 
y(0) 
y(0) = D 


Cc 


First, y can be generated from y by 
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(Note that again the pots are used to generate initial 
conditions.) Usingthis circuit, the complete program 
can be generated by 


D/100 C/100 
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where equation (2) is satisfied by applying the correct 
inputs to integrator 1. 


Next, the use of non-linear elements is demonstrated. 
Consider 


y + Ay + ty = cost,1<A <10, 
y(0) = C, 
y(0) = D. 


The computer program (except for scaling) is 





As another example, consider 


y + f(t)y = 0 
y(0) = A 
y(0) = B 


The computer program is 








The operation is as follows: The clock pulses are not 
applied to DCU 1 until the machine is in compute be- 
cause the output of the RH"' bus is a logical 0 in I C 
and this is an input to gate 1. The FF is initially in 
the R state when the machine is in compute having been 
reset by the FRT' Logic output. Divider 1 reduces the 
clock frequency of 1000 cycles/see to 100 cycles/sec. 
The output of the 0.count of counter 4 is 10 cycles/sec 
and is connected to the input of counter 5. The 7 count 
output of counter 4 and the 1 count output of counter 5 
are connected to gate 2. Thus, whenever t = 0.17 sec 
the output of gate 2 will change from Oto1, which will 
change the FF tothe Sstate. At this time the computer 
will go to the hold mode due to the S output of the FF 
being connected to the problem hold terminal at the 
patchboard. Also, the clock andcounters will be reset 
because of the logical i applied to the OV terminal. 
As soon as the machine goes to the reset mode the 
FRT"' logic output will cause the FF to return to the 
R state. Thus, the logic circuit is ready to be used 


The next example shows how the DCU's (Chapter 8) can 
be used to determine the solution of a differential equa- 
tion for a specified value of the independent variable. 
Suppose it is required to find y (.17) for 


y+ry=0 
y(0) = 100. 


From above, the analog program which generates y is 


-l00 


In order to determine y(.17) the computer is allowed 
to compute for 0.17 seconds and then put into hold. 
The latter can be accomplished by applying +28 volts 
to the "problem hold'' terminal at the patchboard after 
0.17 seconds. This will halt the entire computer and 
in particular the computation carried out by the above 
program. The digital program which will do this is 
shown below: 


PROBLEM 
HOLD 
TERMINAL 


FRT’ 


again as the computer is manually put first in the reset 
mode and then the compute mode. 


ALGEBRAIC EQUATIONS 


The first examples were intended to illustrate the basic 
approach to programming. The following illustrates 
a more serious application of the analog computer. 
Suppose it is desired to solve a set of simultaneous 
algebraic equations which are expressed in matrix form 
by 


AX =C (3) 


where the n x n matrix A and the column matrix C 
are known. For the purpose of illustration it will be 
assumed that A is of rank 2, although the derivation 
which follows is perfectly general. Consider 


K+ AX =C (3a) 
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When the system has reached steady-state (assuming 
there is one) 


X = 0 
BK ee: 


Thus the solution to equation (3) is obtained from the 
steady-state solution of equation (3a). The latter can 
be written 

= 


%y F Ay yy > Ayo%o 1 


Ko + Ag1Xy + AgoX. = Cy (4) 


and the analog program is 





Generally, integrator capacitors are chosen equal to 
0.001 yfd to decrease the solution time. However, 
there is no guarantee that the system of differential 


equations, (4), is stable. A different approach will 
provide an unconditionally stable system of differential 
equations. If 


AX =C 
it follows that 
A'AX = A'C (5) 
where A' is the transpose of A. The solution of equa- 
tion (5) is the same as that for equation (3). Again at 
steady-state the solution of 
X + A'AX = AIC (6) 
will be the solution of equation (5) and therefore equa- 
tion (3). The stability of equation (6) is determined by 
X + A'AX = 0. (7) 


If equation (7) is multiplied by the row matrix X, then 
the result is 


XX + XA'AX = 0. 


It is well known that XA'AX is positive definite (i.e. 
non-negative for all values of X). Unless X = 0, 


XxX < 0 


and stability is guaranteed. The computer program for 
equation (6) is 





In this program it is easy to change individual coeffi- 
cients because each pot is associated with only one co- 
efficient. The program can be simplified so that it 
requires no more equipment than the previous pro- 
gram if the following coefficient changes are made. 


OLD NEW 
a a +a 2 
11 11 21 
A192 441219 * 291%99 


Ao4 A11242 * 494%99 
2 2 
Ao9 ao * AyD 


This simplified program saves equipment but is in- 
convenient to use when the solution is required for 
several different values of the coefficients. 


ARBITRARY CLOSED FUNCTIONS 
Many analytic or closedfunctions can be generated by 
representing them in terms of their generating dif- 


ferential equations. Thefollowing examples illustrate 
the technique. 


f(t) = 1/t, t>a: 


This function has the generating equation 


df 2 
a- 7 
1 

f(a) = > 


The program which solves the generating equation is 





f(t) = A(t + a)"\*t >0: 


df n 





‘dt tea 
(0) = Ad? 


Generating Equation 





Program 


f(t) = In (t+ a, t >0: 





dé __1 
dt t+a 
{(0) = Ina 


Generating Equation 





Program 








f(t) = A cosh(t), t>0: 
f-=f 
f(0) = A 


Generating Equation 


A 


100 
-100 


A cosh(t) 


Program 


Mean value: 


The mean value, ¢ (t), of a function, f (t), over the 


interval [Ty T, + t] is defined by 


T+ 
g(t) = = \ Se f(u)dy, t>0. (8) 
Ty 


The generating equation for ¢ can be obtained by dif- 
ferentiation: 


Se = Str, +1) - o(t)} (9) 


The program is 





Suppose the integrators 1, 2are in the initial condition 
mode until t = T 1, at whichtime their mode is changed 
to compute. Then, theoretically, the program will 
generate the solution of equation (8). However, the 
program willnot work because the right side of equa- 
tion (9)is indeterminate fort = 0. Thus the divider, 
in the program will have an unstable output. Practi- 
cal limitations require the Y-input to be at least 3 
volts for useful operation. Consequently, it is neces- 
sary to choose a different initial condition: 


Tite 
o(e) =+ Jf teas, 
Ty 


This program is 





As before, when t' = T,, the mode of integrators 1, 
2 is changed from initial condition to compute. Initially 
S1 is closed and S2 is open. The output of integrator 
1 will not be ¢ (t) untilt = «. Whent = «, the cor- 
rect initial condition, namely ¢(«), will have been es- 
tablished for integrator 1. Atthis time the comparator 
will be activated, which will open S1 and close S2. 
Also the output of integrator 2 will be «. Thus, the 
generating equation is implemented by the program 
starting when t =«’. The output of integrator 1 will 
be ¢ (t), t>«. 


COORDINATE TRANSFORMATION 


Suppose two coordinate systems have the same origin 
and are displaced by an angular rotation, 6: 





An X-coordinate point, say x, will have the coordinates 


x’ = xX cos 6 
y’ = -x sin 6 


in the x', y'-coordinate system. Similarly, for a 


y-cvuordinate, say y, 


x’ = y sin 0 


y’ y Cos @. 


Thus, a vector, R, in the X, Y-system with components 
X, y, will have the components 


i 


x! x cos 6 + ysin 6 


-xX sin 0 + y cos 8 


y’ 
in the X’, Y’-system. In matrix notation 
R’ = TR 
where R, R’ are the representations of the vector in 


the unprimed and primed coordinate systems respec-~ 
tively, and where 


T = | cos 4 sin 6 
-sin 6 COs 6}. 


In the three-dimensional coordinate system 


positive rotations 6, ¢, \ respectively about the X, Y, 
Z axes correspond to the transformations 


1 0 0 
0 cos @ sin 6 
Q-sin 6 cos 6 


aj 
ul 


cos ¢ 0 -sin ¢ 
Ry = [0 1. 0 
? sing 0 cos ¢ 


cos A sina 0 
R, = |-sinA cosa 0 
0 0 1 


Any arbitrary three- dimensional rotation can be repre- 
sented asa product of these three matrices. The super- 
scripts refer to the axis about which the coordinate 
system is rotated and the subscripts to the name of 
the angle of rotation. 


The transformation of a vector, V, by R* is 


Vv’ = RV 

i = Vy 

Vy = Vy cos 6+ V_, sin 6 
Vy = “Vy sind + V, COs @. 


The program is 





Note that 6 is limited by the range of the sine and cosine function generators. 
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The program symbol used frequently is 





The transformations R), R” have similar programs. 


MATRICES 
A matrix, A, is a collection of elements: 


A fag} tbo 


(in what follows it is assumed that n = m). Thus, if 
n = 3 


a 


11 12 743 
A =|851 ag gg 
431 439 433 


Multiplication by a constant, a, follows the rule 


11 712 09713] | 11 12 13 
Wo, gq Woe 
W431 gq ge 


R 
> 
HI 

R 

© 

R 

> 
| 


ie = a...) = a.. 
ij =e sai 
when it is understood that the a are the elements of A. 


Addition follows the rule 


A+B=a.. + = (a + b).. 
ij ij 

or for, say, n = 2 

a a b b a,,+b a,.otb 
poe al At 1a lene ia Pee at “aie Pye 

291 boo} |Po1 Poo} {891 +Pe1 Aaa tPoa|, 
Multiplication: 

AB ~ 2 Aig? aj 

Thus for, say, n = 2 

44917 + AyoPo1 = 2yyPy2 + AyoPoo 


b + a,,b 


Big + 8gQPo4 — QzPyQ + BgqPoo!, 


01°11 
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Note that AB # BA. The unit matrix, I, is 


be ya ed 


= 0,1 # j 
Thus for, say, n = 3 
1 0 O 
I=/|0 1 0 
0 oO 1 
I, has the property 
IA = AI =A 


Differentiation: 


or for, say, n = 2 





2 : da 1 das 

aa _d i1 12 | _ dx dx 
dx  dx|a a = 

21 22 dao, das 

dx dx 


Transpose: The transpose, A’, of A is 


A’ =. 
ji. 


or for, say, n = 2 


Submatrices: The submatrix, Sij, corresponding to 
the element. ajj in A is obtained by 
deleting the i- th row and j-th column 
of A. Thus for, say, n= 3 


et f 11 “12 713 21 Aap 
13 ~| 221 #22 93 
a a 431 439 
931 339 


Determinant: The determinant, lA|, of A is a scalar 
quantity and can be found by induction. 
By definition 


[a| .  ( 1g. i 83; ‘|, for any i. (10) 
(Note that some particular choice is 


made for i). 


Forn = 1 


Thus, forn = 2, 


el: 


J 


S - S 


1\| Saye ag Ate 


2 
: 12 


(-1)' a, .|s 
1 J 


Ayy%o9 ~ 449494 


or 
2 og 
[al - ae a9 ilSoj| = -421821+@225o0 
Ao 12497499944 


and either calculation leads to the same result. For 
n = 3, 


3 bos 
Al= ¥ (-1)'%a., [S..], 
| | ae ) ij | ij| 
and in particular, any of the equations 


S 1+j 
4] = Zen 215181)| 


s 2+j 
oe a> i[Soj| 


[4| 


I Meo 


- 
JA[ = 2 G1)°"ag,|s. 


j=1 


will lead to the same |A|. Ifi = 1 


[Al = 24181172498) 9413513 
= 841 @p9%35-A93%39) ~ 24 9(@p19337293431) 
+ 843(8912g97A90%31) 
~ 941799733791 17939307 732791%33 + 442293431 
* 4943491939741 3%90%31 
Alternatively, A| can be determined from 





n oe 
|A| ee S (-1)'a,|s.J, for any j. (11) 
i=l Joa 
A determinant can be multiplied by a constant, a, in 
a simple way: 
|vA| = alAl. 


From equations (10) and (11) it can be seen that if a 
matrix B can be generated by m adjacent row inter- 
changes and p adjacent column interchanges starting 
with matrix A, then 


|B] = -1)™"P fal. 


Cofactor: Thecofactor, Ci corresponding tothe ele- 
ment ay in A is J 


or for, say, n = 3 








a a 
ec 4 — is 
31 32 
and 
a a 
coe s ie ie 
31 32 |. 








‘ -1 : . 
Inverse: The inverse, A ~, of matrix, A, is 


-1 


Ci 
A = 

[A] 
For example, ifn = 2 


4444997449791 


Orthogonality: A matrix, A, is said to be orthogonal 
if 


Note that R™ RY, RZ are orthogonal 


matrices. The geometric interpreta- 
tion of orthogonality is rotation. 


Logical Matrix Operations: 


(AB)' = BrA’ 


(AB) i eae 


DECOMPOSITION OF ORTHOGONAL MATRICES 


Frequently the programmer is presented with a com- 
posite orthogonal (rotation) matrix which is a function 
of several angular variables. Usually these angular 
variables represent successive rotations about coordi- 
nate axes. Thesimplest way to generate the computer 
program is to decompose the matrix into a product of 
the RX, RY, RZ matrices and then generate subpro- 
grams for each ofthese matrices. This can be accom- 
plished by following a simple set of rules. These are 
stated without proof and examples are given. R“, RY, 
R42 are said to be type 1, 2,3 matrices respectively. 
Suppose the composite matrix, M, can be represented 


M = A, Ag. Ay - 
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The outer matrices are A,; Ay and the inner matrix 
Mp is defined by 


M, = AyAg.-- Avy. 


Now suppose the simplest element (lowest degree homo- 
geneity) of M is located in rowi and columnj. Then A 

is of type i and Ay is of type j. The angular variables 
corresponding to Aj, Ay will appear in the products 
of columns j, irespectively. The sense of the rotation, 
Aj, (positive or negative) can be found by setting all 
the variables corresponding to Ag..., An equal to 
zero. The sense of the rotation, Ay, is similarly 
found by setting the variables corresponding to A, ..., 
Ay-1 equal to zero. If the simplest element of M is 
zero, then M has the representation 


M = A, Ay 


and there is no inner matrix. If the simplest element 
of Mis not zero, then there is an inner matrix and My 
can be found by setting the variables corresponding to 
Aj, Ay equal to zero. The further decomposition of 
M] is carried out with the same procedure as for M. 


Example: 


Cos ¥ Sin ¥ 0 
M =[-Cos @ Sin ¥ Cos 6 Cos ¥ Sin 06 
Sin 9 Sin v -Sin 9 Cos v Cos @ 


The simplest element is zero and is in row 1 and col- 
umn 3. Therefore there is no inner matrix and Aq, 
Ag, are of types 1, 3 respectively. The products of 
column 3 contain the variable 6 and the products of 
row 1 contain the variable for ¥. Therefore, the var- 
iable for Aj is 6 and the variable for Ag is ¥. If ¥ = 0 


1 0 0 
M = {0 Cos 6 -Sin 6 
0 ~Sin 6 Cos 6 


If @ = 0 
Cos ¥ Sin¥ 0 
M = |-Sin ¥ Cos ¥ 0 
0 0 1 
Finally, 
XZ 
M = R,Ry 
Example: 


Cos¢ Cos¥ + Siné Sing Sinv 
M = -Cosé@ Sin ¢ 
Cos¢ Sin¥ + Sing Sing Cosy 


Sing Cos¥ - Siné Cos¢ Sinv -Cosé@ Sin¥ 
Cosé@ Cos¢ _ sing 
Sing Sin¥ - Sin@é Cos¢d Cosy Cosé Cos¥ 
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The simplest element of Mis Sin ¢@ and is in row 2 and 
column 3. The products of column 3 contain the vari- 
able ¥ and those of row 2 the variable ¢. Therefore 
A,, Ag are of types 2, 3 with variables v, $, respec- 
tively. 


If 6 


@ = QO, then 


Cos ¥ 0O -Sin wv 
A, = 0 1 0 
Sin wv 0 Cos ¥ 


Cos @¢ Sin g¢ 0 


Ag = -Sin ¢ Cos¢ 0 
0 0 1 


If ¢ = © = 0, then 


1 0 0 
Ay = M, = 0 Cose¢ Sin 6 
0 -Sin @ Cos @ 
Finally, 
Y X .Z 
M = Ry R, Ry 


As can be seen from the previous example, it is not 
necessary to find all the elements of A, to determine 
the sense of ¥. Since it is known that Aj is of type 2 
it is only required to evaluate the element in row 1 and 
column 3 with 9 = ¢ = 0. Similarly, if the element 
in row 1 and column 2 is found for ¥ = 6 = 0, then 
the sense of ¢ is determined. 


Example: (Here we use the abbreviations:C = Cosa, 
S_ = Sina) = 
a 
© gl yr€ 8,5, S548 8 ,C, 
M = C §,C.-C C,S.S_-S CC 


o 0° “6° OXY gS O-d 
Te eee ae 


Pa” pov #88 8,5, -C,8,C, C,S,-8,8,C, 
S55 9C g- 8,648, S540 gC ar S870 ,8,C, 


~8,C,8,-C 8, CC, 


The simplest element of.M is C,C, which is in row 3 
andcolumn 3. The products of column 3 contain 9 and 
those of row 3 contain ¢. Thus, Ay is type 3 and Ay 
is also type 3. 9, ¢ are the variables for Aj, Ay 
respectively. Consider the element in row 1 and 
column 2. For \ = ¥ = @ = 0 its value is Sing and 
consequently 


Todetermine the inner matrix, set 9 = ¢ = 0 so that 


Cy 0 Sy 
M, = “857 Cy SCY 
C8, “8. CC, 


The simplest element of Myis a zero in row land col- 
umn 2. Also for ¥ = 0 the element in row 2 and col- 
umn 3 is Sind. The products of column 2 and row 1 
contain the variables), ¥%, respectively. Thus 

X LY 


M, = R, Ry - 





Finally, 


(ee Doe WZ, 
A=R g Ry RUR a 
COORDINATE TRANSFORMATIONS WITH UNRE- 
STRICTED ANGLES 
The coordinate transformation program is the same as 
before, except for the use of Sine and Cosine function 


generators. For unrestricted range we use the vari- 
able 6 and the equations 


Sino = f6Cose dt 
Cose =-féSine dt 


The program is, 


-l00 


100 
cos @ 


The lower two multipliers and summer generate 


E = -(100- Sin26 - Cos“9). 


This error is used to change the damping of the inte- 
grators to ensure orthonormality of Sing, Cosé. 
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CHAPTER 4 
BLOCK PROGRAMMING FOR PHYSICAL SYSTEMS 


A Valuable Aid to Simulation of Lumped-Parameter Physical Systems 


The analog computer, or differential analyzer, is a 
powerful tool for obtaining dynamic solutions to differ- 
ential equations. Generally the equations to be sol- 
ved relate to a physical system of particular interest 
to the computer programmer, in which case the pro- 
grammer wishes to simulate his physical system by 
means of an electrical network whose defining equa- 
tions are analogous (hence the name analog computer) 
to the subject physical system. It is desirable then, 
to arrange the computer mechanizationto have a one- 
to-one correspondence between the problem board and 
the physical system, so that the programmer "sees" 
his system rather than an abstract electrical network. 


The following discussion will illustrate a method of 
programming in which the physical elements are 
treated as blocks rather than as sets of differential 
equations. Certain useful basic building blocks will 
be developed, and methods for combining the basic 
blocks to form any desired physical system (simula- 
tion) will be discussed. Simulations generated by the 
methods to be described will have the advantage that 
any or all physical parameters appear as single con- 
trols inthe computer so that any parameter may be 
changed without affecting any other setting. In addi- 
tion, any or all variables appear as computer outputs 
and may be individually monitored. 


A lumped-parameter system, by definition, consists 
of a number of interconnected discrete elements whose 
individual or collective responses to impressed forces 
or stimuli are of analytical interest. The forces and 
responses are related by a mathematical operator. 
Block programming starts with the selection of basic 
definition and rules. 


MATHEMATICAL OPERATOR (0) 

The operator, or combination of operators, describes 
the functional relationship between quantities. The 
most common operators are the following: 


Summer 


Constant Multiplier (potentiometer, commonly 
referred to as "pot"’) 


Inverter - -1 Sign Changer 
1 
Integrator - -> (La Place Transform Notation) 
Differentiator - s (This operator is never used in 
an analog simulation if it can 
be avoided, as it usually can). 
Variable Multiplier 
Arbitrary Function - f(x) 


Analytical Function - Trigonometric, Log, 
Hyperbolic, etc. 


Inverse Operator - (0)? S, + are inverses 


ELEMENT 


The element is the basic unit of the block representation 
of a system. It hasa pair ofterminals or nodes, with 
a value associated with eachnode. It alsohasa trans- 
ference quantity between nodes, or "through" the ele- 
ment whichis functionally related by some operator to 
the node-pair value. The transference quantity may 
also be referred to as the branchtransference or trans- 
mission. 


The elements required to depict the passive lumped- 
parameters of most physicalsystems are very few in 
number. Inparticular, five basic elements will suffice 
for many linear electrical, mechanical, and thermal 
systems. These arethe Summer, Constant Multiplier, 
Sign Changer, Integrator, and Differentiator. 


BLOCK 


A block is the computer mechanization related toa 
physicalelement. The input/output variables of acom- 
puter block are voltages and currents, but the operator 
relationship is the same as the simulated element. 


Table I gives examples of single elements and their 
corresponding blocks. Table II gives afew commonly 
occurring two-element combinations. Combinations 
of more than two elements can be generated as shown 
in the following examples. Usefultables of more com- 
plicated blocks or transfer-function simulations appear 
in many publications. 


DRIVING SOURCE 


A driving source is a source of energy or power for a 
network of elements. 


SYSTEM STABILITY 


A stable system is one whose responses are bounded 
(finite amplitude limit) for any finite input. An un- 
stable system is one whose responses are not bounded. 


NETWORK 


A network is an interconnected set of elements and 
sources. Elements may be joined together at their 
nodes, providing the connected nodes have common 
dimensions and value, and provided the transference 
quantities have common dimensions. 


This discussion concerns elements and networks in 
which the followingtwo rules apply: (Kirchoff's laws) 


1. Thealgebraic sum of the node values around 
any closed path (loop) in a network is zero. 


2. Thealgebraic sum of the transmission quan- 
tities at any node is zero. 





; Korn and Huskey, "Computer Handbook", Chapter 2, 
McGraw-Hill, New York, 1962. 


4-1 





TABLE I. Single-element Computing Blocks 


PHYSICAL SCHEMATIC MATHEMATICAL . 
ELEMENT SYMBOL RELATIONSHIP PROGRAM BLOCK 


Summer + x _l>— 
Y = - (X, + X,) xX 


Constant 
Multiplier 
(POT) 
Inverter 
= -X b 


Y 
2 
x——({K)——Y 
x 
oo 











Differentiator 
Y=sxX 







X 






Y 

Y 
Y 
-Y 





Resistor 


Torsional 
Spring 










Heat 
Conductor 










Viscous 
Damper 






TABLE I (Continued) 


PHYSICAL SCHEMATIC MATHEMATICAL 
ELEMENT SYMBOL RELATIONSHIP PROGRAM BLOCK 


Electrical 
Capacitor 


Viscous 
Damper 
Mechanical 

Mass 


Thermal 
Mass 


Electrical 
Inductor 


Mechanical 
Mass 





COMPUTER PROGRAM MECHANIZATION a) assigning aposition inthe mechanization 
COMMENTS diagram to all necessary variables and 
operators. 


b) accomplishing the necessary intercon- 


1. Operational amplifiers operate with the junc- nections so that each output has its pro- 
tion at virtually zero potential; the currents per inputs. 
into the junction through all the input and 
feedback paths must sum to zero; and there c) checking to see that all input/output rela- 
is asign inversion between the amplifier in- tionships in the mechanization are satis- 
put and output quantities. fied. 

2. Anyclosed-loop ina problem mechanization EXAMPLES OF NETWORK PROGRAMMING 


simulating a stable physical system will gen- 


erally have an odd number of sign inversions. The following examples illustrate some of the tech- 


niques used in applying the block programming meth- 


3. In mechanizing an equation, set the highest ods. Examples are given for linear elements only. 
order derivative only on the left of the equals Non-linearities can be included simply by replacing 
sign to obviate the necessity for differen- the potentiometer corresponding to the non- linear ele- 
tiation. ment with a multiplier, function generator, or other 

suitable computing element. Backlash, hysteresis, 

4. The mechanization is accomplished by: stiction, etc can be included as required. 
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TABLE Il. Two-Element Computing Blocks 


MATHEMATICAL 
SCHEMATIC RELATIONSHIP PROGRAM BLOCK 
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MASS - SPRING - DAMPER SYSTEM 


X,,V Xo, V: 
rl B 2742 
F MOTION IN X 
K DIRECTION ONLY 
2 
OO OO 


Spring K,;: F5, = K, (xX, - 0) 





Element - force equations 


; ; 7 Initial values att = 0. 
Spring Ky: Fj, = Ky (X, - X,) 

X, (0), X, 0), v, ©), vy (0) 
Damper B: F, = B (vo -v_) 


D 1 


Forces acting on masses 


F = -F , + Fo,t+ F 
Ss Ss 


Ml 1 2 D 


Fa * “Fs2 - Fp t+ F 


Steps to mechanize 


1. Assume xy; -Xo, “Vy> Vo and F are available 


2. Draw blocks for M,> My 


3. Draw blocks for K,, K,, B 


V2’ 


4. Make input connections to mass blocks to produce xX, -Xo, “Vy> Vo 


| | -Xo | 





aaa aes eae ee ee 
bo 






| spring | massi loamPerR | sprinc2! mass 2 
fi as eee eet See ee ee 
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The above mechanization makes each parameter avail- 
able as an individual adjustment. The programming 
is easily accomplished by generating the forces im- 
posed on the elements and then operating on the forces 
to obtain velocity and displacement. Any or all of the 
initial conditions, X,(0), X9(0), vz(0) and vo(0), may 
be zero, or may also be variable. 


Note that the block representing the damper contains 


two summers. The first summer produces (vy - V9), 
and the second summer acts only as a signchanger. 
The two-amplifier, one-pot combination can be reduced 
to a more simple two-pot combination by connecting 
v4 and vo through B pots to the mass block inputs. 
The less complex alternate is shown below. This con- 
figuration requires fewer amplifiers, but variations 


in M,; Mo, or B necessitate two or three adjustments. 





* Do not connect to -v 


1 


to get pot setting of B instead 


of E. A coefficient pot output cannot be connected 


M, 


to another coefficient pot input. This is a practical 
limitation of the computer hardware. 


Rigid Bar Supported on Springs 





o. forg bd OE 
oe sI 
Fi = -K,Y 
F, = -K, Y, 


Ff, - F,4, + Fé 


Rigid bar of mass M, moment of inertia about C.G. of 
I. Assume 6smallsuchthat @~sin@. Vertical motion 
only. Yo is displacement of C.G. 


; : F,+F,+F 
0 sM 
Y = Yo - £49 





MASS 


SPRING 2 IROTATIONAL | 


SPRING IT 


One-Dimensional Heat Transfer 


K- THERMAL CONDUCTIVITY 


M-MASS 


Cc 


SPECIFIC HEAT 





je ele 
3°3 
| THERMAL 
MASS 





Sa ee 
| qT | 
To | 

| 

| 


THERMAL 
MASS 


MAL 
UCT- 
E 
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Electrical R-L-C Networks 





RESISTOR 





PAPAGITOR) REstisice eee eon Neen 


ro 


RESISTOR 


E 


C; 


oD 
pe 
oD 
— 
aw H 
oD 
eX 
m 
a | sr 
N t 
we a 
tl i] 
J ma ie 
' 
—= 4 a 
re re 
al Ni] 
IS HQ jn 
ul tt 
ul om co” 
1 
rf 
| 
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TWIN-T NETWORK 











se, as a ae dl 
9 Ri . 3 SC 
aes as Es 
3 R, 
-(I, +I.) 
g * Is. 
(Eq - E,) = sC, 
I, > I, 


Ey = (I, - I, +15 +1) Rg 


E, = (E, - E,) + Ey; E, = E4- (E, - Eo) 





| 
a | 
Sy 


~y-< 
2|- 
= 


oN 


o 
eee 


-(I,-IptIz+Ig) 


C)R3 
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CHAPTER 5 
SCALING 


The analog computer is voltage limited to a practical 
range of operation of +100 volts as a maximum and +1 
volt or zero as a minimum. It is usually not possible 
tonumerically equate analog voltages to physical vari- 
ables. That is, the magnitude of physical variables 
must be scaled to fall within the useful range of the 
machine. This procedure is called amplitude scaling. 


The analog computer is also limited as tothe speed with 
which it will solve a problem. Practical bounds on 
solution time are 100 seconds as a maximum and 100 
milliseconds as a minimum. Events in the physical 
world usually occur in time intervals which fall out- 
side these limits. Thus, the simulation of real world 
phenomena with the analog computer is ordinarily 
faster or slower. The procedure for relating computer 
time to physical time is called time scaling. 


Correct scaling is animportant factor in reducing sim- 
ulation errors. An error analysis of any particular 
program is about as complex as the problem which is 
being solved by the program. Consequently, except 
for unusual cases, an analysis of error is not made. 
An estimate of the accuracy of the program can be 
found from check cases for whichthe answer is already 
known. Error is reduced by scaling the program so 
that all potentiometer values and amplifier gains are 
reasonable. Frequently, in the process of scaling a 
problem, it will be found that some parts of the model 
are not significant and can be eliminated. Scaling is 
a good check on reasonability. 


Scaling is anart, not a science. A good deal of exper- 
ience is required to become proficient. Scaling canbe 
done inseveral ways. Thetechniques discussed in this 
chapter are those in most prevalent use among analog 
programmers. 


If X is a physical variable corresponding to computer 
voltage, V, then a scale factor, a, is chosen so that 


Ve=a-xX 


will fall within the practical range of the computer. It 
is good practice to make V as large as possible. If 

is physical (real) time corresponding tot,, the com- 
puting time, then these two can be related by a scale 
factor, N: 


ie m  ,mam 
te = Nt, dt, = N ot 


It follows that 


frat, = Nf (Foat,. (1) 


N is chosen so that the problem solution time on the 
analog computer falls within reasonable limits. Prac- 
tically, this is accomplished by a choice of N which 
will not result in prohibitive amplifier gains. 


The general procedure for scaling is: 


1. Generate anunscaled program having a math- 
ematical structure which agrees with the pro- 
blem (model) to be solved. 


2. Identify the location of all voltages corres- 
ponding to physical variables. 


3. Associate a scale factor with each of these 
voltages (e.¢. V,. = aX). 


4, Label the program with the scaled physical 
variables (i.e. a, X) rather than the voltages. 


5. Estimate maximum values of the physical 
variables. 


6. Choose the amplitude scale factors andthe 


time scale factor, N, so that all pot settings 
and gains are reasonable. 


In order to accomplish step (4) it is necessary to know 
how input and output scale factors are related for 
various computer elements. 
POTENTIOMETER 

Physical equation: Y = bX 


Scale factors : Vv. = aX, Vi= ae 


y 
“y 
Scaled equation : Vy = ae Vig 
Program ax ( ) a, ¥ 
ba 
a 
x 


Note that ascale change can be made withapot. Assum- 
ing either a, or ay has been established previously, 
the other is choseh so that ba,/ay is a reasonable pot 
setting. The pot can be usedsolely fora scale change. 
In that case b = 1. 


SUMMER 

Physical 

equation: Z = X + Y 

Scale 

factors : V_=a_Z, V_=a X, V_ =a Y 
Zz Zz x x y y 

Scaled a, a, 

equation: Vv, a Vy + ay Vy 
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Now, if only one-gain inputs for the summer are used, 
then it must be true that a= a =a and the pro- 
gram is y 


A,X 


a7Z 
—ayyY 


Thus, it is seen that input scale factors must be the 
same for a one-gain summer and that a change of 
scale cannot be made through asummer. However, 
adifferent input scale factor can be used fora10-gain 
input. The program for this case is, 


-ayX 
—avlY 10 


where it must be true that a. = 10a. 


INTEGRATOR 


As with the summer, the input scale factors (except 
for the 10-gain input) must be the same. However, a 
change of scale can be made with an integrator. Both 
amplitude and time scaling can be done. Scaling for 
an integrator is derived below, where C is the value 
of the capacitor in ufd, and R is in megohms. 


Physical 

equation: Y =f Xat , 

Scale 

factors : V_= a Y, V_=a.X,t_ = Nt 
y y x ok c p 

Scaled a, 1. 

equation: vy" a, Z NAVxtt, 


Thus, the integrator gain, ae must be 
aR 
RC a (2) 
and the program is 
— xX NayY 
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MULTIPLIER 





Physical 
equation: XY = Z 
Scale 
factors: V=akxX, V_=a_Y, V_=aZ 
x x y y Z Z 
Scaled 
equation: V_V 100a 
x y Zo ow V 
100 aa Zz 
xy 
Program: 
ayX a,Z 
ayY 
It must be true that 
100a,, 
= 1 
aa 
xy 


From this requirement itis seen that the relationship 
among the scale factors is the same as that for input and 
output multiplier voltages: 





DIVIDER 


Physical X 
equation: YY" Z 


Scale 
factors : V_= aX, V_=a_Y, V_ =a Z 

y y Z 
Scaled Ns aa 
equation: 1005— - T00a = Vv. 

y x 
Scale 
factor a, 
constraint: a_ = 100-— 
Zz a 

y 

Program: 
a,X a,Z 


SQUARE ROOT 





Physical 

equation: Y = /X 

Scale 

f : = = 

actors : vy aX, Me ax 

Scaled a 

equation: Vy= y x 10/V, 
10 /a,, 

Scale 

factor 


constraint: = 
in 4, 10 |/ a. 


Program: 


kaon 


FUNCTION GENERATOR 


Physical 
equation: Y = f(x) 


Scale 
factors : Vi=ayY, V. = aX 


Scaled 


equation: Vi=af ae 
y7 afl) 


EXAMPLE; 


Mass-spring-damper system 





Physical 
equations: Spring Force, ES = KX 
dX 
D = po 
amper Force, Fy B at 
p 
Mass 2 -(F_+F,) 
Acceleration, GEM 8 
2 M 
dt 
p 
‘ -4 
M=1;B=4x10 


dx 
at, (0) = 0, X (0) = .95 


>K=6.4x10 


dt 


Estimates: (=| =8x 100: Xmax. = l 
p/max 


The unscaled program is: 





The scaled program is obtained by assigning literal 
scale factors, determining their numerical values, and 
calculating the coefficient potentiometer settings. 





The = and factors appear 
N2 N 
due to the inherent gain of N 


through each integrator, from 
Eq. (1). 


In general, the determination of the numerical values 
involves atrial and error approach. Any convenient 
parameter can be picked for the starting point. Ifa 
selected value leads to unreasonable potentiometer 
settings or unreasonable amplifier gains, then new 
values may be required. Reasonable pot gains are 
between .1 andi. Reasonable amplifier gains are 
between 1 and 20. Amplifier gains should be integer 
values. Non-integer values can be obtained from a 
potentiometer, in combination with gain-of-10 inputs, 
if necessary. 


1. Determination of a. In order to utilize the 
full dynamic range of integrator 2, set 
aX max = 100 (volts). 


a = 100 
5-3 
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Potentiometer 1 setting. 


= .95 (3) 





Determination of oe At this time, a trial 


value for N maybe selected. Analog solution 
times of .1 to 100 seconds, or computer nat- 
ural frequencies from .1 to 100 radians per 
second are the preferable operating ranges. 
Anestimate of the physical time or frequency 
ranges can be used to select a value for N. 


For this example, a value of N = .01 appears 
reasonable. 


b 100 (volts) (4) 


N (dX/dt,,)max. 


Looking ahead, note that integrator 2 will 
require a gain of ax . The gain of N is in- 
herent in the integrator from Eq. 1. There- 
fore, the zm value for the integrator will be 
7. It is highly desirable that this factor be 
an integer, preferably either 1 or 10. There- 
fore, b = 100 will be used. The effect of this 
is to force integrator 1 to utilize something 
less than its fulldynamic range, but the effect 
is not serious inthis example. As a general 
rule, it is not possible to scale so that all 
amplifiers work over their full range. 





Potentiometer 2 setting; 


p2 = x = (6.4x10-%)c (5) 


N aM 





The value of c is chosento assure that a rea- 
sonable value for P2 is obtained and to assure 


that 2 is an integer. 


c = 100 will satisfy both conditions. 


P2 = .64 
Potentiometer 3 setting: 


cB _ 04 


P3 = NpM ~ 


This setting for potentiometer 3 is below the 
desired minimum potentiometer setting value 
of .1. Inthis example, however, it is not pos- 
sible to have a large setting for potentiometer 3. 
This follows from a consideration of the phy- 
sical problem. The damper, B, was chosen 
tohavea very slight effect on the mass-spring 
system, so it necessarily follows that the 
portion of the analog program related to the 
damper willhave onlya slight effect. It should 
be noted that abnormally low potentiometer 
settings will frequently arise in simulations of 
systems containing negligibly small elements. 
It should be recognized that the low settings 
stem from the actual system characteristics. 
The programmer, therefore, need not spend 
time in futile attempts to improve the scaling. 


The final program with numerical values is: 


GENERAL COMMENTS ON TIME SCALING: 


The following points concerning time scaling are of 
particular importance: 


1. 


If all integrator gains are simultaneously 
changed bya factor k, the time scale changes 


from N to and all computer frequencies 


k > 
increase by k. The amplitude scaling, how- 
ever, does not change. All numerical scale 
factors, all potentiometer settings, and all 
computer problem voltages remain the same. 
This property of amplitude invariance under 
time scale change allows the programmer to 
easily speed up or slow down the computer 
solution. 


The most convenient way of changing all inte- 
grator gains simultaneously is by changing 
the capacitors by a fixed factor. Integrator 
capacitors can be simultaneously changed by 
means of the computer Time Scale control*. 


An event which occurs in a fixed real-time 
interval can occur in several different com- 
puter time intervals depending on the choice 
of integrator capacitors, that is, depending 
on the time scale factor. 


Events of identical character which occur in 
different real-time intervals can all be made 
to occur in the same computer time interval 
by appropriate time scaling. 


1 On the SD 10/20 and SD 40/80 analog computers, a master Time Scale switch is a standard feature. This per- 
mits fast and convenient switching of capacitors that are patched -- up to 1000:1 time scale (x10, x100, x 1000, 
depending on patch panel connections). 
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CHAPTER 6 
COMPUTER OPERATION 


This chapter is based on materials from the Systron-Donner 10/20 and 
40/80 Computer Operating Manuals. The descriptions of the computer 
mode controls, computer logic, and various computing modules permit 
the reader torelate all applications and programming examples given in 
this handbook to an existing computer. For example, Chapter 8 offers 
detailed information on the use of digital logic. The modules needed to 
implement the applications described in Chapter 8 are shown in this 


chapter. 


Not all available S-D computing components are described here, only 
those which illustrate and serve to explain a basic function. For the 
latest, up-to-date availability of all computing elements and modules in 
the SD 10/20 and 40/80 series, please consult Systron-Donner or your 
nearest SD engineering sales representative. 
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SD 10/20 Computer 


The Model SD 10/20 is an all solid state analog com- 
puter with an operating range of +100V and expansion 
capabilities for up to 24 operational amplifiers. The 
principle features of the computer are full +100V com- 
puting range, all solid-state design, compact modular 
construction of all computing components, removable 
problem board which mates directly with computing 
modules to eliminate all problem board cabling, a 
control center with readout selectors for amplifiers 
and potentiometers, high-speed reed relays and solid- 
state switches, complete logic capability for repetitive 
and iterative operations, 0.01% accuracy of computing 
resistor and capacitor networks, a panel voltmeter 
having a null mode of measurement with 0.02% (of full 
scale) resolution. The amplifiers are protected against 
accidental grounding of their outputs as are all power 





supplies. 
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All computing components in the Model SD 10/20 are 
modular plug-in units. .They are identical to those 
used in the larger SD 80. This high degree of standard- 
ization between SD 10/20 and 40/80 computers maxi- 
mizes. usefulness of computing modules. Since all 
modules are standard plug-in units, the computer user 
is offered a wide choice of computing components. 
As new modules become available, the SD 10/20-—40/80 
series computers can thus be constantly up-graded to 
include the new modules. Furthermore, a user may 
start with an SD 10/20 and later use the same 10/20 
modules in an SD 80. This common design feature 
results in important cost savings. As the need for 
more computer modules increases, the expansion cost 


from an SD 10/20 to an SD 80 is minimized. 


+ C&C Cees 
+ CC CCK OOS 
CRCCCCCS 





Panel Meter Ranges: 1, 3, 10, 30, 100, and 300 volts, and 
+ null. Full scale accuracy: 3%. Null position provides 
0.02% F.S. resolution with reference potentiometer having a 
+0.05% linearity at 25°C. 


Compute Time Selector — compute time continuously vari- 
able from 5 msec to 10 sec. Reset time varies from 5 msec 
to 5 sec, depending upon coarse steps of compute range. 


Meter Range Selector — with positions for 300 v, 100 v, 30 v, 
10 v, 3 v, and 1 v. Serves also as sensitivity adjustment for 
+ null. 


Function Selector — for rapid choice of: + null, — null, 
Meter, External (connects selected bus to external jack). 


Address Selector — address of all amplifiers and potentio- 
meters. 


Null Reference Potentiometer — provides high accuracy 
readout using null method with 0.02% F.S. resolution. Lin- 
earity is +0.05% at 25°C. 


Function Switches — provide manual switching flexibility 
in problem solutions. 


Mode Selection (lighted pushbuttons): 


Hold — places problem solution on all integrators into 
hold position. 

Compute — applies problem voltages to all integrators. 

Reset — applies initial condition voltages to integrators. 

Rep-Op — places integrators into a repetitive operation 
cycle. Compute time variable from 5 msec to 10 sec. 

Bal/Pot Set — disconnects junction and grounds the 
input resistor summing junctions of all amplifiers. 
Each amplifier is converted to a gain of 2500 for 
precision monitoring of junction offset. 

Pwer On/Off — energizes and de-energizes computer. 
Oven — indicates + 28-volt oven power is on to maintain 
constant temperature of computing capacitors. 

OL Hold, OL Reset — lights up when any amplifier is 
overloaded. When depressed, computer goes into 
Hold; when released, normal operation is resumed. 

Time Scale — activates relays in each integrator module 
to change computing capacitor. (x 10, x 100, x 1000, 
depending on patchpanel connections.) 

Slave — permits operation of computer control circuitry 
from a second console. 


Coefficient Potentiometers — up to 24, available in groups 
of 6, featuring 10-turn wire-wound with lockable counting 
dials. 


SD 40/80 Computer 


The same compact design and high performance com- 
puting capability found in the small SD 10/20 computer 
is also available in the large SD 80 computer. With an 
expansion capability to 84 amplifiers (126 special 
order), 125 potentiometers, and a full set of digital 
logic control modules, the SD 80 provides a problem- 
solving power capable of handling large, complex 
problems normally associated with the simulation of 
dynamic electrical, mechanical, thermal and similar 
systems. 


Yet, the large-capacity SD 80 fully retains the 
versatility of a desk-top computer. Mounted on a desk 
or table equipped with lockable casters, it can be 
easily moved to new locations. In addition to creating 
a highly compact and rugged computer, the all solid 
state design increases reliability, simplifies mainten- 
ance, and eliminates the need to operate the computer 
in an air conditioned laboratory environment. 


Programming a problem on a fully expanded SD 80 
is virtually as simple as operating a small SD 10/20. 
Visual computer circuits in symbols everyone under- 
stands show at a glance how a circuit is patched. 
This feature saves time, reduces patching errors, and 
makes patching a task that the non-expert can perform. 


To allow maximum operator control over a computer 
with a capacity as large as that of the SD 80, all com- 
puter controls are centralized in the left wing which 


is movable to any convenient position to suit the 
operator. Major SD 80 controls that speed and simplify 
problem checking and solving include: 


¢ A 199-point digital address pushbutton selector 
system for rapid access to amplifier outputs and 
potentiometer arms. 


¢ Static and Dynamic checks by simple pushbutton 
control to verify patching and test time scaling 
rapidly. 


¢ Instant alarm controls, both visual and audible, 
with automatic transfer to HOLD at the occurrence 
of an overload. 


¢ Individual control of compute and reset times for 
precise adjustment in Iterative and REP-OP modes. 


¢ Time scaling the problem with a single switch on 
the Control Wing. 


¢ One pushbutton to check any amplifier’s balance 
without repatching and to set any potentiometer. 


The problem board is removable, just as it is with 
the SD 10/20. There is no restriction as to grouping 
of the modular patch panels in the overall problem 
board layout. In both the SD 10/20 and 40/80 compu- 
ters, digital logic control is provided by modular plug- 
in units that slip into the central patchbay. They thus 
form an integral part of the computer and the main 
problem board, not a separate added-on element. 





1. Digital Voltmeter — 5 digit readout, automatic ranging and 
polarity. 


2. Address Selector Panel — with an amplifier and potentio- 
meter address capability of 199 points. 


3. Control Panel — for selection of: 

Op. Modes: Hold, Compute, Reset and Rep-Op. 

Master Time Scale Switch. 

Compute Time: 5 msec to 5 sec. 

Reset Time: 5 msec to 5 sec. 

Slave: for remote control. 

Oven and Power pushbuttons. 

Check Mode switches for Static and Dynamic checks. 

Overload Controls for: 1) audio alarm, 2) transfer to hold 
mode, 3) removal from Hold mode. 

Amplifier Balance and Pot Set: changes all amplifiers to 
gain of 2500 for simplified monitoring of junction off- 
set; and for proper loading of potentiometers. 


4. Function Switches 
5. Coefficient Potentiometers, No. 1-5 


6. Null Meter Reference Potentiometer 


7. Removable Patchboard. — An SD 80 board can hold up to 
42 individual computing module panels. All modules 
are mounted in a universal patchbay directly behind 
the removable problem board. Nomenclature on patch 
panels matches textbook circuit diagrams. Operator 
can thus ‘‘see’’ the circuit while patching a problem. 


8. Built-in Digital Logic Modules, consisting of Flip-flop, 
Gate, and Time/Event Control Modules. 


9. Potentiometer Wing can hold up to 6 panels for a total of 
125 potentiometers (first 5 pots are on Control Wing.) 
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Computer Logic 





Logic Buss 
F. R. (R. 0.) 


Each operating mode switch energizes a logic buss 
which is then distributed to the logic relays in the com- 
puting modules (primarily to the integrators). The logic 
levels in the Analog sections of the computer are: 


+28V = 
OV = 


logical 1] 
logical 0 


True = 
False = 


R.R. (R. 0.) 
A logical one energizes a buss and places the com- 

puter in that mode of operation. All contradictory modes 

have their switches interlocked to prevent two modes of 

operation from being energized at the same time. 


The table that follows outlines the logic busses and 
operating modes. F.R. 


In most cases the nomenclature for each logic buss 
is self-explanatory. The following explanations should 
clarify the remainder of the logic nomenclature. 


Table 1. CONTROL LOGIC TABLE 


COMPUTER OPERATING MODE 


LOGIC 
US ees PROBLEM O.L. STATIC | DYNAMIC | TIME SCALE 

F.R.(R.O.) ae +28 2 
PeRe. ee ee 
oe eae ++ —|—— 
Se pot 2 ft 

Pu ous {4 te poet 
Rn Bes __} 8 70 ee 
serie SS ee ee 


Time Scale 


(1) Blank indicates either 1 or +28 condition can exist, two values indicate cyclic operation. 


(2) Manual Hold is actuated by pushbutton switch on Logic Control Panel. 


(3) Problem Hold is actuated by application of +28V or greater Problem Hold Trunkline on Problem Board... 


(Note: Will not be damaged by any xoltage in computer range) 


(4) O.L. Hold is actuated by operating ‘°O.L. Hold’® switch with an O.L. condition (either temporary or permanent) or any amplifier in 


the computer. 
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(R. T.) Forward Reset (Real Time). 


+ 
NR 
oc 


Description 


Forward Reset (Repetitive Operation). 


This is the normal Reset buss and also 
is the Reset buss energized in the reset 
portion of Repetitive operation. 


Reverse Reset (Repetitive Operation). 


This buss is the logical complement of 
forward reset. It is used for complemen- 
tory integrators and track/store circuits 
in iterative computation. Energized in 
the Compute mode of operation and in 
the Compute portion of Rep-Op. 


This buss 
is used to reset non-repetitive operation 
integrators. It is energized only when 
the reset switch is engaged. 


Forward Hold is the normal Hold buss. 


It is energized in the Hold mode of 
operation or if an overload occurs when © 
the Overload Hold switch is engaged. 


Reverse Hold. Reverse Hold is the logi- 


cal complement of F.R.R.T. It isused for | 
complementary integrators and Track/ | 
Store circuits in iterative computation. 
Energized inthe compute mode of opera- 
tion and in Rep-Op. 





Use of the Mode Controls 


The mode controls are of the push button type with 
interlocks to prevent two contradictory modes from being 
energized at the same time. 


All control in the computer is accomplished in the 
integrators. Thus, when a mode is energized it opens 
or closes a relay in each integrator. Reference to the 
illustration, Integrator Block Diagram, on following page 
(6-6), while reading the following explanations will 
assist in understanding the control system. 


Reset 


The Reset or Initial Condition (IC) mode enables 
the integrators to be set to their initial values before 
starting computation. When the Reset relay is energized, 
an input and a feedbackresistor are associated with the 
integrator amplifier making it, in effect, an inverter. 
Applying a voltage to the input resistor {in this case the 
IC terminal) will cause that input toappear at the output 
and to simultaneously charge the integrator capacitor 
to the same voltage. 

Compute 


In the compute mode there are no energized relays. 
The logic is such that if the computer is not in Hold and 
not in Reset, it is in compute. 


Hold 


In the Hold mode a normally closed relay contact is 
operted which disconnects the input. signal from the 
amplifier. As there is no discharge path for the inte- 
grator capacitors the output will remain at its last com- 
puted value. In this manner outputs may be read or re- 
corded and the problem resumed at will by returning to 
Compute. 


Rep-Op 


In the Rep-Op (Repetitive Operation) mode the com- 
puter is automatically switched between the Reset and 
Compute modes of operation. This mode enables the user 
to obtain repetitive solutions which are easily monitored 
on an oscilloscope. 


The amount of time the computer stays in the compute 
mode is determined by the setting of the Compute Time 
switch located at the lower left of the control center. 
The outer dial of the switch may be positioned in any of 
three ranges: .005 to .1] sec, .05 to 1 sec, and .5 to 10 
sec or in the External (EXT) position. The inner dial 
adjusts the compute time over the range selected. 


The compute and reset times may be adjusted arbi- 
trarily by observing the solution to the problem and ad- 
justing the times until the desired display is obtained. 
If more precise compute and reset times are required the 
Compute-Reset cycle may be observed on an oscillo- 
scope. Connect the oscilloscope to any of the FOR- 
WARD RESET terminals on the integrators or to Test 
Point 6 at the rear of the computer. A rectangular wave- 
form will be observed. The positive (+28V) portion is 
Reset, the negative (OV) portion Compute. The Reset 


and Compute time controls may then be adjusted to give 
the desired times. 


If the Compute Time switch is placed in the External 
(EXT) position external circuits may be used to generate 
the Rep-Op cycle. The drive signals must be compatible 
with the logic used in the computer (+28V = True, OV = 
False) and should be connected to the R (reset) and 
C (Compute) terminals located on the problem board. 


OL Hold — OL Reset 

When engaged, the Overload Hold portion of this switch 
will cause the computer to automatically switch to the 
Hold mode at the occurrence of an overload in any. ampli- 
fier. When the overload is found and eliminated, the 
switch is re-engaged toreset the overload hold circuitry. 

If an overload occurs, the light within the OL Hold 
/OL Reset will begin to blink on and off and will con- 
tinue to do so until the overload is removed. 


Time Scale 

The Time Scale switch determines which of the 
integrator capacitors connected to the A and B terminals 
on the integrator modules is in the circuit. With the Time 
Scale switch not engaged, the capacitor connected to 
the A terminal is in the circuit. When the Time Scale 
switch is engaged the capacitor connected to the B ter- 
minal is in the circuit. Any two of the three available 
capacitors may be patched to the A and B terminals. If 
time scaling is not required the capacitors may be patch- 
ed directly to the output. 


Special Trunklines 


There are six special trunklines located on the pro- 
blem board of the Computer. These are located on the 
top left row of modules. The first module contains the 
M, X and Y special trunk lines and the third module 
contains the H, C and R special trunklines. The special 
trunklines are used as follows: 


M: This line is connected to the panel meter or the EXT 
jack as determined by the position of the Meter Select 
Switch. If the Meter Select Switch is in the METER 
function the M line is conriected to the EXT jack; if 
the Meter Select Switch is in the EXT function, the 
M line is connected to the panel meter. 


X &Y: The X and Y lines are connected directly to 
the X and Y terminals on the left side panel of the 
computer. These lines are intended for use with an 
X-Y plotter but may be used to connect signals to 
any external device or they may be used to connect 
external inputs into the computer. 


H, C&R: The H (Hold), C (Compute), and R (Reset) 
lines are used as external input lines to the Control — 
Logic circuitry. A logical 1 applied to any of the 
lines will place ‘the computer in that mode of opera- 
tion. To use these lines the computer must be in 
REP-OP and the Compute Time switch must be in 
EXT. 
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Integrator Modules 





INPUT GAIN f RESET © 


POT. 





MERON MEH 


Integrator Block diagram illustrating flexibility 
of circuit and patching convenience. 


Model 3320 


gg) Dual integrator with two uncom- 


mitted operational amplifiers. 


Model 3320, Dual Integrator 
The Model 3320 is divided into three sections: the 


top and bottom integrators and a logic section common 
to both. Basic patching of the integrator requires only 
a patch from the end of the integrating capacitor to the 
A terminal of the time scale relay. By patching the end 
of another capacitor to the B terminal a second time 


scale is set up and is controllable by the TIME SCALE 


switch on the control panel. 


Initial conditions (IC) are applied at the IC terminal. 
The source of the IC voltage may be a pot, amplifier 
output, the reference voltage, or from an external source. 


The logic section consists of a 7-terminal matrix 
near the center of the patching block. The nomenclature 
is the same as that explained in the Logic Control sec- 
tion of this manual. The two center terminals of the matrix 
are connected to the coils of the Reset and Hold relays. 
For normal and repetitive operation, patch from the F 
(forward) terminals to the relay coils. For reverse and 
iterative operation, patch from the R (reverse) terminals 
to the relay coils. For nonrepetitive operation, patch 
from the F(RT) terminal to the reset relay coil. (In this 
mode the integrators will reset ONLY when the reset 
button on the control panel is pushed.) Note that the 
logic patching controls both integrators in this module. 





To use the integrator as a summer, patch from the 
input side of one of the resistors to the output of the 
integrator. This will connect a resistive feedback around 
the amplifier; the other resistors may be used as summer 


inputs. 
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Model 3329 

Quad integrator, dual multiplier, 
operational relay. Contains four 
integrating networks, two four 
quadrant multipliers and one 
operational relay. Also available, 


less multipliers, Model 3329A. 





Model 3329 — Quad Integrator, 
Dual Multiplier, Function Relay 


The Model 3329 has four integrating networks without 
operational amplifiers. The operational amplifiers neces- 
sary for integration are normally patched into the Model 
3329 from an adjacent Model 3325. The input resistor 
network of the Model 3325 then becomes the input to 
the integrator. Relay logic connections are in the lower 
left hand section of the Model 3329 patchboard. Nomen- 
clature is the same as that used on the Model 3320, Dual 
Integrator. Each Integrator Network has a Hold (H) and 
Reset (R) terminal associated with it, these terminals 
must be individually patched to the Relay Logic termi- 
nals. The Squid (multiple) patch cords should be used 
when more than one integrating network is to be used 
at the same time. 


The Dual Multiplier section of the Model 3329 may 
be patched for division, squaring and square root func- 
tions as well as for multiplication. For multiplication 
or division operation, three operational amplifiers are 
required. In the division and square root functions arbi- 
trary limitations as to polarity and relative amplitude of 
the inputs are imposed. For example, in division the 
Y input must always be positive and equal to or greater 
than the X input. Note that in square root an external 
diode must be patched into the circuit. 


Patching and operation of the Function Relay in the 
Model 3329 is identical to that of the Function Relays 
in the Model 3322A. The relay is energized by a nominal 
+28 volt signal and will not be damaged by a signal of 
up to +100 volts. 





Summer Modules 


Model 3321: Dual summer with 
two uncommitted operational 
amplifiers. 


Model 332] — Dual Summer 


The Model 3321 contains two operational amplifiers and 
the components necessary for summing operations. Basic 
patching consists of connecting one end of the feedback 
resistor to the output. The top summer has three input 
resistors, the bottom summer two input resistors. There 
are four spare summing resistors located in the center 
section of the module which may be used as either input 
or feedback resistors with either of the summers. 


The module also has two 20 kohm_ resistors and 
two uncommitted diodes associated with each summer. 
These may be used to patch the summer as a limiter- 
comparator or for any of the common diode circuits used 
in computation. 


A function switch termination is provided on each 
summer module. The function switch is a single pole- 
double throw type switch and is physically located near 
the control center. Normally only the first three summing 
modules will have switches terminated in them. 


The center section of the summing module has ter- 
mination for a diode function generator. Use of these 
terminations and of the diode function generator is ex- 
plained in that section of this manual. 


Model 3325: Quad summer, four 
summing networks, with four 
uncommitted operational 
amplifiers. 





Model 3325 — Quad Summer 


The Model 3325, Quad Summer is identical in per- 
formance to the Model 3321, Dual Summer. As illustrated, 
patching of the Model 3325 is, however, considerably 
different than that of the Model 3321. Basic patching for 
any of the four summers in the module consists of patch- 
ing the SJ’ terminal to the J terminal and patching the 
right end of the feedback resistor to any of the output 
terminals. Any of the resistor input networks may be 
used with any amplifier. Thus summers with six, nine 
or twelve inputs may be patched. 


The Model 3325 also contains four spare resistors, 
four diodes and four trunk lines. In addition, it has 
terminations for a Function Switch and for a Diode 
Function Generator. 


The Model 3325is also used as a source of amplifiers 
for the Model 3329. When used in this manner, the Model 
3325 furnishes the amplifier and input resistors for the 
integrator circuit. 
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Model 3322A — Dual Inverter, Dual Function Rel ay 


The Dual Inverter contains two operational amplifiers 
for use as inverters or highgain amplifiers in comparator 
and limiter circuits. Operation as an inverter requires 
two patches: (1) From the summing junction to the junc- 


¥_ 100 «100 4 ? tion. (2) From the input side of one of the resistors to 
the amplifier output. The inverter may be used as a 
Ps summer by patching additional resistors into the sum- 


om ming junction. (Note: Potentiometers used as inputs 
a to the Model 3322A should be set in the Reset mode of 


operation.) 


The Model 3322A also contains two function relays. 
Each is a two form C (2FC) relay corresponding to a 
double pole-double throw (DPDT) switch. The relays may 
be energized by any voltage from +28 volts to +100 volts. 
The energizing source is typically the output of an amp- 
lifier or a logic signal, but any source within the +28 
volt to +100 volt range will energize the relays. Any 
voltage in the range plus and minus 100 may be applied 
without damage. 


As in the Dual Summer, the Dual Inverter contains 
uncommitted diodes, limiter resistors, terminations for 
a Function Switch and Diode Function Generator ter- 
minations. The use and operation of these elements is 
the same as that described in the section on the Model 


3321. 





Model 3322A: Dual inverter, dual operational relay. Two 
uncommitted operational amplifiers with 


two operational relays. 


ne ae 


a ee 
*e 


FUNC. SW. FUNC. RELAY 





Block Diagram, Model 3322A 


Model 3323 — Dual Inverter, Dual Multiplier 





Model 3323: Dual multiplier, two four quadrant quarter 
square type multipliers with two uncom- 
mitted operational amplifiers. 


Multiplication Operation. The multiplier/divider is 
based upon the quarter-square principle, in which the 
following relationship holds: 


V4 [ (X+¥)2 - (X-Y)2] = XY 


In +100 volt computers, the multiplier output voltage is 
always scaled to -0.01 XY. 


The circuitry consists mainly of two diode squaring 
networks, one based upon positive polarities and the 
other upon negative. For four-quadrant operation, volt- 
ages representing X, ~X, Y, and ~Y must be present at 
the inputs to the squaring networks regardless of the 
polarity of input voltages. 


The center section of the module contains the termi- 
nations for the multipliers. The section is divided into 
two identical halves, one for each multiplier. The multi- 
plier will also perform division, squaring, and square 
root operations 





CONNECTIONS FOR MULTIPLICATION 


Division Operation. Division is accomplished by the 
terminal connections shown in the figure below. The 
concept is illustrated by the circuit diagram in the same 
illustration. If the output of amplifier 1 is designated Z, 
the output of the multiplier must be+.01YZ. By reference 
to the basic discussion on operational amplifiers (para- 
graph 2.8.), it follows that X = -.01YZ, or Z = -100X/Y. 
Since +.01YZ and Z must always be opposite in polarity 
to X, it also follows that Y must always be positive in 
order to maintain this relationship. 


If Y occurs in the problem only as a negative voltage, 
apply it to the -Y multiplier terminal and patch the input 
and output terminals of the inverting amplifier to the -Y 
and Y terminals, respectively. Z still equals -100X/Y. 
The absolute value of Y must always be larger than X; 
otherwise, the output would tend to be larger than 100 
volts, causing an overload. 


O Y (ALWAYS +) 


INSERT 


CHANGE ONLY THESE 
CONNECTIONS TO CON. 
VERT FROM MULTI- 
PLICATION TO DIVISION 





CIRCUIT FOR DIVISION OPERATION WITH INSERT SHOWING 
HOW TO CONVERT FROM MULTIPLIER PATCHING. 


Squaring. Squaring is performed as a special case 
of multiplication where X = Y. Only one input inverter 
amplifier is required. Patch together the R and Y input 
terminals and the -R and -Y terminals. Leave the other 
connections as for multiplication. The output voltage 
taken at the amplifier output terminal represents -0.01X 2 


Squareroot Operation. The squareroot operation is a 
special case of division where Z = Y. The relationship 
Z = -100X/Y of division becomes Z2 = -100X for -100 
<X<0, Z=10VX , for0<X <+100, Z=-10VW%X. 


lf X is always negative, patch together the R and Y 
terminals and the -R and -Y. If X is always positive, 
patch R to -Y and -R to Y. Leave all other connections 
as for division. 

The Inverter sections of the Model 3323 are similar to 
those described for the Model 3322A. The two basic 
patches required for inverter operation are: From the 
summing junction to the junction and from the input 
side of one of the resistors to the output of the amplifier. 


6-9 


Model 


Model 3324 — Dual Inverter, Quad Electronic Switch 





3324: Electronic switch, four high speed elec- 
tronic switches with two uncommitted 
operational amplifiers. 


+100 -100 
eo ©., 





Electronic Switch Block Diagram 


(Connections shown for comparator amplifier) 


The center section of the Model 3324 contains four 
electronic switches. There are two normally open (NO) 
and two normally closed(NC) switches. The normal volt- 
age required to excite the switch is a nominal +28 volts; 
however, voltages of up to +100 volts will not damage 
the: switch. The current required to excite the switch is 
less than three milliamperes at 28 volts. 


The Electronic switches were designed and are nor- 
mally used to replace the reed type relays in the inte- 
grator modules. The Electronic switches may be used as 
normal switches as long as the input voltage and current 
through the switch limitations are observed 


Patching to the Model 3320, Dual Integrator is shown 
in the accompanying illustration. Patching to the Model 
3329, Dual Integrator is similar. Also included in the 
Model 3324 are two .001 mfd. capacitors and four .01%, 


10 K ohm resistors for use in high speed Sample-Hold 
circuits. 


The inverter sections of the Model 3324 are identical in 
operation and performance to those of the Model 3322A. 
Patching is identical to that of the Model 3322A and the 


instructions for that Model may be used. 


ELECTRONIC INTEGRATOR 
SWITCHES 





Integrator with Electronic Switches 


Digital Logic Modules 
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Model 3327 


Twelve three-input AND 
Gates with four relay drivers. 


Model 3326 
Twelve RST Flip-flops with 


four relay drivers. 


Digital control, the new way of multiplying the efficiency 
of an analog computer, can be included in the SD 40/80 as 
well as the small 10/20 computers. The SD hybrid computer 
expansion system is comprised of three types of compact plug- 
in modules. 


The advantages gained by digital control in an analog com- 
puter are of far-reaching significance. Here are some important 
new advantages made possible by SD’s Digital Logic Control: 

1. Track and hold operation by individual integrators. 


2. Sub-routines can be flexibly programmed at different speeds de- 
pending on decisions made by logical equations. 


3. Program statements can be arranged into a flow chart quite sim- 
ilar to those used in digital computation. 


The flexibility gained through this interplay of analog/ 
digital equipment results in: 
1. Better and greater problem-solving capacity. 


2. Ability to solve a wide range of problems that before could not 
easily be handled by an analog computer. 


3. Speed. Problem solving time is greatly reduced. 


Through digital logic control, sub-routines start and termi- 
nate when the corresponding binary control variables change 
state as logical functions of: 

1. External control (switches, relays controlled by external devices). 

2. The states of timers or sub-routine counters. 

3. Analog-comparator decisions. 


The interplay of binary control variables and analog com- 
putation results in a special hybrid analog-digital structure. 
Relays or electronic switches implement analog sub-routine 
changes under control of digital (binary) control variables 
and constitute the digital-to-analog interface of the computer. 
Analog solutions, in turn, can modify digital control. 


* When programming involves the use of digital control, logic levels for the digital logic elements are the 
opposite of those for the mode control and function relays. 


* Refer to Chapter 8, ‘Basic Operation of Digital Logic Elements,’® for complete operating descriptions 
of the Systron-Donner digital logic control modules. 


MASTER RESET 
a 
| ee at 


| c 


Time-event Control. 1 kHz 
clock with three decade 
dividers, three decade counters 
and one gated relay driver. 


Model 3326A 
Six RST Flip-flops, six AND 


Gates and four relay drivers. 





Additionally, in the SD system it is possible to insert a 
delay of the Hold-mode command to the Complementary 
Integrators. This delay is important because it enables the 
C Integrators to store the final values of their inputs in a 
problem solution where the C Integrators are tracking rapidly 
changing problem variables. 

The Reset-Compute-Hold modes are, of course, settable to 
any desired intervals, and it is this new degree of freedom 
imparted to integrators that illustrates the importance of 
digital logic control. 

The combination of the three SD digital logic modules 
(Flip-Flops, Gates, Time/Event Control) results in a most 
flexible hybrid analog-digital structure. The operator can 
easily set all Reset-Compute-Hold intervals of integrators with 
the Time/ Event Control Module. 

A typical example that illustrates the use of SD’s digital 
logic is mode control of iterative integrators. In iterative 
operation (IO), results obtained during or at the end of one 
solution of the problem are used to change parameters or the 
circuit configuration (switching) for the next solution. The 
following figure shows how in IO, integrators are paired into 
normal and complementary (opposite) logic to implement 
iterative solutions: 


R 
Normal \ H 


c 


R 
Complementary 1 H Pu 
c 


Mode Duty Cycle for Integrators in Iterative Operation 
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Function Generator Operation 


Principle of Operation: The Function 
Generator approximates the curve of an arbitrary function 
with a series of connected straight-line segments as 
shown in Figure IA. This graph represents the sum of 
the individual line segments shown in Figure 1B: Each 
line segment corresponds tothe output voltage of adiode 
netwark in the function generator. The point along the X 
axis where each line begins corresponds to the selected 
bias or ‘‘breakpoint’’ voltage. The diode networks are 
arranged so that six conduct when the X inputis positive 
(0 to +100 volts) and six when the X input is negative 
(0 to -100 volts). For better accuracy, the input voltage 
may be biased to use all 12 networks above or below 
zero. Two function generator units can be operated as 
a single-channel of 24 segments for more accurate 
simulation. 


The slope of each line segment corresponds to the 
ratio of output voltage per volt-change of input and can 
‘be varied between the values of +2 and -2 by adjusting 
the current flow through the diode network. For steeper 
changes of line slope at any point, two or more networks 
can be operated in parallel by shifting to the same 
breakpoint. 


When the individual outputs of all diode networks are 
added together by an operational summing amplifier, the 
total output over the entire range of X-inpuf will repro- 
duce the original: straight-line approximation. No output 
is produced when X = 0; therefore, a constant voltage, 
Yo, must be added at the output summing amplifier, when 
Y is not equal to 0 at X = 0. 


Set Up Procedure: A function may be set up on the 
function generator either by the preparation of a table of 
values to which the individual slope and breakpoint 
potentiometers are adjusted, or by adjusting the potentio- 


meters so that a trace of the output voltage on an XY 
plotter will match a prepared graph as closely as pos- 
sible (visual method), The second method may be faster 
but is limited in accuracy by the read-out instrument. 





A. Straight Line Approximation of a Curve. 





B. Outputs of individual Diode Bridge Network th 
of which produces the Straight Line Asproxinericn: 
shown in Figure 1A, 


Figure 1. Simulation of a Function with Diode Bridge Net- 
work Outputs 


TABLE 2. DATA FOR SETTING UP FUNCTION Y = COSX 


X (DEGREES) Y (COSINE) X (VOLTS) 


* Arrows indicate sequence of adjustments. 
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CONTROL* 


Y (VOLTS 
BREAK SLOPE 








Y (VOLTS) 





SREAKPOINT: -6 -S§ -@-3 -2 #14243 +4 48 46 


SCALING: X: 1 VOLT = 1° (INPUT) 
Y: 100 VOLTS = 1.00 (OUTPUT COSINE VALUE) 


Figure 2. Typical Graph prepared for setting up the Function 
COSsxX 


Tabular Method (for Model 3352) 


General: The maximum slope change for all segments, 
except the first, Kx, segment, is +2.5 : 1. The first, 
Kx, segment has a maximum slope change of +5: 1. 
Breakpoints may be stacked to achieve greater slope 
changes. 


Step 1: Set all Break-point/slope change (polarity) 
switches to the correct combination of polarities as de- 
termined by the function to be simulated. (See previous 
discussion on tabulating a function.) 


Step 2: Set all Slope potentiometers to approximately 
center position. Set all Break-point potentiometers to 
one end as determined by the arrows associated with 
the polarity switches. (E.g., if the polarity switch is in 
the +BP, +AS or in the -BP, -AS position the break- 
point potentiometers should be set to the counter-clock- 
wise end.) During this procedure the output amplifier 
may go into overload. Although no harm is being done, 
the computer may be placed in the Pot Set/Bal mode to 
prevent overload. The computer must be returned to 
Reset before proceeding with the adjustments. 


Step 3: Set the Kx/Yo switch as determined by the 
function and by the following instructions: 


a. The Yo potentiometer sets the value of Y when X 
equals zero. This value may be set anywhere within 
the range of plus to minus 100 volts. 

b. The Kx potentiometer sets the slope of the first line 
segment. The Kx slope is adjustable within the range 
of plus to minus five (5). 

c. If the Kx/Yo switch is placed in the Off position 
both Yo and Kx will be zero. 


d. If the Kx/Yo switch is placed in the Kx position Yo 
will be zero, Kx will be adjustable over its range. 

e. If the Kx/Yo switch is placed in the Yo position, 
Kx will be zero, Yo willbe adjustable over its range. 

f. If the Kx/Yo switch is placed in the Kx & Yo 
position both Kx and Yo are adjustable over their 
ranges. 


When the switch has been set, the Yo adjustment 
should be made. With X = O, the potentiometer may be 
adjusted for the proper value of Y, at the output of the 
Model 3352 (O,). The Null voltmeter or a precision ex- 
ternal meter should be used for this and all succeeding 
voltage adjustments in this procedure. 


Step 4: The Kx adjustment should be made next. 
Apply the value of X at the first break-point (BP) to the 
input of the Model 3352. (If a potentiometer is used to 
obtain the X value its output should be applied through 
an inverter to the Model 3352. This will prevent errors 
due to the diode function generator loading the potentio- 
meter.) Adjust the Kx potentiometer until the proper value 


of Y is obtained at the output of the Model 3352. 


Step 5: The first break-point (BP) is now adjusted. 
Adjust the first break-point potentiometer until the Kx 
value previously set is offset by approximately 150 milli- 
volts. (0.15 volts) This offset compensates for the diode 
characteristics and will help produce a more accurate 
function. 


Step 6: Set the input to the value of X at the second 
break point. Adjust the slope potentiometer associated 
with the first break-point until the output is equal to the 
value of Y at the second break-point. 


Step 7: Adjust the B.P. potentiometer associated with 
the next segment until the value of Y set in the preced- 
ing specification is offset approximately 150 millivolts. 
(See discussion in Step 5). 


The remaining slope and break-point potentiometers 
are set alternately as indicated by the large arrows on 
the panel, by repeating Steps 6 and 7 until all segments 
have been adjusted. 





CHAPTER 7 
LOGICAL ALGEBRA 


It is necessary tohave a rudimentary understanding of 
logical algebra in order to use logic elements. The 
fundamentals of logical algebra are presented for this 
reason. 


LOGICAL VARIABLES 
A logical variable, L, has two values: 
1 (true) 


0 (false) 


ADDITION (OR) 


The logical sum, Lg = L, + L,, is defined by the fol- 


lowing table: 


2’ 





This is commonly called the OR function and has the 
alternative symbology 


MULTIPLICATION (AND) 


The logical product, L, = L, « L, is defined by the 
; 3 1 2 
following table: 





This is commonly called the AND function and has the 
alternative symbology 


COMPLEMENT 


The complement of Lis denoted by L. Loosely speak- 
ing, L is the inverse of L. That is, whenever L = 1 
or 0, then L = Oor 1, respectively. In manipulating 
logical formulas, the two obvious relations below are 
often useful. 


L+L=1 


b| 
i 
° 


L . 


COMMUTATIVE LAWS 


L, + L, = L, + L 


ASSOCIATIVE LAWS 


L, + (Ly + L 


3) = 


1 1 3 
Ly * (lg * Lg) = (Ly Ly) * bg 
DISTRIBUTIVE LAW 
L,° (Ly + Lg) = Lj Ly + Ly Ly 
COMPLEMENTED SUM 
Lt Ly ~ Ly Ly 





L, °° bt, = Lb, + by 


7-1 








DE MORGAN FORMULAS 
SE. SL 


me 
iw) 
— 
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These can be obtained by complementing the comple- 
mented sum and product formulas. 


EXCLUSIVE SUM 


The exclusive sum, @, of logical variables is defined 
by 


@ Lig = ALy + Ly) (Ly > Ly) 


This is also called the exclusive OR. Ly ® Ly is the 
shaded area of the ''set'' diagram below. 
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In other words, the exclusive sum is the sum L. + L 
minus the product Ly - Lg. A more convenient ex 
pression for @ can be derived. 


—_—_— __. 


Ly @ Ly = (Ly+hy) » (Ly: Ly) 


(Ly+Ly) + ©, +L.) 


+ 


Ly: L 


CHAPTER 8 
BASIC DIGITAL LOGIC ELEMENTS 


Digital logic can be included in an analog computer. 
It is used for making program decisions based on events 
occurring during the analog simulation. It is also used 
for program mode control. The application of these 
elements is discussed in other chapters of this hand- 
book. This chapter describes the basic operation of 
the logic elements. 


OPERATION OF DIGITAL LOGIC IN AN ANALOG 
COMPUTER. 


Logic levels for the elements are 


+28 volts false (logical zero) 


O volts = true (logical one) 


Most of the elements have a logically complemented 
output in addition to the normal output. All of the ele- 
ments except the drivers and the dividers have status 
lights. When the light is on, the element is in the 
"true" state. 


All of the elements, with the exception of the drivers, 
are designed so that their outputs can be connected 
together without damage. This feature is very impor- 
tant because it permits an OR operation by connecting 
two or more outputs to the same point on the patch- 
board. (See driver exception below. ) 


However, when two outputs are connected together in 
this way, they cannot be used as separate logical var- 
iables. To avoid this difficulty diodes can be used. 
Thus, if Ly, Lg are two logical outputs and Lj + Lg, 
Lj, Lg are all required for logical inputs elsewhere 
in the program, the circuit below can be used 


Ly Ly 
Ly + Lo 
Lo Le 


Each of the logic elements will drive any reasonable 
number of others. However, the logic elements are 
not capable of activating analog control relays or func- 
tion relays. For this purpose drivers are supplied. 
The driver acts as both a power buffer and logical in- 
verter. It should be noted that the logic levels for 


mode control and function relays are the opposite or 
complement of those for digital logic. That is 


0 volts false (relay de-energized) 


+28 volts = true (relay energized) 


Thus, the driver not only delivers the current neces- 
sary to operate several relays but also provides the 
logic level inversion required. To use a driver, the 
output of the logic element is connected to the driver 
and the output of the driver to the relay input. These 
connections are made readily at the patchboard. It 
should be noted that the driver is a special element 
and its output cannot be used as an input to the flip- 
flops or DCU's. It can be used as a gate input. 


Important: The driver outputs must never be connected 
to any other digital logic output. Con- 
necting a driver output to ground (logical 1) 
will destroy the output transistor, that is, 
if connected to a gate, it should not be 
OR-ed with any other element. 


The remainder of this section describes the basic oper- 
ation of the digital logic elements. 


DRIVER 


The driver is used as a power buffer and logical in- 
verter from digital logic elements to analog mode con- 
trol, function relays, and electronic switches. Its 
logic levels are 


INPUT OUTPUT 
0 (false) 0 volts 
1 (true) +28 volts 


The program symbol is 


—$>—- 


GATE 


The logical gates are ofthe coincidence type. With the 
exception of input number 1, they operate as AND gates 
only on the inputs connected. That is, if no logical 
variable is connected to one of the inputs at the patch- 
board, then this input is essentially de-activated. Note 
that this is the same as connecting a logical 1 to the 
input. 


8-1 





Input number 1 is different because it must be connec- 
ted; otherwise the gate will not operate. However, if 
only input 1 is connected, the gate will act as an in- 
verter if its complemented output is used. 


The program symbol for the gate is 


2 (NORMAL) 
(INVERTED) 


where the input numbers are shown. 


The truth tables below will clarify the operation of the 


gate. 
INPUTS NORMAL INVERTED 
1 2 3 OUTPUT OUTPUT 
none any any 0 1 
INPUTS NORMAL INVERTED 
1 2 3 OUTPUT OUTPUT 
0 none none 0 1 
1 none none 1 0 
| INPUTS NORMAL INVERTED 
| 1 2 3 OUTPUT OUTPUT 
| 0 0 none 0 1 
| 0 1 0 1 
1 0 0 1 
1 1 1 0 
INPUTS NORMAL INVERTED 
1 2 3 OUTPUT OUTPUT 
0 none 0 0 1 
0 1 0 1 
1 0 0 1 
1 1 1 0 
INPUTS NORMAL INVERTED 
1 2 3 OUTPUT OUTPUT 
0 0 0 0 1 
0 0 1 0 1 
0 1 0 0 1 
0 1 1 0 1 
1 0 0 0 1 
1 0 1 0 1 
1 1 0 0 1 
1 1 1 1 0 
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The gate is shown on the patchboard as 





The true state response to inputs 1 and 2 is delayed 
by 50 usec, while input 3 is not. As will be seen in 
Chapter 9, covering logic applications, this permits 
the gates to be used in connection with flip-flops to 
form a shift register. Also, this feature frequently 
obviates the necessity of using one-shots for delays 
in a logic program. 


Each gate has three inputs on the patchboard. If the 
logical product of more than three variables is desired, 
the additional elements can be connected in through 
externally patched diodes, as shown below. 





The additional elements can be connected to any of the 
three patchboard inputs, except number 1. 


The gate can also be usedfor pulse-shaping. It is ON 
(logical 1 at the output) when no input voltage is greater 
than 4 volts and OFF when any input is greater than 
6.5 volts. When the net input is between 4 and 6.5 
volts, the output corresponds to the last state of the gate. 


GATE-DRIVER 


This is a special element found in the TIME/EVENT 
CONTROL MODULE. It isapower buffered, logically 
inverted AND gate which performs the combined func- 
tion of an AND gate and driver. The program symbol 
is 


Ly 
Le 
L3 


Co 


where Lj, Lg and Lg and eg are respectively logical 
input variables and the output voltage. This gate is 
the same in logical operation as the AND gate except 
there is no delay introduced at any input. When Ly: 
Lg:Lg = 1, the output of the gate is 28 volts; other- 
wise it is zero. 


The primary use ofthis gate is for implementing ana- 
log mode control changes based on logic outputs. For 
example, in this way a multiple decade counter can be 
used either to stop the analog computer after a pre- 
scribed number of runs or in conjunction with a clock 
to specify precisely the length of time the analog com- 
puter is inone of several modes. This is accomplished 
by patching the output of a gate-driver to either the C, 
R or H inputs on the patchboard. 


DECADE COUNTING UNIT (DCU) 


The DCU is either a decade divider or adecade counter. 
It contains ten internalbinary states which correspond 
to integersfrom0to9. Note that this is different than 
a binary counter which counts in the binary number 
system. There are two types of DCU's in the comput- 
er. Both of these are in the TIME/EVENT CONTROL 
MODULE (also referred to as Counter Module). One 
has all the integer binary states available as a "'ten 
line” output. This element is called a decade counter. 
A status light is associated with each of the integer 
binary states. The other type has only the decade out- 
put available and is called a decade divider. 


There are two master resets for eachgroup of 3 decade 
counters together with 3 decade dividers. These are 
the 0V and +V terminals. The OV input is directly 
coupledtothe clock. As long as alogical one is applied 
to this terminal the clock will be stopped. A state 
change from zero to one will cause all the counters 
and dividers to be reset. As long as a logical one is 
maintained at this terminal, all other reset signals are 


inhibited. 


A state change from logical 1 to 0 at the +V terminal 
will cause all the counters and dividers to be reset. 
(This is called a trailing edge trigger.) The +V input 
does not inhibit other reset signals, and it does not 
turn off the clock. 


All counter or divider inputs (including reset inputs) 
require a change from a logical 0 to a logical 1 for 
operation. Thus, the application of a logical 1 to an 
input will cause it to operate once. It will not operate 
again until the logical one is removed and reapplied. 
This is true alsofor the master resets for each count- 
er/divider group. 


The program symbol for the decade counter is 


N 


The counter outputs are the ten integer counts from 0 
to 9. Whenever a logical 1 is applied to the reset ter- 
minal, the counter will return to count 0. Whenever 
a logical 1 is applied to the input terminal, the counter 
will advance to the next higher count. The counter is 
cyclic: when it is in count 9 and receives an input 
pulse, it returns to count 0. 


Note that the counters can be used in series. This is 
accomplished by connecting any output from the first 
counter to the input of the secondcounter. Then, when- 
ever the transition from alogical 0 to alogical loccurs 
at the output of the first counter, this transition from 
0 to 1 is applied to the input of the second counter. 
This causes the second counter to advance one count. 
Normally, the second counter receives its input from 
the "0" digitof thefirst. Then, whenatransition from 
count 0 to count 1 occurs in thefirst counter, the out- 
put from count 0 goes from a logical 1 to a logical 0 
and the second counter will not be advanced in count: 
its inputis activated only withalogical changeof state 
from zero to one. Note that when the first counter 
goes from count 9 to count 0, then the count 0 output 
goes from 0 to 1, and the second counter is advanced. 


The program symbol for the decade divider is 


- DIVIDER OUT 


The divider has ten internal integer counts inthe same 
way as the counter. However, these are not available 
as logical outputs. An input to the reset terminal will 
cause the divider to return to count 0 internally. The 
output willonly return to state Oifitis in state 1. Ifthe 
logical variable applied to the input terminal changes 
state from 0 to 1, then the divider will count inter- 
nally. The divider is cyclic: count 0 follows count 9. 


With a sequence of input pulses, the output is alogical 
1 except when the internal counter is incount 8 or count 
9. Thus, when the internal transition from count 9 to 
count 0 occurs, there is a transition from a logical 0 
to 1 at the output which can be used as input to other 
logic elements. However, it should be remembered 
that if the divider is reset during the time when inter- 
nally it is in count 8 or count 9, the transition from 0 
to 1 will occur at the output then. 


CLOCK 

Each TIME/EVENT CONTROL MODULE provides a 
1000 cycle clock output at the patchboard. 
Maintaining a logical 1 at the OV master reset termi- 
nal will stop the clock. 


The program symbol for the clock is 


o= 


The logical cycle of the clock is repetitively a logical 
1 for 800 » sec and logical 0 for 200 » sec. 
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FLIP-FLOP 


A flip-flop (FF) is a one-bit memory. It has a set (S) 
and a reset (R) state. The set state corresponds 
to true (logical 1) and the reset state to false. The 
program symbol is 


—| 


There is a status light connected to the S output so 
that when the FF is in the S state the light will be on. 
The FF can be placed in the Sstate by applying a logi- 
cal 1 to the S input. It can be placed in the R state by 
applying a logical 1 to the R input. Applying a logical 
state change from 0 to 1 to the T input will cause the 
FF tochange to the opposite state. The flip-flop Rand 
S inputs are direct-coupled. Each can be a-c coupled 
by inserting a 200 pF capacitor in series with the input 
so that these inputs willbe activated by only logic state 
changes from 0 to 1 and not by logic levels. 


Internally, the application of a logical 1 to the S input 
turns off the R output, and a logical 1 at the R input 
turns off the S output. The R and S inputs are level 
sensitive in contrast to the T input which is activated 
by state changes. Consequently, if a logical 1 is ap- 
plied to both the R and S inputs simultaneously, both 
the R and S outputs will be a logical 0, and a state 
change at the T input will have no effect on the state of 
the FF. If both R and S inputs are a logical 1 and then 
both are changed to a logical 0 simultaneously, the 
final state of the FF will be either R or S with nearly 
equal probability. If the change to 0 does not occur 
simultaneously, the input which sees the logical 1 last 
will control the state of the FF. 
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If either the Ror Sinput is energized by a logical 1 and 
if a pulse or logical variable state change is applied to 
the T input during this time, the FF will act like a one- 
shot mono-stable multivibrator. Thatis, it will make 
a transition to the opposite state for about 5 usec. 


Two master reset terminals are located at the patch- 
board for each FF module. One will put all the FFs 
in the R state for a transition from 0to 1. The other 
reset is level sensitive. As long as a logical one is 
applied all the FFs are constrained to be in the Rstate. 


COMPUTER CONTROL LOGIC OUTPUTS 


The analog control logic outputs in the integrator mod- 
ule must not be used as an input to any logic element 
except a gate. Therefore, special F-RT and R-H out- 
puts are located in the FF and gate modules. These 
can be used as an input to all of the logical elements. 
They are identified as F-RT' and R-H'. The logic 
levels for both F-RT' and R-H' areconsistent with the 
digital logic elements. There is a logical inversion 
between F-RT and F-RT' as well as R-H and R-H'. 
For example, when F-RT' is connected to a flip-flop 
R input, the flip-flop will be in the reset state when 
The F-RT bus is energized. When the F-RT bus is 
de-energized the FF is free to respond to signals at 
its T and S inputs. 


The F-RT and R-H busses are logical complements, 
which is the reason for their selection in logic control. 


LOGIC SWITCH 


There is a manual logic switch in the counter module 
which is convenient for single-stepping programs. The 
output of this logic switch (ungrounded terminal) can- 
not be connected to any other logic output. The reason 
for this is because the output of the logic switch is 
shunted internally bya large capacitor. Consequently, 
if this element is connected to another logical output 
it will cause degradation of the rise-time of that ele- 
ment. 


CHAPTER 9 
CIRCUITS FOR SIMPLE LOGICAL FUNCTIONS 


Several fundamental circuits which are used frequently 
for logical programming are presented below. In all 
cases the trigger is a negative-going voltage. 


ONE-SHOT 


A one-shot circuit will produce a logical pulse of pre- 
determined length whenever its input is triggered. 
That is, before the trigger, its output willbe a logical 
zero. When the trigger occurs its output changes to 
a logical one for a fixedtime interval. A one-shot can 
be constructed with logic elements in several ways. 


5 psec: 


The circuit below will provide a 5 ysec logical pulse 
when triggered. 


Lo 


If a second trigger occurs during the pulse, then the 
pulse will last for 5 usec after the second trigger. In 
other words, after the one-shot is turned on by the 
trigger, it willstayon until 5 usec after the last pulse. 
The sensitivity to multiple triggering canbe eliminated 
with a gate: 


Ly 2) 


The logical pulse from this circuit will be delayed 50 
psec due to the delayed input of the gate used for L,- 


R 
sb 


S Lo 


50 usec: 


The circuit below will provide a 50 usec pulse when 


triggered. 
C : 


This is due to the 50 psec delay at input 1 of the gate. 
Its behavior for multiple triggering is the same as for 
the 5 psec one-shot. 


Variable: 


S Lo 


This one-shot will produce a logical pulse of fixed 
duration. The pulse duration depends on the choice 
of R, C. Its behaviorfor multiple triggering is some- 
times erratic. The sensitivity to multiple triggering 
can be eliminated by gating. 





IN 
COUNTER 
R 





This one-shot will produce a logical pulse whose dur- 
ation is determined by boththe clock frequency and the 
output line number fromthe counter. The pulse dura- 
tion for the circuit shown is n times the clock period. 
Obviously several dividers may be inserted in series 
with the clock if desired and several counters may be 
used. This circuit is insensitive to multiple triggering. 


HALF-ADDER 


A half-adder is acircuit for implementing an exclusive 
sum, ® : 


Ly 








Although digital inverters are shown in this circuit, 
they are not usually necessary because all logical ele- 
ments except for the counters have both complemented 
and non~complemented outputs. 


SHIFT REGISTER 


A shift register is a binary register which is shifted 
one bit each time a trigger is applied. The triggered 
half-adder is needed for this circuit shown at right. 


If Ly ® L2=1, then whenever a trigger, T, is ap- 
plied, Lo will be a logical pulse whose length is equal 
to the duration of T. Note that T is applied to input 3. 


If the duration is less than 50 usec then changes of Ly, 


Lg after the application of T will have no influence on 
Lo due to the built-in delay for gate inputs 1 and 2. 


The circuit for a shift register is 








Notice that the inverters for the triggered half-adder 
are not required since complemented inputs are avail- 
able. An initial setting for this register can be entered 
by resetting the flip-flops individually, or with the com- 
mon reset, and subsequently gating a logical one into 
each S input as required. This circuitry is not shown 
because it requires little explanation. Whenever a 
trigger, T, is applied which is less than 50 usec (a5 
usec one-shot is a convenient source) the following 
state transfers will be made 


L, —-FFI1 


FF, —FF2 


FF —>FFn 
n- 
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Thus, when each trigger occurs, the current value of 
Lj is entered intothe first flip-flop. The last flip-flop 
can be connected to the first flip-flop, instead of Ly, 
to obtain a register with an end-around shift. 


UP-COUNTER (BINARY) 


The circuit for the binary up-counter is (shown for 
three bits) 


FFI FF2 FF3 
R R 
Li T T T 
S S S 


The counter is initialized with the common reset. Al- Number of 
ternatively, an initial value canbe entered in the same FF NUMBER State Changes 
manner as for the shift register. This circuit is acti- 3 2 1 for Ly 

vated by a logically positive change of state. It will —— 
count the number of positive changes of state in the 


binary system. Each time FF n-1 returns to the R 0 1 1 3 
state, the state of FF, is changed. The table below 1 0 0 4 

shows the states of the FF's for successive logically 
positive changes of state of L. 1 0 1 5 
1 1 0 6 
Number of 1 1 1 7 
FF NUMBER State Changes 0 0 0 8 

3 2 1 for L, 
0 0 0 DOWN-COUNTER (BINARY) 
: The circuit for the binary down-counter is (shown for 
0 2 three bits) 


INITIALIZE 





INITIAL VALUE 


This circuit is activated in the same way as the up- 
counter. The table below shows how it counts. (It is 
assumed that the initial value is the binary number 111. ) 


FF NUMBER INPUT 
3 2 1 STATE CHANGES 


Ke OF OC Oo OF FY Ff & F&F 
- oOo Oo Fk FY CO CO KH 
ep OF OF OF Oo - 
aN moar WO NY KF OC 








TIME INTERVAL MEASUREMENTS USING DIGITAL 
. LOGIC 


The digital logic modules, which are part of the Systron- 
Donner Analog Computer complement, can be used to 
obtain measurements of time intervals with a resolution 
ofupto one millisecond. The overall timing accuracy 
depends on the number of digits displayed. It can be 
as good as .01%, which is a significant improvement 
over the typical 3% accuracy obtained through the use 
of an oscilloscope. Also, this particular capability 
of using digital logic modules to make accurate time 
interval measurements can eliminate the need for a 
frequency counter, which is generally an extra cost 
accessory. 


The basic time interval measuring program is given 
in Figure 1A, and the timing relationships are shown 
in Figure 2. This program measures the duration of 
a logic signal, L. Initially, all flip-flops and DCUs 
(Decimal Counting Units) are inthe Reset state. When, 
and only when L changes from 0to1, Flip-Flop 1 is 
toggled into the Set condition. Gate 3 is opened and 
the clock timing signals are sent to the Time Interval 
Indicator. This is a cascade of DCUs of sufficient 
quantity to measure the total time interval. At the end 
of the time interval being measured, L returns to the 
0 state. L from the inverted output of Gate 1 resets 
F-F1. F-F1, in turn, toggles F-F2 into the Set state 
and turns off Gate 1. The Time Interval DCUs stop 
counting and hold the elapsed time display. The Set 
output of F-F2 inhibits the output of Gate 2 through 
Diode 1, preventing any further changes in L from 
toggling F-F1. F-F2's Set output also opens Gate 4, 
and clock signals are sent into the Display Time DCUs. 


The Display Time Generator is also acascade of DCUs. 
Outputs from the Display Time Generator are logically 
summed in Gate 5. When the count reaches the value 
corresponding to the desired display time, Gate 5 is 
changed to the 1 state and toggles F-F3. F-FS3 is 
connected as a5 microsecond one-shot. That is, F-F3 
goes into the Set state for 5 microseconds before 
returning to the Reset state. The Set output of F-F3 
resets all flip-flops and DCUs, and the circuit is ready 
to measure another time interval. The automatic 
Display Time Generator can be dispensed with, and 
the system can be manually reset by using the logic 
switch on the clock module (Model 3328). 


If L changes state before the end of the display time, 
nothing happens due to the Inhibiting action of F-F2 and 
D1. If the reset action at the end of the display time 
occurs while L is in the 1 state, the system remains 
in Reset. Interval measurement will not commence 
again until L changes from 0 to 1. This is illustrated 
in the Timing Sequence diagram, 3rd interval time. 


To summarize, this program takes no action until L 
changes from 0 to 1. It then measures the current 
time interval, holds the display until automatically or 
manually reset, and then waits for the next state change 
from 0 to 1 in L. 


An alternate program is shownin Figure 1Bfor meas- 
uring the time interval between two independent logic 
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Signals. Lg initiates the timing action. Gate 2 acts 
as a buffer so that Lg is not affected by the inhibiting 
action through D1. Lp terminates the time interval 
when it goes into the1 state. The interval terminates 
at that time independently of the state of Lg. If Lp 
causes trouble by being inthe 1 state during part of the 
interval to be timed, it can be connectedto a one-shot 
Similar to F-F3 and then into the Reset input of F-F1. 


The program of Figure 3 allows measuring of the Nth 
time interval after reset. It is identical tothe program 
of Figure 1 except that a pre-set event counter is in- 
serted between L and the toggle input of F-F1. The 
event counter can be made of DCUs, or a cascade of 
flip-flops. 


The period of N intervals can also be measured, by 
inserting the pre-set circuit between L and the Reset 
input of F-F1. This is useful for measuring equal- 
period intervals with an effective resolution of 1/N 
milliseconds, or for averaging readings over several 
cycles. 


Figure 4 shows how to convert analog signals to logic 
levels. The logic gates can be operated from analog 
inputs; the output will conform to logic states. How- 
ever, the gates have hysteresis. That is, the gate is 
off for inputs greater than about +8 volts. Once opened, 
the gate will not turn on until the input becomes less 
than about +4 volts. Negative inputs to the gates cause 
no action. Also, the gates are not damaged by input 
voltages within the +100 volt computing range. 


Gate B is turned on at the start of the time interval 
period, thereby toggling F-F A into the Set state. 
Gate B turns off during the interval, but this has no 
effect on F-F A. When Gate B turns on again, F-F A 
resumes its Reset state and the timing is terminated. 


Figure 5 indicates the method used to convert inte- 
grator mode logic bus signals to logic levels. Gate 
input 1 must be grounded or connected to a logical 1. 
When the mode bus is energized, it turns the gate off. 
L then is taken from the inverted output of the gate 
since this is a logical 1 when the mode logic bus is 
energized. 


Display of the less significant digits of a time interval 
(when the more significant digits are known), can 
appear on the Model 3328 which allows only a 3-digit 
display. A3-digit display of the less significant digits 
will suffice when the more significant digits are known. 
For example, to read an11 second interval to a reso- 
lution of 1 millisecond requires a 5-digit display. If 
the first two digits are already known, the last 3 can 
be indicated on the clock. The counters are cyclic, 
so overflows of the more significant-digit counters do 
not affect the count in the less significant digits. 


As another example, the period of al radian per second 
sine wave can be measured to 1 millisecond resolution, 
and .015% accuracy, without displaying the most signif- 
icant digit since the period is knowntobe 6. + seconds. 
The display would indicate the .283 part of the 6. 283 
second period. 
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Figure 1A. Time Interval Measurement Program 
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Figure 1B. Alternate for Interval Between Two Logical Events 
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Figure 2. Timing Sequence 
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Figure 3. Alternate To Measure Nth Interval After Reset 
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Figure 4. Program to Convert Analog Signal to Logic Signal 
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Figure 5. Program to Convert Integrator Mode Logic (+28V) Bus to Logic Signal 


CHAPTER 10 
CIRCUITS FOR SIMPLE LINEAR AND NON-LINEAR FUNCTIONS 


Several simple operations are used frequently tor pro- 
gramming. They include division by a constant, dif- 
ferentiation, time delay, and a variety of sectionally 
linear transfer functions. These are described below. 
(It should be remembered that all electronic switches 
in this section must be used as‘J inputs in series with 
the input resistor. ) 


DIVISION BY A CONSTANT 


The circuit 


IM 
e, 
Co 
O.1M 
eo 
@ 
The symbol 
2 
e e 
I 10 0 
eo 


satisfies the condition 


ae, + e; + 10e, = Q. 


(The input gain labels are usually omitted if they are 
unity.) The transfer function is 


This circuit is useful for dividing a sum of variables 
by the same constant, a. Sometimes it is more con- 
venient to set a pot to a, if ais anatural parameter 
intheprogram. Otherwise the pot would be put between 
the amplifier output and the following input and would 
require the calculation of a7! for each new setting. 


This configuration is alsouseful for obtaining gains ot 
greater than 10 without requiring more input resistors 
or more amplifiers. 


DIFFERENTIATION 


Programs are usually written to avoid differentiation. 
This is done because differentiation inherently stresses 
noise and because the circuits used may lead to insta- 
bility. However, in some instances it is difficult to 
avoid forming derivatives. This is especially true for 
on-line data processing or editing applications. For 
this reason, it is helpful to be familiar with available 
differentiation techniques. 


A single amplifier program which can be used is 


R 


Cc 
e, — / Co 


which has the transfer functicn 


©o 1 

ec we peace. 
r is usually about 10K producing reduced gain at fre- 
quencies above (2nrC)"* cps. Since C is generally 
large, the circuit tends to load unduly the preceding 
computer element. Also, a large C may lead to sta- 
bility problems inthe amplifier driving the differenti- 
ator. The smallresistor, r, helps prevent drive cir- 
cuit instability. 


The more sophisticated program 





has the transfer function 





£0 _ _Ss 1 + Ts 
. "EE: .. 
1 12 l-a 
1 feed 
ee 
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where T = CRip. Rin, the input resistance at ej, is 
usually 1 Ma. Thus, expressing C in microfarads, 
T (sec.) = C(uf). 


Ideally, 


a=1- K,K,C 


but in practice@ is adjusted toa slightly different value 


due to amplifier and potentiometer phase shift. The 
error of the circuitis strongly dependent on a causing 
its setting to be critical. C must be selected so that 


a is large but not close to unity. This is due to the 
stability of the loop composed of amplifiers 1 and 2, 
and pot 3. A good choice for C is 


0.3 
C =) 
Ky Ky 
The circuit will give good performance for 10< K,K9 
<1006. The lower limit is necessary in order to limit 


the size of C and thus avoid stability problems. It fol- 
lows that the setting fora is roughly 0. 7. 


The test circuit for adjusting ais (the block denotes 
the differentiator circuit), 


DIFFERENTIATOR 


OSCILLATOR 


The oscillator is set to the highest frequency of interest 
and is adjusted until the error is minimum, using a 
scope. In practice a>0. 7. 


An excellent low-frequency differentiation circuit is 





It has the transfer function 


& .ils + 1 


Go rete % 
.00ls° + .Ols + Ky 


pot 


and at low frequency 
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This circuit works well for 


Ky 
1=—- <= 4100. [If desired, 
Ko 


pot 1 can be adjusted precisely with the previous test 
circuit, in which the integrator input is set to 


Ko 


Ki 


TIME DELAY (Continuous) 


Sometimes a program for time delay is needed for 
solvinga problem. The time-domain definition of time 
delay is 


e,(t) = e,(t - T). 


This has the transfer function 


oO 


-sT 
=e . 


| 
me |O 


There is no convenient way of generating a continuous 
time delay withonly analog computer elements. How- 


ever, various approximations can be made; most of which are due to Padé. The simplest (first order) is 


wn 


e 1 - 
ose 


CT” eg 


bora] nofr3 
yg 


which works reasonably well for a delay of ¢<0.6 where ¢ is expressed in radians. Since ¢= wT, where w is the 
frequency, 
0.6 


Ww 


T = 


That is, at very low frequencies, the approximation will provide relatively large time delays with reasonable 
accuracy. The program is 


e| e 


2 
; 


In practice, the highest significant frequency in the simulation is estimated in order to determine an upper 
bound for T. 


A second-order approximation which is reasonably accurate for ¢$ <1.8 is 








The program is 





A fourth-order approximation which is reasonably accurate for ><7.5 is 


1 3 2 1 3 1 4 
e 1 - 9(Ts) + 39(Ts) - g4(Ts) + T680(T) 


lo 





oO 
— 


1 1 ina\3 1 4 
1 + 5(Ts) + oa(Ts)* + gq(Ts)" + zegqg(Ts) 
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_ The program is 





DEAD- ZONE 1 


The dead-zone operation relates ep & by the sectionally linear function 


slope = -—m 


The analog program is 


+100 


-100 


d In the following circuits note that mode relays and electronic switches are operated by logical zeroes. The 


circuits have been shown with electronic switches for mode control. Function relays can be used instead of 
electronic switches, if high-speed operation is not required. 
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Pot 1 is set so that its diode will not conduct until e, 2 a. Pot 2 is set so that its diode will not conduct until 


e, = -b. Thus, for -b < e, <4, neither diode will conduct and e,= 0. Also if pot 3 is set so the total gain is m 


then 


oO 
I 


-m (e, - a), e,2a 


= 0 , “bse,~a 


-m (e, + b), e; =-b. 


Note that pot 3 cannot be set to m exactly, since the source impedance due to the pot-diode networks is not zero. 
If digital logic is used, the program is 


ELECTRONIC 
SWITCHES 


(a2 Os 
E> KE KE 





-100 
mb 
100 
+100 & 
ES3 DEAD ZONE 
U3 
L, U 
N.O. =Normally Open 
C2 OU = 
N.C.=Normally Closed 
me, L, L, Lo U, Uy ES, ES» ES, ey 
-mb < me, < ma 0 1 1 1 0 1 O O O 0 
+, >ma 1 0 1 0 1 1 C Cc oO -m(e, - a) 
-, <~mb 0 1 0 0 0 0 Cc O Cc -m(e, +b) 
Pots 2 and 3 are used to set the breakpoints for the ES, is open when Ly = Ly = 1, which is equivalent to 


functions. Their outputs are used both to set the com- 
parator references and to provide a biased input to the 
amplifier as required. The amplifier output will be 
zero when all three switches are open. Logical vari- 
ables L,, L, are defined by 


-b se, = 4. 


1’ “2 
Since pot 1 is set to m, eg will have the desired func- 
L, = sen (e, - a) (output of C1) tional relationship to e,. Normally closed and nor- 
mally open contacts of the electronic switches are used 
L, = sgn (e, +b) (output of C2) to provide the logic level inversion required. 











SATURATION 





The analog program is 


The saturation function is +100 


Co 


-100 


Pots 1 and 2 are set to limit at -a and b respectively. This circuit gives a "'soft" limit. That is, €9 does not 
remain constant at the saturation level, but continues to increase with increasing e; due to the non-zero feedback 
resistance when the diodes are conducting. 


The program with digital logic is 


m SATURATION 
e| C Kuo) FUNCTION 
aoe 
100 
+100—{_ ) (10> eo 
De 

















eo 
-100 dyes 
3 
U 
. b 
-m slope 
<7 Co 
~@a 
Li 
ca 
L2 
e, L, ly Ly U, UU, Uz ES, ES, ES, ey 
b a 
58 <= 0 1 1 1 1 0 Cc O O -me, 
+, ey a 0 1 0 0 0 0 O Cc O -a 
m 
-, @< > 1 0 1 0 1 1 O O Cc b 
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The operation of the electronic switches is governed by 
the logical equations 


3. =b < < a 

ES; ~ Ly Lo m <i. a 
=_ a 
: _ ib 
ES 7 Ly (ey ~ m 


The normally closed or open contacts of the electronic 
switches are usedto provide the logic level inversions 
where required. The logic program will provide hard- 
limiting;i.e., the slope after saturation will be zero. 
This is not true for the purely analog program. 
BINARY (bang-bang) 


The binary function is 


Co 


The analog program is 





Pots 1 and 2 are set to limit at -a and b respectively. 
As in the "soft"’ saturation function, the output in the 
limiting region has a non-zero slope. 


The program with logic is 





—100 


The NC (normally closed) contacts of the switch are used. 
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ABSOLUTE VALUE (variable slope) Co 


—m 
The absolute value function with variable slope is 7 m, 


The analog program is 








The program with logic is 





BACKLASH 


The backlash function is 
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The analog program is 





Pot 3 is set tot where m is the desired slope. Pots 2 - Sufficiently large C to ensure that the inte- 
1 and 2 are set to limit at -b and a respectively. The grator drift is permissible. 

limiting network and amplifier 3 comprise a dead-zone 

circuit whose output is an integrator rather than a sum- The logic program for this function is obtained by re- 
mer. Integrator drift will be about 50pv/sec/pfd, so placing the analog dead-zone circuit with the equivalent 
it is necessary to consider two factors when selecting logic program. 


the integrator capacitor: 
BISTABLE FLIP-FLOP 


1 - Sufficiently small IC time constant to permit An 
adequate tracking The bistable flip-flop function is 


The analog program is DITHER 


The dither function is 






BANG —-BANG 


Pot 1 adjusts the value of ey which results in a change 
of state, namely C. 








The analog program is 


nip phen re ek ee 


ob 


calle te 


emt Tl eget 





In both the analog and the logic programs, the gain of integrator 1 determines the dither frequency. 
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CHAPTER 11 
SIMULATION OF CONSTANT COEFFICIENT TRANSFER FUNCTIONS 


The simulation of control systems requires programs 
for constant coefficient transfer functions of the type 


T(s) « P(s) 


Q6) 


where P, Q, have polynomial form: 


n : m 
P(s) = >, 348 , As) > bjs ; 
i-0 j=0 


Ifn > mdifferentiation is implied. Differentiators are 
treated inChapter 10. Here, it is assumed that n < m. 
For physically realizable systems the aj, bj are all 
real numbers. It is well known that any polynomial 
with real coefficients can be represented as a product 
of quadratic factors. Further, in manypractical 
problems the physical model is composed of dominant 
second or first order transfer functions. For these 
reasons only cases for which m < 2 are considered. 
Programs for specific higher order transfer functions 
can be derived easily with the approach outlined in 
Chapter 3. The materialhere is intended to comprise 
a dictionary for secondor first order programs. With 
this information, generating a program for a complex 
system diagram with many blocks is easy. The pro- 
gram is written for each block and the system is rep- 
resented by making appropriate interconnections be- 
tween blocks. 


In developing transfer functions using only pots, am- 
plifiers, and integrators, the following fundamental 
characteristics of feedback or closed-loop systems 
are of interest. 


General Closed-Loop System 


—Q- SUMMING POINT SYMBOL 
+ 





The transfer function of this system is 


0__G 
E, 1-HG 





where G(s) is the forward path transfer function and 
H(s) is the feedback path transfer function. 


Ee 
es ‘dia dl Eo 
or 


The transfer functions of this system are 


go OE 
E, 1- HG, Gy 
St ie 


That is, the response to inputs other than E1 is equal 
to the normal response divided by the transfer function 
preceding the actual input. 


Specific transfer function programs are given below. 


T(s) =a: 


EO) 4 
Ey 


Ee) £ 
a 


T(s) a 











. Ts) 








as + 1: 
Ts) = <5 
T(s) = AS, 
T(s) = $39° 
Y 
to 
s + ags+8 
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E eS 
I a - Eo - 
OE, aS+l 


E Eo -g8 
Eo = 
E, Sta 








Eo_ Sta 
° E, St+B 





Tls) = 5 : 
s°+ @s+B 








For a, g > 9, a” > 48, another program is preferable. 


Let 6 =, ~2 = 62-8. Then 


T(s) = (s +6 + ye + 6 -A) 


and the program is 


S+taA 


The program for this function uses types (1), (2), & 


(3) of T(s) =——"_, shown above, with the 
s + as +B 
addition of two pots andone or two summers. In each 


of (1), (2), (3), pot 1 is set with y= 1 








2 
T(s) = “5 es 
S +as + B 


For convenience T(s) is simplified 


Ts) = 1 + “7s + 8) 


Ss +as + 6 


andthe programs are (where againextra circuitry has 
been addedto (1), (2), (3) above andin each pot 1 is set 
with Y = 1). 





c~ ~~ 
2 o 
b e 
u 
Oo} — o| — 
i | LW ij | WW 











CHAPTER 12 
CONTROL SYSTEM SIMULATION 


A modelof a control system is usually represented by 1 - derive a program for each block 
a blockdiagram The blockdiagram contains constant 


coefficient transfer functions, non-linearities, and 2 - interconnect the blocks. 


frequently, special differential equations. The easiest This method is called block programming and is simple 
way to generate the complete computer program for using information developed inChapter 4. Some illus- 
the model is to trative examples are shown below. 








MODEL: 
wese 1 
s+a “S+ B 





PROGRAM: 








_ MODEL: 








s+l SS ae 


; 
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ABSOLUTE VALUE 
CIRCUIT 


PROGRAM: 





Note that no pot is required for the parameter, B, be- 
cause this parameter is used to translate the dead- 
zone function. 
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CHAPTER 13 
FUNDAMENTALS OF VECTOR ANALYSIS 


The development of vector operators given below is 
within the framework of Euclidean (x, y, or x, y, 2) 
space. Due to the invariance property of vectors and 
vector operators under transformation to other coor- 
dinate systems of engineering interest, itis reasonable 
to proceed in this manner. Sometimes we will show 
the representation of the vector operator under discus- 
sion in other coordinate systems for the purpose of 
illustration. Indeed, the main motivation for using 
vector notation is its non-dependence on a particular 
coordinate system. 


DEFINITIONS 
(A - B) is the scalar (dot) product of two vectors. 
A-B = |A|- |BlCos (A, B) 


(A x B) is the vector (cross) product of two vectors. 


If 
C = AxB 
then 
|c| = |A|- |B] Sin (A, B). 
Note that 


(C x (Ax B)) = 0 
Vis called the vector operator del, 


y 42, 2, 2 
ax ay oz 


(Here ordered n-tuples are used instead of the cus- 


tomary ij k unit vector notation, i.e. i eT 


ax’ dy Ox 
dee j) 
ay) : 
If ¢ = ¢ (x, y, z) is a scalar-valued function, then 


_ 296, 96, ag 
Vo = ox’ dy’ oz 


is a vector-valued function also known as grad ¢. 
Vhas different forms in other coordinate systems. 


Example: In cylindrical coordinates 


For a vector V = Vx, Vy, Vz; 


pO Vig ns ON ys ON ge OM ae. 
WN Se go gas eee tee oe 
y Zz Zz x x y 


is a vector-valued function known as curl V. 


2 92 sd 9” 
V‘vVzyv = 9. ¢ ao to is known as the 
x y Z 


Laplacian. 


SUMMARY OF OPERATOR ALGEBRA (greek letters 
are scalar-valued functions, capital letters are vec- 
tor-valued functions) 


Vv: (aV) a(vV-V) + V- (Va) 


Vx @V) a(VxV) + (Va) xV 


Vv: (VxU) = U-: (VV) - V~- (vxU) 


V(V: U) = Vx (VxU) +(V°vU+Ux(vRV) + 


(U-vy)V 
Vx(VxU)=(U:-WV-(V:- YU-U(V: V) + 
viv U) 


Vx (Va) = 0 


V-(VxV) = 0 


vx (vx A) = V(v- A)- v2A (this defines V” 
which is not a true vector operation). 


AxB =- BxA 


A’ (BxC) = (AxB): C (triple scalar product) 


DIFFERENTIAL OPERATIONS 


The remainder of what follows will be derived for two 
dimensions but may be readily extended to three. 


Divergence: 


If f (x, y) is a scalar-valued function then the family 
f = c of curves in the plane are level lines for f. If 








f is thought of as a potential then these curves are 
equipotential lines. 


_ af of a _ 
df= 5, + 5, dy = 0 


Vf -: (dx, dy) = 0 


where (dx, dy) is a vector element of arc along the 
curve f = c. (This is so because we constrain the 
values of x, y so thatf = c). Thus we see that 


vf 1 (dx, dy), 
and the family of curves obtained from the differential 
tangent field, vf, are orthogonalto the family of level 
lines and represent a flux (or flow). Vf: N S a is 


calledthe directional derivative of f in the N direction 
(N is a unit vector) at the point (x, y). The meaning 


av av 
is obvious. V:V = a + an is a scalar-valued 
x 


function known as the divergence of the vector-valued 
function,V. It is commonly written as div V. 


Consider a small rectangle in the x, y-plane 


[ye 


(x,y) 


We think of V as a flux, and determine the net emer- 
gence of flux from the rectangle, 


Vy(x+&,y tk) 


(xth, yt+k) 
-Vx(x,y+l’) Vx(x+h,y+O) 
(x,y) 


-Vy(x+€,y) 


We have 


Outflow = [oe +, ytt)- V,05 +e) k 


+ [Myers yek)- V(x +e", ») h 


unit area — h 

Vi (xt+é, y+k) - Vi(xte', y) 

E> Se Se mre Henney 
k 


Outflux V,(x+h, y+t) 7 Vy. (x, yr t') 


+ 


In the limit as h, k=» 0 we have 


Outflux _ ove 


— + 


unit area ax = Vv: V= div V 


g| 2 


Thus if ¢ is a conservative potential field 


Vv: Vo = 0orv¢ = 0 


At.a point x, y; v2 ¢2 0 implies a local "07" © in the 
case of incompressible flow. Notethat with compres- 
sible steady flow V“ ¢~0, in general, without the 
presence of sources or sinks. In transient physical 
processes ¢ is a measure of the change of energy 
from one form toanother. For example, in heat trans- 
fer processes heat is stored locally by the distributed 
heat capacity of the medium in which the process occurs. 
Thus, if T is the temperature at (x, y), the local con- 
ductive heat transfer is k VT where k is the thermal 
conductivity, and the local heat stored by the heat ca- 
pacity of the medium is cAT where c is distributed 
heat capacity. The loss of thermal flux (heat flow) 
from the conductive process equals the heat flux re- 
quired for heat storage andthis time-varying relation- 
ship is given by 


vee Bek 
V:-kVT =c at 
or 
2p =f OT _1 2T pittusi 
VT -k at ~K oat (Diffusion 
equation). 


where K is the thermal diffusivity constant. 


Inthe case ofthe small oscillations of auniform mem- 
brane, if w(x, y) is defined to be deflection, then the 
potential energy stored by deformation of the membrane 
is a function of Vw while the local kinetic energy of 
motion due to the mass of the membrane is a function 

£ ow 
of , 
one state to the other is described by 


Thus the time-varying transfer of energy from 


(wave equation) 


where c is the wave propagation velocity. 
Curl: 


If we think of the vector-valued function, V, as a ve- 
locity field, we can interpret curl V in a convenient 


way. Consider the tangential velocity around the sides 
of a small rectangle. 


~Vx (x+&,y +k) 
(xt+h, yt+k) 


—Vy(x, y+’) Vy(x+h,y+6) 


(x,y) 
Vx(x+€,y) 


circulation : Si, ee 
The Sniba res (spin) around the boundary is given by 


Vy& +h, y+t) - V(x, y+¢é') 
wD ee 


h 
V,& +6", y) - Vi +e, y +k) 
ee ee 


In the limit as h, k—~0 we have 


. ; 8V oV 
circulation, . _ y x _ es 
anicaren (spin) = a Vx V = curl V. 


To visualize the result in three-dimensions, we think 
of an infinitesimal paddle wheel located at the point 
(x, y, z). This paddle wheel is angularly oriented so 
that its shaft has the greatest rate of spin. Using the 


righthand rule, the shaft then points in the direction 
of VxV and its angular velocity is proportional to the 
magnitude of VxV. 

If VvxV = 0, V is said to be an irrotational field. 

If V = vVxA, V is said to be a solenoidal field. 

If V = V: ¢then V is irrotational. (VxV- ¢ = 0). 
If V = vxAthen Vhas zero divergence (V: VxA = 0). 


VECTOR TRANSFORMATION THEOREMS 
Theor 1: (Divergence or Gauss theorem) 


Let R be a region over which a vector valued function 
v is defined, then (differentiability assumed) 


fvo-v-fvy-n 
R 8R 


where nis aunit exterior normalto aR (boundary of R). 
Theor 2: (Stokes' theorem) 
Let R be a 3-dimensional region over which a vector 


valued function V is defined, and let o be a 2-dimen- 
sional subregion of R, then (differentiability assumed) 


fn -cvevy= fv + as 
o 30 


where N is aunit positive (sense to be determined by 8c) 
normalto o andds is a direction vector which lies in ac. 





CHAPTER 14 
PARTIAL DIFFERENTIAL EQUATIONS — PART 1 


This section treats the class of partial differential Helmholtz: 
equations which can be represented 


> 
vs KV®@ = Bd, ¢= P(Ky,° °°; xy t) 
.OOl 
where the operation, =, is defined by one of 
=o = 0 (Laplace equation) + 
ae 


=o = k (Poisson equation) 
=o = ko (Helmholtz equation) 
Ed = C(X) ze (Diffusion equation) ae caer is added toprevent algebraic loop insta- 


2 Diffusion: (N is the time scale factor 
=o = M(x)? 2 (Wave equation) 
at 


d 
The analog computer is used to integrate =, in other 
words to generate gl This will be discussed before 
considering the generation of V- K V¢. The block pro- 

-H¢ 


grams for these cases are U 





C(x) N 
Laplace: 
> 
ool Wave: (N is the time scale factor) 
dp 
-H¢ 
| 
The high-gain amplifier forces Zp = 0 when connected MO)N2 
tothe remainder ofthe program. The capacitor is used (x)N 
to reduce high-frequency gain which may result in 
instability. 
Poisson: 
H¢ 
? ; nth 
With the blocks for =¢@ available, all that remains is 
to find how to represent V- KV@¢; the Laplacian of ¢. 
OO! The form of V will depend on the coordinate system 
(seeChapter 13). To illustrate the concept we will 
k derive the Laplacian for one dimension first and then 
Too develop it for both two-dimensional Euclidean and Rei- 
mannian manifolds. The extension to three dimensions 
—~#H¢ +100 will be obvious. For simplicity it will be assumed that 
o K is not a functionof time. For thecase where K is a 
function of time, the pot setto Kis replaced by a multi- 
The high-gain amplifier ensures =¢ - k = 0. plier. 
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ONE DIMENSION 


V->- KV = afew) & — ¢(x, 0]. 


The derivatives will be replaced by finite differences: 


Since 


f(x + Ax,t) - f(x,t) 


— x i, t) = lim x 


Ax~>0o 


then, for small Ax 


a = f(x + 4X ty = {(x+ Ax, )e - f(x, t) 


Adopting the notations 


f(x,» t) = fs for any n 


we have, for finite differences 


xXx,+xX t 
2 * X> ) 
ag (“Ba * 


7 f(x, t) - f(x,, t) 





Ax Xo - Xy 
For any n 
X41 7 xy Kil, t 
Af es Af nts’ : 
Ax "An " 
Af 1 
: nts : fx ap t) - f(x, t) - 
ox +1 . a 
= n+1_ en 
Xn ~ xy 
Thus 


Ax, 1, t) ~ Wx,» t) 


KV = K(x_1)- 
[ are +5 nts Xl - x 
M11 ~ 

= K 1 
ny +1 7 *n 
Similarly 

% ~ %y-1 
[Kye], 97> Kl ¢ TR 
n-5 2 n n-1 
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Using the same approximation 


[xve],, 41 - [Kv¢], 1 
“+5 ae 
Lv: Kvg], = = Ef, 
n x.1-x 1 
nts n-9 


Therefore, by substituting the different equations pre- 
viously obtained, we have 


Rus [tn+4 - $, | Kn-3 [tn - $. 1| : 


x 1 n-1 


The program which generates Zo, is 


-P g-) +o, —P nai 





If Xi *n Ax, for alln, then the program reduces to 


—Py_| ?,, —Paat 








K,1 Kt 
Pp, -4, p, «5, P, = 
1 ax ’?°2 Ax ’ "3 Ax’ 


Finally, we have the program for all five types of 
equations above. Two are given below as examples. 


Helmholtz: 





Wave: 





ONE-DIMENSIONAL BOUNDARY CONDITIONS The programs are shown below for the Wave equation, 
but are typical for all Eq. 
Boundary conditions of two types are commonly en- 
countered: 
Type 1: 





Ha, t) = f(t) [Type 1] 


K(a) V9 (a, t) = f(t) [Type 2] 


Sometimes both types are specified simultaneously as 
a mixed boundary condition. These are approximated 
by 


io f(t), x, = oO 


K 1 v1 = f(t), x1 = 
nig. tg aes 


TWO-DIMENSIONAL EUCLIDEAN MANIFOLD 


Type 2: 
A two-dimensional Euclidean manifold is commonly 
called an x, y-coordinate system. Here 


a Q 0 a 
ve Kv = (2 £) > (Kev) 3% Kt, $2) 


Pp 
and 
| V° KVod = 7 Sy By K By zo 
M ,NoAX 
As for the one-dimensional case, finite difference 
approximations are made for the derivatives: we will 
adopt the notation 
= x 
$n, mi “ n’ 7) 





f (t) 
AX 
Sate m P+, m Py m n-= ™ ry m 1-1, m 
x 1-x 1\ *p1%n x 1-x 1\*n“*p-1 
mt+y n-5 nty N-5 
Mixed boundary conditions are represented by 
K 
2 ny m9 % m+17 ?n, m 
s y “y 
Yn+3 a 1 m+l1 ’m 
—f Ch) 
K, m- - 
tek 2 n,m y m-1 ae 
y iby i Ym ~ Ym-1 myn 
=o. oe 
M,, aan 
if X17 *_, = Ax, Vag n= Ay, for alln, m, 
then the equation reduces to 
K 1 K_ 1 
nt5, Mm N-5> m 
fa(t) Ax (ns, m_ on, m| * Ax (0-1, m” *n, mn) 
Ax 
K 


K 1 1 
n, mt+z n, M9 
se Ay (6, m+1— on, mn) Ay (, m-1- on, ma} 


This discussion can be extended to higher-dimension 


cases by treating each of the boundary points in the 
same fashion. = Ax AY=o, tA 
3 
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The program is 





TWO-DIMENSIONAL RIEMANN MANIFOLD 


Let the metric for the manifold be given by (ds = arc 
length) 


ds* = gij dx'dx! . 


In what follows we will use summation notation (i.e., we 
sum on repeated indices). 


g = |gij| = g,,G° (no sum on i) 
sothat G” is the cofactor of g,; inthe determinant | gi |: 


We define 


gi _S" 
- 


Note that 


It can be shown, that if ¢is a scalar function, 


2 1 a ij a¢ 
eo AEA) (ve oll 28 
Ve ax? ox’ 


ox 


Thus, we have a general expression for V ¢ in non- 
orthogonal curvilinear coordinates. Wecouldtreat the 
non-orthogonal case. However, that treatment is much 
more complicatedthan for the orthogonal coordinates. 
Since the latter is of greater practical interest, we will 
restrict our attention to 





Bij 0, i ~ j, 


gl= 0, i= j 





Thus, in three dimensions, 


3\2 


2 1,2 2,2 
ds" = 1, (dx y+ Boo (dx yy + B53 (dx ) 


& = 811 899 833 


ii 1 ‘ 

= — (no sum on i) 
ii 

We denote 


ox’, x, x): x” = AY a=i, j, k, A, are constant 


by ¢. . 1. 
a Cas 
Consider the diffusion equation 


Vv: KV = 


We will need to replace the space derivatives with 
finite differences. 


1 3 ij a 
V°KV¢ =— -= (K Ve gu2e ) 
Ve ox’ ax! 


= £q K=K(x!, x2, x) 
~~ (Kvegl 22) = ve x 
ax ax! 


Now 


V811833 
V 899 

V8 11829 
V 833 


22 
Veg” = 


33 
Veg” = 


and 


& | 
B 
eo 

«1S 


re 
=) 
jaa 
-_ 
5) 
oo 
oo 
© 
ND 


aw 
n 
2 
ay 
a 
= 
nN 
ND 
a8 
it 
a 
hI 
@ 
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When approximated with finite differences 


KV822 V833 ag | 


ax? Vey ax ax! ax’ ax® 
axtax® [| KV8 a2 V833 
: 1 
a ve11 i+1/2, k, j 


2x2 [KVE_2 V833 
‘l%+1,jk %ik] ~ [od VER 
814 


- e gs _ . -¢. ° 
i-5, i, k i-1, j, k ‘i, j, k 


and similarly for the other terms. Let us consider 
the one-dimensional projection to determine the sig- 
nificance of these terms. From the definition of the 


metric, clearly 
1 
Ax V811 
2 Ve. 
AX VBo0 
3.7 
AX 833 


are lengths in the xi, e x directions respectively. 
We denote these by 1}, i” 1° respectively. ax! Ax 
ax? Ve is the volume V = 1/1°1°, 


In general, the value of the lengths and the volume for 


increments in x}, x’, x? depends on location. The one- 
dimensional projection of the problem is 


12 itgoik i+ 1,j,k ‘i,j,k 
[28 4 
i 1, \4-t, 5%, 5, i 


i- 2? ); 








NG pe Pig, 










o DEFINED HERE 


Thus, for example, we see that, in the Diffusion 
equation 


1713x 1213x 
I i and I 
i i- 5,1. 1 i+1/2, j,k 


are lumped conductivities defined at appropriate loca- 





tions and V; j k= is a lumpedcapacity, also defined at 
949 


an appropriate location. The extension to the non-pro- 
jected problem in three-dimensions is obvious. Itis 
interesting to note that the theoretical result agrees 
with our intuition. 


ASYMMETRIC LATTICE 


The set of space points over which the operator 
V- KVo 


is definedis calleda lattice. In arbitrary regions one 
cannot finda metric in any simple way. Thus, the ap- 
proach developed for the Riemann manifold is not ap- 
plicable. Ifa Euclidean system is used, then the bound- 
ary must be composed of acollection of orthogonal line 
segments. This has two disadvantages: the lattice 
density must be high on the boundary inorder to repre- 
sent the boundary. The lattice density on the boundary 
determines the lattice density in the interior of the 
region, leading to an exorbitant equipment requirement 
for the program. The asymmetric lattice representa- 
tion of V - KV will allow one to avoid both of these 
difficulties. 


Consider a lattice whose only restriction is that a set 
of perpendiculars to the branches (one for each) can 


@ DEFINED HERE 


. K wa 
; we x 
COORDINATE 


itl, j,k 


(2223) iad, isk 


(CROSS-SECTIONAL AREA) 


g@ DEFINED HERE 


be found such that the perpendiculars in each polygon 
enclosed by a branch loop intersect at only one point. 





Referring to the figure above, 


1. Let a point function ¢_ be defined over the . 
lattice . 


2. Consider one point of the lattice where we 
arbitrarily choose % and denote by e the ¢'s 


at branch connected neighbors 


3. Define £ to be the length of the branch between 


lattice points associated with o Pi 


4. Define r, to be the length of the _|_ to t 


r 


5. Define the symmetry ratio Y; = >= 
i 
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Now consider the representation of v6 at the lattice 
point 0 given by 


0(¢,) yA oh > (6,~ 4) 
a 


Where qa is the index set of the neighbors. We will 
choose the Y, and A_ suchas to minimize error. Now 
for the continuous medium (where we denote the loca- 
tion of lattice point i by S;) each (S,) is given in terms 


of o(S,) and the derivatives of ¢ at S, by 


8) = #8.) +E + @& 2 +y, 2" 48) 


i ox  Viay 


We have assigned aCartesian metric to the domain so 
that Le is the pair (x,, y;)- Thus the representation 


0(¢,) is given by 


ce 
1 8 8 \n 
Ly, u n! (x; ox i ay? 
a n=1 
or 
5 + yy. « 2 y 
n! i ax i dy 
n=1 a 
Thus 


0(¢,) ~ x YX %xo : Xo Yi yo 


a 2 a 2 
v2 2% x0 * 2 24 %yyo0 


+ 24 Ji 1) 


yO + higher order terms 


Now we require 
> Y;5; = 0, LY,y; = 0 (1) 
a a 


(1), above, is the x and y projection of }> ¥i4 ‘ 
a 


re ae 
a a 


TT ; projected on x, y is zerosince the polygon formed. 


a 
by the t is closed. Thus, (1) is automatically satis- 
fied. Evidently our choice of ¥; was Suitable. 
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ei ae 2 
By 3 LX » By 2 Yi¥; > Bxy DY 
a a 
Then 0 can be expressed 


0(¢,) = By bxo * By byy0* Pxy *xyo 


+ higher order terms. 


We want to reduce this to canonical form. The result 
above is independent of an orthogonal transformation 
on (x,y). Note that there exists an orthogonal trans- 
formation such that Bay = Q. Henceforth, without loss 


of generality, we will assume Pay = 0. 


Thus, 


0(¢,) = By Px t Pte higher order terms. 


B,+8 
0644) = (25D (exg + byyo) 


B,+B 
+ APe007 


dyy a} + higher order terms, 


We see that if By is significantly different in magnitude 
from By? then Ois apoor approximation of v. Thus 


we require the polygon formed by the perpendiculars 
to have relatively good symmetry properties. The 
higher order terms can be reduced by increasing the 
lattice density. Onecanget an intuitive feeling for the 
symmetry properties of the polygon by observing that 
the system has a physical interpretation in terms of 
moments of inertia as shown below. 





Y. = Mass at ¢. 
i i 


All that remains is to find the normalizing constant 
A. We choose 


+ 
A _ Px” Py 
Oo 2 
1 2.2 2 2 
Ay 4 LY, ty 0 = 22. ¥i% 


Geometrically, 





Interpretation: 


If £4 is a branch from % to dis then the lumped flux 
from polygon i to polygon j corresponding to V¢ is 
Ts 
— 
1. [+ . 

1) 1) 
tance. A, is the area associated with lattice point j. 


: 4] so that ¥j, = ; is a lumped admit- 


MacNeal, R.H., "An Asymmetrical Finite Differ- 


ence Network", Quarterly Appl Math,October 1953) 


-14-9/10 








CHAPTER 15 
BASIC ITERATIVE PROGRAMMING 


REPETITIVE OPERATION 


In Repetitive Operation (RO) the entire problem on the 
computer is solved repetitively. This is useful for 
determining the influence of parameters on the solu- 
tion. Potentiometers and switches are operated man- 
ually while the computer solves the problem repeatedly. 
The RO mode cycle is R, C, R, C, etc. The length 
of the R (initial condition or reset) and C (compute) 


R-H! 
(©) 


DIVIDER 
COUNTER 


mi 


COUNTER 


7 


COUNTER 


n3 


OV 
MASTER 
RESET 


intervals can be set approximately by the controls at 
the console. Usually theC interval is adjusted so that 
the problem solution time is long enough for one to 
observe the significant part of the solution. The R 
interval is set sothat all integrators return to correct 
initial values. The RO controls are calibrated only 
roughly. For a more precise setting of the IC inter- 
vals, the program below can be used with the RO 
selector switch at the console set to EXTERNAL: 


COMP TIME 


7X 


MASTER 


RESET 





F-RT' 
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The output fromthe F-RT' terminal resets the divider, 
counters, and clock, because the logical function 
(GATE 1) + F-RT' is connected to the OV reset ter- 
minal. Note that RH=F-RT'. The clock is off as long 
as the F-RT' terminalisa logical 1. The F-RT' ter- 
minal is a logical1 whenthe main machine mode is R. 
The FF is reset whenever 


(GATE 1) + F-RT' = 1 


When the FF is reset allthe integrators will be in the 
R (IC) mode since the R output of the FF is connected 
through a driver tothe R terminal. Each time the FF 
goes to the S state, the divider and counters will be 
reset momentarily due to the connection of the in- 
verted output of Gate 2 to the +V terminal. The FF 
will be set when Gate 2 = 1. The logical equations 
for the gates are, 


G1 


i} 
2 

58 
i) 

5B 
oo 


G2 


Thus, the reset and compute intervals are, 


R = 100m, + 10m, + m, 


100n + 10n +n 


Cc 3 2 


1 
times the period of the divider output. 


Note: The connection to patchboard C terminal turns 
on mode switch compute light during C portion of cycle. 


ITERATIVE OPERATION 


In iterative operation (IO), results obtained during or 
at the end of one solution of the problem are used to 
change parameters or the circuit configuration (switch- 
ing) for the next solution. The term JOhas an obvious 
derivation from mathematics. In IO, integrators and 
DAC's are usedfor memory. Thus, inparticular, there 
is a need for integrators with opposite or complemen- 
tary logic. These are called complementary integra- 
tors. They must be in the R(IC)-mode when the pro- 
blem is being solved and in the Hold-mode while the 
normal integrator is in the R-mode. The C-mode can 
be used instead of the H- mode if the integrator has no 
compute input, but only an initial condition input. 


The chart below shows the relationship between inte- 
grator modes 


Normal Complementary 
Integrator Integrator 
IC Cc 
H H 
C IC 


15-2 


The mode duty cycle for integrators in IO is 


Cc 
NORMAL {H 
Ic 


Cc 


COMPLEMENTARY 4 H 
ic 


The insertion of the hold interval avoids timing prob- 
lems between mode switches for normal and com- 
plementary integrators. Italso has another important 
purpose as shown below. Suppose a complementary 
integrator is tracking a problem variable, x, 


and x is rapidly changing. Then if the C-integrator 
goes to compute when the N-integrator goes to initial 
condition, the C-integrator will not store the final 
value of the N-integrator due to the time constant of 
the IC-circuit: 


N-OUTPUT 
C-OUTPUT 


N-MODE 
COMPUTE Ic 


This problem canbe solved if the C-mode command to 
the C-integrator is delayed: 


N- OUTPUT 







C-OUTPUT 
N-MODE 


COMPLEMENTARY 
C- COMMAND 


Without the insertion of the hold interval this cannot 
be accomplished. In some problems, both the N-, and 
C-integrators are used for computation, in sort of a 
bang-bang fashion. For example, ifthere are two sub- 
simulations, each of which depends on the final values 
of the other, then both types of integrators are used 
for equation solution, not merely memory. In this 


case delayed commands for both types are necessary. 
If all the variables are sufficiently slow, then the need 
for delayed commands is obviated. The implementation 
of mode duty cycles and delayed commands will be dis- 
cussed presently. 


The programmer needs aswitch at the console to start 
IO. This is the RO mode switch when the RO selector 
switch is in external. It follows that there must be a 
"safe'' mode when not in IO. This is the standby or 
S-mode and is activated by the reset switch for itera- 
tive integrators. Consider the circuit 


It is plain that the N-integrator must be in IC in the 
S-mode. The C-integrator cannot also be in IC in the 
S- mode because the two integrators would form a posi- 
tive feedback system resulting in instability. The C- 
integrator cannot be in compute since it may have an 
initial compute input and consequently would integrate 
to overload. Thus, the C-integrator must be in hold 
during the S-mode. Finally, there may be non-IO 
integrators in the program and these must remain in 


The patching diagram for IO integrator control is 


F F 
IRT el 
| 
R F 
RRO ROE PROX RO F 
RESET RESET 
R F R F 
o-— Loans o— — —o— —-o 


HOLD HOLD 


NON -IO WITHOUT 
HOLD 


NON-IO WITH HOLD 


The F-RT terminal is energized when the machine 
mode is Reset, putting N-integrators in IC and C-inte- 
grators inHold. The F-Holdterminal willbe energized 
when the Hold button is pushed, putting both N-integra- 
tors and C-integrators in Hold. There is no conflict 
withthe RO C-, R-terminals since these terminals are 
disconnected from the patchboard in any mode except 
RO. When the machine mode is RO, the problem Hold 
or H-terminal on the patchboard will be connected to 


Compute or Holdduring the iterative cycle. The fore- 
going can be summarized by 


NON-IO C 
INTEGRATOR 4H 
(NO HOLD) IC 


C 


NON-IO 
INTEGRATOR 
(WITH HOLD) jc 


N-INTEGRATOR 


| 
C-INTEGRATOR © ; 
: ; 
es 88 
oO 


MACHINE 
MODE R 


>t 


F F 
IRT RT 
RROW RO F R oe 
-=— = —-o 
RESET RESET 
R F R F 
o— i oN o- — — 
HOLD HOLD 
NORMAL IO COMPLEMENTARY IO 


the F-Hold terminal. Consequently, to put both the 
N- and C-integrators in Hold, it is necessary that 
the H-terminal be energized and that neither of the 
RO C-, R-terminals be energized by C, R inputs 
fromthe patchboard respectively. The diagram above 
shows patching for various integrators with mode 
relays. If electronic mode control is used, then the 
relays are replaced with electronic switches as shown 
in Chapter 6. 
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_ The logic program for basic IO control is shown below. 


DIVIDER 


AMAZCOO 
TE 
a 
=z 
oz) GS 


n 


TO RESET COILS 
- OF DELAYED 
COMPL. { 


YY 


C +V 
6 lhe |_| Las 
u et Ul} en 
ee err 
E ot 
HOLD 
Pe 
_[ 1 ° 
COMP 
ry = aw 
K-ts 5 5 a 
FS2 RESET 
4 | »s 
DOWN 4 UP | 7 
~ Ov F-RT' 


Fig. 1. Iterative Operation Logic Control 


Assume FS1 (Function Switch No. 1) is up and FS2 is 
down. Whenthe machine is inthe Reset mode, the out- 
put of the F-RT' terminal willbe alogical 1. This out- 
put is used to reset the divider and counters and to 
place both FF1 and FF2 in the R-state. The clock is 
off. FFland FF2formatwo-bit binary counter. (The 
IO cycle H, C, H, IC corresponds to the counter out- 
put 00, 01, 10,11.) Gate H is ready. Whenthe RO 
mode switch is activated the decade counters will start 
andbothN, C-integrators willbe inHold. At 10H, + Hy 
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the coincidence at Gate H will reset the decade counters 
and advance the binary counter to 01. (The use of the 
output of Gate DH will be discussed later.) The output 
of the binary counter will ready Gate C, thereby putting 
the N-integrators in Compute, and the C-integrators 
in IC by energizing the C-terminal through its driver. 


At 10C, + C, the coincidence at Gate C will reset the 


decade counter and advance the binary counter to 10. 
The output of the binary counter will ready gate 1 and 


put all integrators in Hold. At 10H, + H, the coinci- 


dence at Gate H will reset the decade counters and 
advance the binary counter to 11. The output of the 
binary counter will ready Gate R, put the C-integrators 
in Compute, and put the N-integrators in Reset by 
energizing the R-terminal through its driver. At 


10R, + R, the coincidence at Gate R will reset the 


decade counters and advance the binary counter to 00. 
The circuit is now ready togo through the same cycle 


again. 


If FS2is up, thenthe IOcycle will start with C (i.e., it 
will skip the first and only the first H interval). 


Delayed Hold commands are provided for the delayed 
complementary integrators by gate DH and the associ- 
ated FF. The logical equation for gate DH is, 


DH = H- DH, : DH 


1 2° 
Note that 


1< 10H, +H, 


10DH, + DH 
so that the FF will be toggled before the end of each 
Hold interval. This causes the Reset period of the de- 
layed complementary integrators to bracket the Com- 
pute interval of the normal integrators. 


Delayed Hold can be obtained by other means. Delayed 
Hold can be generated for both N-, and C-integrators 
by delaying the deactivation of the C, R inputs, res- 
pectively. The circuits for obtaining these delays are 
shown in the following diagrams. The inverters are 
those associated with the preceding gates. 


A simpler but less flexible IO logical control program 


DIVIDER 


1S 


COUNTER IN 
| | | 


FRT' 








If FS1 is down, the logic switch found on the counter 
module can be used to single-step the computer through 
the iterative program. 


If the machine Hold switch is depressed, all integra- 
tors, IO and non-IO, will go into Hold. This happens 
because the R, C inputs at the patchboard are discon- 
nectedfrom the F-RO and R-RO terminals inthe inte- 
grator modules and the main machine mode takes pre- 
cedence. Theclock input to the divider is also stopped 
because the RH terminal output is a logical 0. This 
saves the state of the divider, counters, and two-bit 
binary counter. Consequently, the iterative program 
can be put in Hold at any time while preserving its 
status. To return to the iterative program after being 
in Hold, it is only necessary to activate the RO mode 
switch. 


Sometimes it is desirable to control either the C or the 
IC-interval or both by an event in the program. This 
is accomplished by using comparator outputs instead 
of decade counter outputs as inputs to the appropriate 
gates. 





VY VV VY 


OV 
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This program will permit only an IO cycle composed 
of 10 time intervals. As before, the RH" output is 
used to gate the clock. 


MODE CONTROL NOTATION 


When delayed commands are to be used for an integra- 
tor it will be denoted by: 


[I> o | 


depending on its type 


Integrators can be controlled separately by pairs, not 
under general IO control. 


In this case the following notation is used: 


[> 


R=LOGICAL EXPRESSION (IC RELAY) 
H=LOGICAL EXPRESSION (HOLD RELAY) 


For example, R = L, . L,- The drivers are still used 


as logical andpower buffers. If a resistor rather than 
a capacitor is used for the feedback element of an inte- 
grator amplifier, then the mode control relays can 
still be used. The notation is 


>Papb«D 


R=LOGICAL EXPRESSION 
H=LOGICAL EXPRESSION 


Thus, a summercanbe made to act as a sampler ina 
sampled data system by controlling the hold relay. 


MEMORY CHAIN 


A memory chain program is 


f(t) 


- [i lp 


This circuit will store a sequence of values of f(t). The 
first integrator requires a delayed command to obtain 
the final value of f when the program which generates 
f goes into Hold. The chaincanbe any length desired. 
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ACCUMULATOR 


\W/ 
f(t) ( ) fer e 
| 
K 


This program will compute the sum of n successive 
values of f: 


_L 
fo” K 


Ms 


‘ f (t,). 


i 


The pot is for scaling. The output, ey asa function 
of time is 





Note that the complementary integrator has two IC in- 
puts. The second is obtained by patching a resistor 
between the output of the normal integrator and the IJ 
terminal of the C-integrator. 


AUTOMATIC RESCALING 


Sometimes it is necessarytorescale a problem during 
the computer solution in order to meet the accuracy 
requirement imposed on the simulation. This occurs 
because of the wide range of variables involved. One 
good example ofthis situation can be found in the sim- 
ulation of a nuclear reactor start-up, where the neu- 
tron level may rise through six decades of power. A 
technique is presented here for automatic rescaling. 
Although only a simple example is used for illustra- 
tion, the technique is applicable to the reactor start- 
up simulation and many other similar problems. 


Suppose it is necessary to generate a function, f(t), 
where 


- at 
f(t) = Ke 


a> o 


K = amplitude scale factor, 
n 


and as the function exponentially increases and reaches 
some limit, f(t) wax’ it is required to perform ampli- 


tude rescaling and proceed with the simulation. The 
following program automatically fulfills these require- 
ments. 





The integrator required for the generation of f(t) is a 
normalintegrator. The compute interval is controlled 
by a comparator. The function switch is in the Up 
position only at the start of the simulation and the ini- 
tial IC of the integrator is entered from potentiometer 
1. The switch could be replaced by a function relay 
energized by the FR bus. The IO logic program (Fig. 1 
or equivalent) goes through its cycle and waits in the 
Compute mode until an external signal commands it to 
go into the Hold mode. Note that the comparator out- 
put is used instead of the counter output to activate 
Gate C (Fig. 1). This signal occurs when the output 
of the integrator reaches a limit set on potentiometer 
2 and the comparator is triggered. A complementary 
integrator with a delayed command learns this value 
and feeds it to the IJ terminal of the normal integrator 
through potentiometer 3 and a resistor R. 


As the logic goes to the IC mode, a new initial condi- 
tion is established and the entire process is repeated. 


The amplitude scale factor K, is given by the follow- 


ing equation 


n 
Kk {0.058) n=0, 1, 2, 3---- 


where R is in megohms, 6 is the setting of potentio- 
meter 3 and n is the number of scale changes. 


FUNCTIONS OF A DEPENDENT VARIABLE 
This technique is not highly accurate, but is useful if 
precision is not a requirement. Suppose we want to 
generate 

f(x), x = x(t) 
If f is the solution of an ordinary differential equation 


L(D) f(t) = c(t) 


Where L(D) is a differential operator, then 


tx) - § L7 Dyedtyat. 


COMPUTE INTERVAL 
GATE OF LOGICAL 
CONTROL PROGRAM 







f (t) max 
100 


-100 

We see that time can be used as a dummy integration 
variable. Suppose g(x) is the solution of 

M(d) g(x) = b(t). 


That is 


x 
g(x) = \ om! (p)b(t)at. 


Then f {g(x)} can alsobegenerated, using time as a 
dummy variable, from the simultaneous solution of 


It 


g(x) 5 
f(g(x) ) { L”* (D)c(t)at 


eel 
| mt (ypet)at. 


g(x) 


The program for the generation of f(g(x) ) is, 


COMPUTE GATE 





|cd> |= f(g(x;)) 


The comparator is energized when eit) = x, so that 
t = g(x). Therefore, the output of the f circuit at this 
time will be 


f(g (xy). 








Example: f(x) = Ke 





COMPUTE GATE 


-!00 





The comparator will activate when t = x. The output 
of integrator N1 is Ke so that the output of integra- 


tor N2 is Ke *. If these IO integrators compute at 
high speed, say with 0.001 »fd capacitors and elect- 
ronic mode control, the output of N2 is a closed stair- 


Example: f(x) = x , x20 





case approximation to Ke *") for slowly varying x(t). 
This IO program can be used as a subroutine for a 
non-IO main program. 


The same general comments apply to the following 
examples. 


COMPUTE GATE 





15-8 


PERIODIC FUNCTION GENERATION 


Suppose a non-IO program requires a periodic func- 
tion. We can use an IO integrator and a function gen- 
erator to provide the periodic function. The function 
sub-program is 


t-RATE 
oo pn sie 


Here the N-integrator is under regular IOcontrol. The 
main program will see the output of the DFG as a per- 
iodic function because its non-IO integrators are in 
Hold during the time the N-integrator is in IC and Hold. 
In fact, the Hold can be eliminated if there is no other 
IO calculation to be performed. The non-IO integra- 
tor must have the logical control 


so that it will be in Hold, not IC, when the IO integrator 
is in reset. 


NEWTON'S INTERPOLATION FORMULA 


A 
f, = {(x,) 
f [X, | = f, 
f, - fie 4G. 
f Xo ma] F 
f |x,,x - f |x_,x 
1’°2 1] 
f [Xo ¥ 1 ¥2 | “ites he 
and generally 
f [Xor ¥4 ae x] 
f [1 Xp" ey x, | -f Xo Xr X,-1] 
eee 
n Oo 
We have 
f(x) =f [®o] + (x-x,)f [For] wooe- (1) 
f [Fo x)= f [x9 “l + (x-x,)f£[x, XX] o---- (2) 
and generally 


£[ XQ»: . Xp X]= f [*or" ‘ 2%, | + 


(x-x, )f [For "9X x] “ 


Substituting (2) in (1) 


eetht cae Gon 
+ (x,) (x-x,) f [x x14] 
and 
f(x) = f [*6] + Ox,)f lk, x, pee 
+ (xox Jt" Gx, 4)f es, | + E(x) 
where 


E(x) = (x-x,)° *s (x-x,)f lx," ti x x| 


If f has polynomial form, then for sufficiently large n, 
E=0. A good approximation for f, in any event, is ob- 
tained by assuming E=0. The result is called Newton's 
interpolation formula. This formula can be used to 
provide a continuous representation of a function de- 
fined at n points. The intervals Xi x,-1 need not be 


equal. However, to simplify our derivation we will 
assume that these intervals are equal. It should be 
remembered that the continuous memory to be devel- 
oped can be used for unequal sample intervals with dif- 
ferent values for coefficients in the computer pro- 
gram. We let Xoo o and assume all Xi-Xiy = a. 


Then 


(E- li fy 
f [*o ey] = , E fy e 


and 


a 6(E- 1)", m-l_ 


es IT es 
f(x) =1,+ > = jf (w-x,) 
m=1 mia 
so that 
x n (E-1)™4, nt : 
@-t+)) sare wh G-) ®) 
m=1 


For a specific number, k, of E"f, a representation of 
£€) can be obtained for (3). For n=3, we have 


(E- v", 


£) =f, + (E-1)f, ) + —z 9 & @-1) + 


(E- a 
——-@ C- 1) G- 2). 
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collecting terms 


- | (E-1)"4, aaa : 
£®)=£,+ |(E-1, - —g—"+—3—] @ 
2 3 _ 
. ey i 3(E-1) ‘| a2 
2 6 a 
3 
4 ee %3 
6 a 
1 3 11 
f£(x) = i. + Gf, - 5 fy + 3f, oe f,) x 
1 5 2 
+ (-5 fy + 2fo 9 fy + f/x 


+ @e,-42, + 42, -41)x 


( 3 
6°93 22 21 6°70 


_ 1, _3 
f(x) = £, + Gf, - {, + 3f 


11 
3 fo + 34, - Gf) * 


2 
x 
+ (-f, + 4f, - Sf, + 2f.) aT 


3 


x 
+ (f, = 3f + 3f, - f,) 31 


CONTINUOUS MEMORY 


In the analog computer we can retain a continuous rep- 
resentation, f, of f by sampling f at n points (not 
necessarily uniformly spaced). f can be generated 
later on command with the interpolationformula. The 
computer program which implements this operation 
from four uniformly sampled values of f is shown 
below. 





Comparators canbe used in conjunction with one ramp 
integrator to de-energize the R-relays of amplifiers 1, 
2, 3, 4 in the proper time sequence. The integrators 
are switched from the initial condition mode to the 
compute mode when it is desired to "playback" f(t). 
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If, instead of f(t), it is desired to playback f (T-t), the 
summers and switches are added to the integrator 
chain as shown below. 


tek |e fed <P lena || 


The procedure is the same as before except that the 
switches are energized at t=T. 


Although only the circuit for four uniform samples of 
f is shown, the technique can be extended to the case 
of n non-uniform samples. The pot coefficients can be 
calculated previously with a digital computer or desk 
calculator. 
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CHAPTER 16 
SAMPLED-DATA SYSTEM SIMULATION 


Sampled-data systems result from the utilization of 
telemetry, digital control, multiplexed controllers 
and other devices of a similar nature in a control 
system. The analog computer is a useful tool for the 
study, evaluation andsynthesis of both linear and non- 
linear sampled-data systems. 


THE SAMPLER 


The single-rate sampler is a switch whose periodic 
duty cycle is shown below. 


At 
CLOSED cg loa 
o-o— 
OPEN | | 
to tot T to t2T ! 


It f(t) is sampled by this switch, the output function 
is denoted f*(t). 


a ——— F*(t) 


f(t) #*(t) 





In most sampled-data systems, At is small com- 
pared to T. In classical analysis At is assumed to 
approach zero with finite area, so that f*(t) = f(0), 
f(T), ---. The introduction of this assumption may 
or may not significantly alter system behavior. When 
the analog computer is used for analysis, the sampler 
can be simulated exactly. The H-relays for the inte- 
grator amplifiers can be used to represent samplers. 
The duty cycle of these relays can be automatically 
implemented with the digital logic. 


The sampling operation produces sidebands. 
L oe] 


7 dL Fs +jnw,), 


n= -°% 


Lf* = 


(L indicates Laplace transform) 


where wo, = ar In order that the sidebands donot 


T 
overlap, we must have 


a ‘a 
Ff = 0, |w| >: 


(The operator F is the Fourier transform. ) 


If there is sideband overlap it is not possible (with 
techniques consideredhere) to recover f(t) from f*(t). 
Of course, in practice, it is usually not necessary to 
recover f(t) exactly. Consequently, judgment is re- 
quired to determine the extent of allowable sideband 
overlap. This is one question which is conveniently 
answered with the analog computer. 


FUNCTION RECOVERY 


An ideal filter for function recovery has the Fourier 
transform 


Te) = Awet) 


A(w) = 1, lw] sw, 
= 0, Jw]>w, 
gd(w) = -Kw. 


Thus, if F(w) is the Fourier spectrum of f(t), and 


Fw) = 0, |] > 


Cc 


then the Fourier transform of the recovered function, 
F'(¥), is 


F'() = T()F*(v) = F(u) e 


so that 
1 
' ss _ - 
(t)= 7 f (t - k). 
T(w) represents an ideal low-pass filter and a "cardi- 
nal data hold" ifk = 0. 


The impuse response of the ideal low-pass filter is 


w, Sin [Y. (k - t)} 


T(t) = T w, (k - t) 
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The filter is not physically realizable due to the pre- 
cursor (response for t<0). Ideally wy < w, < (w, - w) 


where W, > 2 Wy- shouldbe sufficiently greater than 


2 oa so that an adequate practical filter can be found. 


The foregoing applies to piecewise continuous func- 
tions. Special filter techniques such as clamping 
can be used for discontinuous functions 1. Hold cir- 
cuits are useful approximations of the ideal low-pass 
filter. 


Zero-order hold filter: 


Let the zero-order filter satisfy the relationship 


f(nT +t) = f(nT), O<t < T. 


That is, the hold circuit senses the value off at 

t = nT and holds this value until a new value is 

sensed at t=nT+T. Suppose f(t) = Up {t) so that 
1 

F(s) = 3 . 


Then 


If the transfer function for the zero-order hold cir- 
cuit is T(s), then 


T(s)F*(s) = + 
-sT 
= 1 1-e* 
Tis) = sF*(s) — s 
Thus 
-j T 
_, _ l-e ats 
a4 joT = =jwTl 
ag [Pe 
T(j#) = Te =< 
2G 
2 
1 


M. Schwartz, "Information Transmission, Modulation, 
and Noise, "' McGraw-Hill, 1959. 
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T(jo) = Te 7 = Gel) 
wy 
2. 
ta Sin cy 
G(w) = 
eh 
o) = -o 
G(w) 






CARDINAL HOLD 





W, 2We IW. 


A computer circuit which simulates a practical zero- 
order hold filter is shown below. Integrator modes 
are controlled by the IO logic. 


F*(s) 





-Ts 
I-e 
F*(s)| 5 


If an ideal hold filter is to be simulated, then the main 
computation must be interrupted while the memory 
loads. 

First-order hold filter: 


For the first-order hold filter, we have the predictor 
f(aT +t) = f(oT) += [£mT) - faT-TY] , O<t<T. 


If f(nT) = 1 forall n, then the output o(t) from 
the filter is 


o(t) = 1 + 


Thus, if the transfer function for the filter is T(s), 
then 


O(s) = T(s) F*(s) 
1 
Pitta pace 


os) = € + 4 )a- eS?) 
s 2 
Ts 
1 1 -sT 7 
~ Te) = @ +45) 0-08) 
Ts 
It can be shown that if T(jw) is expressed 
Tia) = Gel 
then 
Sin oe - 
G(w) = Not 41 2 
wT 
2 
-1 
go) = -wT + tan ~ w 


G(w) 






_—— CARDINAL FILTER 





Ws two S3ws 


K-th order hold filter: 


The predictor for the k-th order hold filter is 


: (™ m 
fmt +t) = So oS & f(nT) , o<t<nT 


m = 0 


where A is an operator such that 


Af(nT) = = f(nT - T). 


Practical low-pass filters: 


Recovery can also be accomplished with a variety of 
low-pass filters. The analog simulation of the sampled- 
data system usually will provide a means of deciding 
the filtering requirements. Various types of filters 
may be tried in the simulation. 


RIPPLE 


Various responses in the sampled-data system may 
contain components that are either completely or 
partially hidden by the sampling process. These are 
called ripple and may have a bearing on system 
Stability. In classical analysis, the system response 
between sampling instants is examined either by the 


introduction of fictitious samplers or the use of the 
modified Z-transform. The analog computer is an 
effective tool for evaluation of ripple since the con- 
tinuous time response of all elements of the system 
is available for observation. 


PULSE TRANSFER FUNCTION 


In what follows, it will be assumed that all samplers 
are synchronized and have the same period T. In 
order to find the pulse transfer function of an element 
of a sampled-data system, consider the simple sys- 
tem below 


une 





We have 
C(s) = G(s) R*(s) 
C*(s) = # > C(s + jnw_) 
C*(s) = # > G(s + jnw,) R* (s + jnw,) 


=-0 


Note that R*(s + jnw,) is periodic with period W.- 


(Consider the sideband structure. ) 


Thus 


i} 


R*(s + jnw,) R*(s), for alln 


C#(s) = R*(s) * ye G(s + ine)... Equ. 1 
Nn=- © 
C*(s) = G*(s) R*(s) 


From the relationship between the Laplace and Z- 
transforms, it follows 


G(z) = oe 


Where G(z) is the pulse transfer function and is the 
Z-transform of g(t), the impulse response. For con- 
venience, we will use F(s), F(z) to denote the Laplace 
and Z-transforms of f(t), respectively, although F(s) 
and F(z) have different functional form in general. 


PULSE TRANSFER FUNCTION OF VARIOUS LOOPS 


For convenience, we define 


z{F(s)} 4 z{ulFG)} 
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The system 
rt) C (t) % 

()}— a > C1) 
satisfies 

C(s) = G, (s)G,(s)R*(s) 

C*(s) = [G, 6)G,(6)R*(s)] rs 

C*(s) = [S,6)G,(6)] *R*(s) 
Example: 

r(t) 

Example: 






(see equation (1) above). Thus 


GZ} - Z {G,(s)G,(s)} 


It is convenient to denote Z {G,(s)G,(s)} by G,G,(2). 
It is important to note that G,G,(z) # G, (2)G, (2). 


The computer block programs for the simulation of 
various simple sampled-data systems together with 
the closed-loop pulse transfer function are shown in 
the examples below. In each case the sampling func- 
tion is implemented with mode relays under IO logic 


1 1 


control. In allcases G = aa H = a 


(The symbol f for an integrator means it is not under 
10 control. The normal IO logical program is used 
for the sampling amplifiers). 


—e cit) 
C(z)=GR(z) 


c*it) 


—e ct) 


G(z)R(z) 


C(z)= 7 GH(z) 
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Example: 





G(z)R(z) 
I+G(z)H(z) 





C(z)= 


r(t) 





Example: 


r(t) —e"(t) 


GR(z ) 
I+GH(z) 


C(z)= 








DESIGN TECHNIQUES 


Root-locus: 


Root-locus plotting techniques apply without modifi- 
cation in the z-plane. The region of stability is 
|z| <1. The interpretation of the locus is somewhat 
different. The family of circles. 


represent constant damping, a, (not damping ratio) 
for the response 


f* = [eno coat] * , 
The family of radii 
Argz = @ = oT 


represent constant frequency for the above response. 


Thus, ¢, the damping ratio is given by 


¢ = Sin 


or equivalently 


t = sin”! 


1 a 
VJ 2 2 
Q@ + W 
oO 
V 2 2 
gc +0 
Pole-zero cancellation techniques: 


Suppose the desired pulse transfer function for the 
controller 
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(where G(s) is the compensating function) is given 


by T(z). Then since 


G _G(z) 
ZI) nO 
T@) = REY 7 +G,G@) 
we have 
T 
GG(z) = 7S, . 


This technique is useful when G. is realizable. 


Digital Compensation: 


A digital controller can be used instead of a com- 
pensation network above. This controller can be 
Simulated with the analog computer if it is suffi- 
ciently simple. If not, then a hybrid system may be 
used with the analog computer simulating the remain- 
ing sampled-data system. 


MULTIRATE-ASYNCHRONOUS SYSTEMS 


Multirate and asynchronous sampled-data systems are 
at best difficult to analyze with classical techniques. 
The analog computer, with its asynchronous control 
structure finds powerful application for such problems. 


CHAPTER 17 
PARTIAL DIFFERENTIAL EQUATIONS — PART 2 


PROPAGATION OPERATOR 


There isa class of equations which can be solved with 
the similarity transformation 


i 
o(R) = (TAA T) 6 (eR) 


where T is a suitably defined integral transform over 
8R (boundary of the region, R) and propagation into R 


is along the x’ coordinate of R. That is, the surfaces 


x! — constant are "concentric" to aR. (T" 1, Xp) is called 
the propagation operator and A the propagation function. 


That TL a*T exists, canbe shownif R is homogeneous 
in the boundary coordinates. That is, usually R can 


be non-homogeneous only in x’. We will use this ap- 
proachtofind Green's functions for the solution of dif- 
ferential equations where we know the required inte- 
gral transforms over oR. 


First we derive Afor a few simple cases for the pur- 
pose of illustration. 


Consider the cases: 


yy xx (1) 
Wy = Uy (2) 
te Ux (3) 


We willuse the Fourier transform and its inverse. We 
claim (for propagation along the x coordinate) 


uy) = TAXG («) 
= = J SG (wel? Faw j 
Then, 
ee = 7 f se Gwe! ¥ dw 
Me > f fel K G(w)elY dw 


In A A~ Gwe! dw 


ro 
Md 

li 
role 
q 
8mm 8 


For (1) above 


1 ¢, 2.2 
Woy + x =o Seo +In“A). 


& e!Y Glu) dw = 0 


for arbitraryG. Asufficient condition for the vanish- 
ing of the integral is 


iIn?A = an 
A = e | » 
and 
u(x,y) = pl e 11% Gy 
-1 
u(0,y) = T ~ G(w) 
so that 
G(w) = (T)u(0,y). 
Consequently 
_ mel = |w|x 
u(x,y) = (Te T) u(0, y) 


and (rte lel %p) is the desired propagation operator. 


For (2) we have 


wy 


G(’) dw = 0 


e j 
s. few - Inde 


so that for arbitrary G 


Aze” 
and the operator is 
ql ee” 
For (3) we have 
2 
ey es 
In A = 











Let 


- -1 
uy) = Tat G, () + TT? Ay Ge) 
_ cl eG, (u) ¢ ri eo JOx G2 (a) 
Then 
jwx 1... -jwx 


-1. J 
u, (x, y) = T “jwe G, (+) - T “jee Gy (a) 


u (0,y) = co} G,(%) aol Gy (a) 


oe: Se 
u, (0, y) = T jG, (») A jw Ge (w) 


G, + Gy (w) is (T) u (0, y) 


G,(a) - G4) = — (1) uy) 


Gg) =F [Mu@y) + Gu, y)] 


Gy) = = [ux @y) - 757) u,0,y)] 


-1 el? 
u(x, y) (x : oe cr) u(0, y) 


2)" 
‘ ( ce 5) u,.(0, y) 


eo 
=F m) u(0, y) 





2 (ot 0 = u, (0, y) 
u(x, y) = Ce coswx m) u(0, y) 
+ (x an (x) u, (0, y) . 


Let's examine one more case 





2, (EI-3) + m 24 - 0. 
ay ay ox 


Let EI, m be constants and define 


m 
k= EI 


so that 


then 





and we see from the last case 


2 
u(x,y) = a cos(— Te) u(0, y) 


k 


w 
—f$—] cn) u,_(0, y) 


SUMMARY FOR POLAR COORD'S (r, 9). 


an aig 
(v)£(e—) = f tthe ae = Fa) 
0 
qo} F(n) = a F(n)e!"® = f(0) 
For v6 = 0 we have 


u(r, 6) = [a r! I ¢r)) u(0, 6) 


For v4 + "4 = 0 we have 
u(r,e) = ( (> cc) sar }t0) u(0, 9) 


HIGHER DIMENSIONS 


We can derive similar results in higher dimension with 
the use of suitably defined transforms on AR. 


GREEN'S FUNCTIONS 


We proceed to use the propagation operator to obtain 
Green's functions. 


Casel: u + u = 0. 
See es yy 


u(x, y) = (r7,-I917 2) u(0, y) 


L c ~|wlx jw 
= L (fe e ¥ ao): 


=o 


oo 


§ u(0, 2) eI ay 


= 0 


= a J u(0, t) as) f e [ox joly- 5) dy 
Now 


f e 141% = IAY- Sa, = 


= oO 


f i” (y-¢) eX gy 
0 


0 
= foes lty-5) a, 


ea oe ee eR 
x+jlo-y) — #HE-y) Pay yy? 
so that 
uy) == f BO) ay 
-2 x +(y-f) 


Case II: u =u 
—— ae x 


0 2, 
u(x, y) = = fe® XIV a, 


- Oo 


-( f u(0, te? Sar 


™ oo 


== fulo, daz 


oe. Oe 9 
fe x JO(- Da, : 


-— 0 


Now 


20 2 
f e © Xcosw(y-t) dw 


=o 


o 2, 
fe? x Jjoly-2) _ 
2 2. 
+ j fe ® * sinw(y-t) dw 


2 2 
= 2 fe® Xcosw(y-6) dw 
0 





2 
-(y-&) 
2V1rx 
so that 
ss ys) 
u(x, y) = fe * wOs) a 
Wr -2 Vx 
Case lll: u =u 
> yy xX 


u(x, y) = (oe coswx 7] u(0, y) 
-1 sinwx 
+ (7 ar 7) u,(0, y) 
Let f(x,y) = (ie coswx 7] u(0, y) 


3 -1 sinwx 
gis y) = (vt x 7) @y) 


Then 
1 r jwy 
f= or S coswx e dw 
J woo, yeas 
1 eo 
ae Jf 00, ae 
f coswx lo 6) dw 
Now 


feosex Jay- Dg. = f COSWX COSw(y-%) dw 


=a OD — OO 


foosex Joya, = Il E (x-y+0)+5 (-x-y+2) | 
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and 
1 oO 
ty) =F ful0,c) [(&-yes) 
+ a(-x-y+6) dg 
f(x, y) -5 [w(o, y-x) + u(0, y+x)| 
Also 
gonad, fenes oY as 
e -jwt 
f u,(0, £) e dg 
-+ f u,(0, o)d¢ JS — [cosuty-x) 
+ j sin wly-t)| dw 
for y-¢>0: 
[SINE cosuty-4) do = ry y= <% 
= ar, y-G=x 
7 0, y-G>xXx 
for y-¢<0: 
{sinux cos(y-4)dw = 7, -y<x 


=xeQr, C-Y=xX 
=0, C-yr>x 


for the first case (y-¢>0) 

(==, y>try-x 
for the second case (y-¢<0) 

( = 7, ytx>O>y 
so in order to meet both conditions 

\- ™, y+x>O>y-x 


and 


ytx 


ay) = =| 


u_(0, c)d¢ . 
y-x * 
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Finally, 


u(x, y) = 5 [w0o, y-x)+u(0, y+x)| 


y+x 
u (0, t)dg . 


i 
2 | 
y-xX 


Case IV: 
Consider Uy = u, inthe region0 =y =7, 0=x. 


wv 


TE(y) =§ £(¢) sin nea = FC) 
0 
a F(n) = ae F(n) sin ny 
nr n=l 
-1 x 2° X,. 
u(x,y) = T A Gi) == 2A (sin ny) G(n) 
_2 2 .X& 2 7; 
wy == 214 (-n“) (sin ny) G(n) 
u =2 3, A* loga (sin ny) G(n) 
x «nel 8 
We have 
logA = -n2 
ane 
A =e : 
-1 “pee 
u(x,y) =T e G(n) 


u0,y) = T Gia) 


G(n) = (T) u (0, y) 


-1 


2 
ux,y) = (Te *T)u(0,y). 


-n x.. 
e sin ny 


a) 
Me 


u(x, y) = 
1 


u(0, ¢) (sin ng) d¢ 


Cay F 


u(0, 6)d¢ 


i 
a |p 
ons 


‘60 2 
>> e ™ * (sin ny) sin ng 
ned 


or to write out a few terms 





T 
u(x, y) -2 e “siny f u(0, £) (sing)dy +e *¥sin 2y 
0 
T 
(f u(0, £) (sin 2a 
0 
-9x . - . 
+ e sin 3yf u(0, ¢) (sin 3g)dg+... 
0 


All of these techniques can be implemented using re- 
petitive operation. The idea is to evaluate at high speed 
the equations for fixed values for the coordinate para- 
meters. Then by changing these parameters slowly 
the solutions can be found over the region, R, of def- 
inition for the equation. 


The program is: 
u(o,ytt) 


7X -) 
O : me 


a_i 
100 
-100 : 
WF = 
| ZN +2 
1 dt 

O) 100 dr \/ 

-100 +100 FR-B 

FR-A = 





A 


P COUNTER 


RESET OV 


EXAMPLES 


u + 
xX 


Lett = t - y 


be: x-ulo,yft) 
aS 


R H 


/\ ZN QA 





o FF MASTER RESET 


Ol 
/ 


u = 0 
yy 
1 oo 
u(x, y) = rs xu) | dg 
-*© x°+(y-&) 
1 7 
u(x, y) = = f xu(0, =D) dt 


1 f xu(0, yt 
u(x, y) =— fA ptt) dt 
0 


u( 
xe 


+f xu(0, y-t) dt 
o x +t 


0 u(x,y) 


Ro 











The FF's act as an up-counter with successive trigger 
inputs to FF1. When the machine is in reset the counter 
output is zero and the integrators are in the standby 
mode. (The ROselector switch is in external). When 
the RO pushbutton is activated, the N-integrators go 
tocompute. As the counter is incremented, the mode 
sequence is 


COUNTER MODE 
3 2 1 FF 
0 0 0 
001 H 
0 1 0 IC 
011 H 
100 
101 H 


Note that (FF3) = 0 when the counter has reached 101. 
When this occurs there are no more trigger inputs to 
FF1 and the program stops. Before this, the trigger 
for FF1 is provided by the clock (Lg) when the machine 
is not in compute and by the comparator (Cj) when the 
machine is in compute. There are only two compute 
intervals. The function relay (FR) is energized (+) 
during the first compute interval and de-energized (-) 
during the second compute interval. 


The analog program is 


-l00 


rome 


-100 


| 
ax 


y 
400 i Y 
om) |N> ae = FR-A > 


Fp 


The first and second integrals are calculated during 
the first and second compute intervals, respectively. 
The comparator is energized when y +t reaches its 
maximum absolute value for each calculation. It then 
terminates the compute interval. The final output from 
the accumulator is the function u(x, y) for the partic- 
ular x, y chosen for the run. 








Uy = ve 
2 v sy 

u(x, y) 5 fe u(0, g) 9) dt 

a 2 

Let t=¢- y 

peste 

u(x, y) = — fe4* uo, z) ae 
4n Nix i 

-t? 

u(x, y) = J f ere u(0, y+t)dt 


4n Nix 0 


-t? 


+f eM* u(0, y-t) dt 
0 


u(x,y) 


SI 


GATE | 


Both the logic program and the principle of operation 
ofthe complete program are the same as for the pre- 


vious example. 
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CHAPTER 18 
CORRELATION ANALYSIS 


PROBABILITY THEORY 
Sample space: 


Asample space, S, is the collection of all outcomes of 
an experiment. 


Probability: 


P(S) is a function defined over S. P maps S into the 
interval 0 < P(S) < 1. Pis not necessarily univalent. 
The integral (or sum) over the range of P isunity. 
In case § is finite, P is said to assign a probability to 
each point of S. In case S is infinite, P is said to 
define a probability density over S. Suppose Sq isa 
nested collection of subsets of S such that 8; C S;, 1 


for alli. The extremum of Sg= S and the infemum of 
S,is empty. The measure of the range of P for the 
monotone sequence Sq is the probability distribution 
for S with respect to the collection Sq. 


EXPECTATION AND VARIANCE 


The expectation of S, E(S), is the measure of SP(S). 
The variance, V(S); is the measure of 


s*pis) - E(s) 
Example: 


Suppose S is given by -~ < x < ~ and 


3 (x - y)" 
P(x) = (27) exp =. 
Then 
E(x) = [= P@) ax = p 
V(x) = J x*Pewdx-1? = a” 


—_ oO 


Let the sy be defined by -~< x < Xv? then 


x 
ox,) = j P(x)dx 


~oO 


is the distribution function. 


STATIONARY RANDOM PROCESSES 


If the statistics of a process are time-independent, 
the process. is said to be stationary. 


Suppose f(t) is a random function defined for 
-~e <t< », 


Let the domain of f be denoted by D. We define an 
ordered countable set {D_} so that 


Do: x, =t< x, x 


, - xX. >A. 
i itl’ “i+l i 


UD. = D 
i 


If f, = {(D,), then the collection i is an ensemble. If 


f is stationary and A is sufficiently large, the ensem- 
ble will be ergodic. (That is, the statistical properties 
of the members of the ensemble are identical.) Hence- 
forth, all ensembles considered will be assumed er- 
godic. 


FIRST AND SECOND PROBABILITY DENSITIES 
Let S be the sample space consisting of all possible 
values for the ensemble f(t, + A) where A<A. 


P(S) is said to be the first probability density. P(S) 
is independent of A, since fis stationary. The ensem- 
ble is ergodic so that 


6-2 


6 
E(s) - Lim xe io 
6 


Lim 


6—~ 2 


|e 


v(S) 


i) 
oO 


6 2 
2 

( £7 (t)dt - E(S) 

6 


Let Sp be the sample space consisting of all possible 
pairs of values (f(t; + A), f. (t, +4 +T)) for the en- 
semble where A+ T < A. PS) is said to bethe 
second probability density. PS) depends only on 


T, not on A. The ensemble is ergodic so that 


Lim 


1 
E,pSp) = 2m@ x; | swt + That. 


Or ee, Oy 


(It can be shown that this equation is valid for T > A.) 
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AUTOCORRELATION FUNCTION 


ES); above, is the autocorrelation function for f 


and is denoted by dee(T). That is, 


6 
_ Lim 1 
p(T) = Lim { sete + Tat 
-6 


The autocorrelation function has the properties: 
2) Pee (o) is the average power of f. 
3) |oe(T)] < a4 (0) 
Proof: 


Lim 


léee(T)| <5 ag J [te] - [een] at 


ees, 


a 


Lim 


= 6—~ 0 


|- 


6 
(4 jl? + 


~6 


bh 
oO 


lee +7)| 2 [few] - jsc+7)[] 24 ae 


6 
< o,,(0) % ae e ( I f(t)| 7 
-6 


| e¢t +7)|| * at 


4) fe Pi ——> ep € P,, (is periodic with period ») 
Proof: f(t + T) = f(t+T +») 
5) Pee contains no phase information about f. 


6) Forevery f there exists a unique Pep? the con- 
verse is not true. 


7) Iff & P, for any », then Lim p(T ) = 0 
Too 


8) If f can be represented by 


f = >. A Sin(wt + ¢,), 


n=-0 


then 


Pe (T) = » Pan 


n=-° 
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6 
where ¢ , =Lim a ( A” sin (wt £ ¢,) 
67> 0 
-6 


: (sm (wt +T + $,) at) 


Proof: 
6 
Tia! ae J ajsin(w,t + ¢,)). 
Bice 26 n n n 
-6 


A Sino, t +T + $,,)at = 0, m=n. 


rt 
9) If f, f'< bd, then dee (T) = - yp (T)- 


Proof: 
6 
Gergs (T) = Lim F ( f' (t)£" (t + T)dt 
67> © 
-6 


67> 


5 
= Lim 5 [sme + 7)] - 
-5 


6 


fs, (t)e"'(t + T)dt 
-6 


1 
= -— =z Lim 3; ’ 
aT? gra 2 


j f, (t)fg(t + Tat 
-3 


CROSS-CORRELATION FUNCTION 


Suppose f, g are two stationary random functions. f, 
g can be decomposed (as above) into two ergodic en- 
sembles f ow 8a’ The union of these ensembles is a 
Let Sr be the 
sample space consisting of all possible pairs of values 
(f(t, + A), gi (t; + A + T))for the composite ensemble, 


composite ensemble of pairs (f, g),° 


where A +T <A (A is the ensemble duration). Since 
the ensemble is ergodic 


E,(S,p) = Lim 


67 0 


6 
5 J f(t)e(t + T)at 


En is the cross-correlation function for the pair (f, g) 
and is denoted by 5(T). 


It can be shown that 9(T) is valid for T >’. The 


cross-correlation function has the properties: 


1) Generally o_(T) # %_(-T) ; 


Consider g(t) ,-extc¢e 


= 0, |t| > 


2) d(T) = dy¢(-T) 


3) Max Op (T) does not necessarily occur for T=0. 


4) | ,(T)| < ¥ %¢(0) $600) 


Lim x ( o7 (teT) dt 
6 
600 a 

AUTOCORRELATION FOR A SUM OF FUNCTIONS 
Henceforth we will denote 

f(t)dt 


ai 1 
Lim 25 


Qe, ©, 


670 


1 
m 


by Mean {£(t)} , or just M{f} wherever convenient. Let 
f(t) = f, (t) + f, (t). 


beg(T) = Mean {[ 4, ) +f, (¢)] [#,(¢+7) +f, +7)] 


= M {£,(t)f,(¢+T)} +M {£,f,(¢+7)} 
+M {ret (t+7)} +M {in(t)f, (t+T)} 
Thus, 


Pee = de. + ¢ + @ + @ 
ff ff, ff, fof, fofo 


The general rule for a linear combination can be de- 
duced easily. 


EXAMPLES OF AUTOCORRELATION FUNCTIONS 


Example 1 (sinewave): 


Let f = ASin (wt + 2), then 


Pee = Mean { a? sin( wt + ~) Sin( wt + wT +y)} 


Since f ¢ Py» 





an 
2 ato 
wA / ; 
dep = “Oe Sin (wt +y) Sin(wt+%T +y) dt 


a 


aw + yy + An 


a2 
Pee = o> Sin(¢) Sin (4 + wT) df 
aw + 
a2 Qn a2 
es = Ap j COSwT dt= > «cos wT 
oO 


Example 2 (general random function): 


Suppose f is defined by 
a) fm6 +t) = f(nd), o< t < 6 
b) P [f(ns) = x] = P(x) 


c) E [f(m6)] = 0 


then 


2 


E [£(c)e(e + T)| = f [xP(x)] | xt E| Jax 


=—0o 


oo 


\ [xP(x)] LyP()] | dxdy, |T|< 6 


—-o 


= ( [ xP(x) | [yP(y) | dxdy, [| > 6. 


-%O 


But \ xyP(x)P(y) = 0, so that 


—- 3 
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Example 3 (binary noise): 


T 2 T 
Pep = a-|2| ) | x P(x)dx = V(x) (1- | ), TS6 Binary noise is a special case of (2), above, where 
= P(x) =2 u Gey ek (x +x_) 
= 0, T> 6 2 °1 Oo 21 an 
That is f takes on values of x_, -x_ with probabilities 
Pet 1 o” “*o 
of 3° 
V(x) 
Then V(x) = xe and 
2 T 
T zs ps eee 
8 3 tee = UI z| 
Example 4 (periodic binary noise): 
Suppose f is defined as in (3), above, except that 
f(t) = f(€+v) for large v. Then Pee is periodic. 
> ff 
x8 
-v -8 o 8 v 
POWER SPECTRAL DENSITY r 2 
: 1 -j 
The operator F denotes the Fourier transform, so that & (%) = Lim aa | g(t)e Jw at 
gg r 
Fg = G(»). oe a 
jor aa is said to be the power spectral density. We 
$0) = Foo e= a j e dT have 
-6 
i f¢ tgg(T) = \ ge" oe 
Lim 5, j g(t)g(t +T) dt -% 
6-0 “y 
r WHITE NOISE 
6 (% = Lim x j g(t)dt White noise has constant power spectral density of ar- 
Be 6+ LA bitrary magnitude (say A). Thus, the autocorrelation 
function for white noise is 
6 
: = Au,(T), 
Lim { g(t+T)eI?t dT "ge 1 
Ores see where u, is the unit impulse. 
i A i. PERIODIC WHITE NOISE 
6 (4) = Lim 55 j g(t) e! ‘dt oe, . 
88 §-00 A If f satisfies the requirements for white noise for some 
-r interval 6, where 6 is large, and if f(t+é) =f(t), then 
6 eo 
; -jot 
Lim e d 
vem g(c) C bs(T) = ay) u, (T-n6). 


n=-«0 


FINITE CORRELATION FUNCTION 


If f, g are defined only foro =t< ©, then 


6 
Soper oa 
Pe = oe ‘ Jr) g(t+T)dt. 


YF can be approximated by a sum: 


N 


= 


et f(m At) g(m At + T) (1) 


%p = Mean 


Inthis case At mustbe small enough that the high 
frequency components off, g are adequately represent- 
ed (see Chapter 16). N must be large enough that the 
mean of the sum is a good approximation to PE" 


CORRELATION PROGRAMS 


Simple computer programs for (1) require that 
At > T. (2) 


To obtain the best approximation, At = T+ ¢, where 
e« is the minimum time required for sampling in the 
computer program. (2) imposes a limitation on the 
magnitude of T. Thatis if T is large, then so is At. 
Thus, for large T, relatively few samples are made 
of f, g. This restriction is not unduly severe, because 
ty cannot be determined accurately for large T to- 


gether with finite histories for f, g. 


If (1) is to be implemented directly, the program 
is 





Integrator 1 samples f and then after an interval, T, 
integrator 2 samples g. Integrator 3 is in hold until 
both samples are made and then it samples the product 
f(t,) g(t, + T). The output of integrator 3 is averaged 


by a Mean Program (Chapter 3) which inturn generates 
PE" If T is to be specified accurately, the integrator 


capacitors must be as small as possible to minimize 
the effect of the IC time constant. 


The total sampling time can be reduced by replacing 
integrator 1 with a quantizer whose input is ftogether 
with dither (see reference). f will be represented by 
a two-bit word (output of the quantizer). 


The table below defines the quantizing function 


BINARY DECIMAL 
OUTPUT OUTPUT OUTPUT 
~-2=x<-1 . 00 -2 
-l=x< 0 01 -1 
O=x<l 10 +1 
l=x=2 11 +2 








Except for scaling the program is 







ES3 





MEAN 








ES2 
PROGRAM $ fg 
g 
R 
f 
permit DITHER 
PROGRAM 
rae ; 
x00 ~100 mit i Ct 
RTS RTS RTS 
ES3 ES2 sal 


ov a 


Ss) oivineR HE} counter 
F-RT’ +V 


i 
ES4 


: R 
ES MEANS 
ELECTRONIC SWITCH COUNTER 


rife 
COUNTER 


The Dither Program is found in Chapter 10 and the Mean Program in Chapter 3. 


Thecorrelation parameter, T, and the reset interval Reference: Korn, G.A., "Hybrid-Computer Tech- 
niques for Measuring Statistics from 


for integrator 1 are given by (10T, + 10T, + T,) P 
and (10 I, + 101, + 1,)P, respectively, where P is Quantized Data’) Sinulation, dspei tse: 


the period of the output of the divider. 


CHAPTER 19 
ON-LINE DATA ANALYSIS PROGRAMS 


An analog computer can be used to advantage for on- 
line data analysis. Programs are given in this section 
which correspond to the more frequent operations on 
incoming data. The analog computer is used for both 
processing (analysis) and editing. The principles de- 
veloped in other chapters can be used for processing. 
The programs here are for simple editing. 


The program which measures error frequency is 


ERROR PEAK VALUE 


Suppose, with the above definition of error, it is re- 
quired to find a sequence of peak values of error, 


Max |f-¢| 
Max |t-¢| : 


The program is 


ABSOLUTE 
VALUE 
DERIVATIVE 












MEMORY 
CHAIN 





ERROR FREQUENCY 


Suppose an error is defined to occur whenever 


[#t) - g(t)| > «- 






COUNTER 
COUNTER 
COUNTER 






NEGATIVE 
ABSOLUTE 
VALUE 
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All the programs indicated will be found in various 

_ chapters ofthis manual. The integrators in this pro- 
gram are under iterative controlfromthe C, R logical 
controlinputs. Pot 1 isset toa small value to provide 
positive zero derivative detection. 


TIME BETWEEN EVENTS 


Suppose two events, E}> Eo are defined by 


f(t) - g(t) 2 0 


Ey: h(t) - m(t) = €9 


h(t) - m(t) 2 


0 


A program which will determine a sequence of time 
lapses between E,, Ey (not Ey, E,) is 






COUNTER 
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CHAPTER 20 
SYSTEM OPTIMIZATION 


System optimization is a procedure for generating a 
response, yg, of a known experimental system so that 
it matches, as closely as possible, some desired res- 
ponse yg. yq maycomefrom tape or may be the res- 
ponse of a known but more complex system. In what 
follows we will assume that either an external device 
generates yy or else we know the computer program 
for it. A measure of the accuracy of fit of yg needs to 
be established. This is called the Criterion Function 
(CF). The CF has different forms depending on whether 
dynamic or static optimizationis required. The exact 
form of the CF in either case is a matter of technical 
judgment and will depend on the particular situation. 
The CF for the dynamic problem will need to take ac- 
count of the transient behavior of both yg and yg over 
some intervalof time. Typical choices for the CF are 


STATIC 


ple 2 
RS a. = Ya) 


DYNAMIC 
f = (" [y,() Z yq(t)] * at. 


(Note that the CF is defined so that it is non-negative) 


The optimization problem amounts to finding maxima 
or minimaoftheCF. Since yg is afunction of its para- 
meters, say Cj,..., Cy, so is the CF. Thus the CF 
can be viewed as a surface over the parameter plane. 
In general the C F will have several maxima and minima. 
We will restrict our attention to finding the nearest min- 
imum from some initial point in the parameter plane. 
Once techniques are developed for the minimum prob- 
lem, they canbe applied with modifications, which will 
be obvious, to the maximum problem. 


We will develop two methods for finding a minimum: 
steepest descent and relaxation. The static and dyna- 
mic optimization problems will betreated separately. 


STATIC PROBLEM 


The family of curves f = constant onthe surface for f, 
the CF, are constant elevation contour lines. Differ- 
entiating 


$f ac =vt. d= 0 
1 2C, j 


df = 


TMs 


where C is the parameter vector -- i.e., 
C = (Cy, Cg,..., Cy] (Remember that the CF is a 


function of the parameters of yg. ) 


We see that Vf is perpendicular to dC and dC is adif- 
ferential vector which is tangential to the space curve 
f = constant. The projection of the curve on the para- 
meter plane is identical geometrically to the space 
curve because it is at constant elevation over the plane. 
The projection of Vf on the parameter plane is perpen- 
dicular to the curve. It is evident that Vf corresponds 
tothe slope of the surfacefandVf = 0 atthe minimum. 
Toget tothe minimum from some point on the surface, 
then, all we need to do is follow - Vf. (If we wanted to 
find the maximum we would follow + Vf.) The pro- 
jection of the path on the parameter plane is 


STARTING 
POINT 





XS f= CONSTANT (TYP) 


If the system is static (algebraic) we can change the 
parameters of y, while following - Vfwithout intro- 
ducing transients in the system. Now if 


2 
f=-5 


Ss Yq) 
then 
vi = (y, - ¥g) V¥,- 


If we lete = yg - yg, then C can be determined by 
(since we follow - Vf) 


C = -e j vy, dt. 


or in scalar form 


oy 
s 
C, =e \se- dt. 
dy 
es oe le 
C.. e Ise dt 








In practice, itis usually advisable touse gain factors, 
K, so that 








dy 
he ks s 
Cy = Ke 5G dt 
dy 
aes Ss 
C. Ke Ise" dt. 


K is chosen sothat the procedure converges but at the 
same time is stable. When Vf = 0, the integrals will 
be stationary andC will be the parameter vector (point 
on the parameter plane) corresponding tothe minimum 
of f. Vyg is generated by differentiating the generating 
equation for Ye with respect to the parameters C.- 


DYNAMIC PROBLEM 


Except under specialcircumstances, Vfcannot be gen- 
erated on aninstantaneous basis. Even if it could, the 
parameters could not be changed during the determina- 
tion of f. Sincethe calculation of f is carried out over 
a finite time interval, varying the parameters would 
lead to a result which is not the CF. 


Some definitions for f will permit Vfto be computed as 
oftenasis f. For example, if 


1 -T 2 
f= (y. - Yq) at 
2 ‘5 Ss d 
then 
_ eT 
vi = i (y, - Yq) Vy,at. 


As soonas we knowthe components of yg, we can use 
a steepest-descent method. After each solution a new 
position in the parameter plane is chosen by taking a 
finite step along -Vf. Now Vy, canbe generated by the 
parameter- influence equation method. The differential 
equation for each component of Vy. is obtained by dif- 
ferentiation of the generating equation for y,. For 
example, if 
C, Vg + Cy Vg = 0, y,(0) = 0, y,,(0) = A, 
Oye OY, 
the equations for ac,’ aC, 
ferentiating first with respect to Cy and then Co- 





can be obtained by dif- 








2 oy dy 
5¢ 7 Cy SS Ge) + Cy a ~ “Ye: 
1 dt 1 1 
dt aC, , aC, 
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a, d s s |. 
3c, * © (Fe) + ©2 Er = 7? 
2 2 2 2 
dt 
a (“¥s(0)) _ 9, 280) « g 
dt. 9G. aCe 


These parameter-influence equations can be solved 
with an analog program. Their non-homogeneous terms 
are obtained from outputs of the generating program 
for yg. It will be noticed that they have identical form 
except for the non-homogeneous terms. Due to equip- 
ment limitations, it is usually necessary to solve them 
one at atime, switching in the correct non- homogeneous 
inputs from the y, program in a suitable way. This 
leads to the use of the relaxation method. 


With the relaxation method we change only one para- 
meter, C;, until the component, (Vf);, corresponding 
to the parameter is zero or nearly so. The pro- 
cedure is applied to all the parameters on a cyclic 
basis. An example of the resulting path in a_ two 
dimensional parameter plane is 


Co 


C, 


This method obviously is ideal for use with the para- 
meter-influence equations. Nothing needs to be changed 
in the influence equation program except the non-homo- 
geneous term. 


In the event it is not convenient touse the parameter- 
influence method, successive values of f can be cal- 
culated for a sequence of values of Cj. This is con- 
tinued until two values of f bracket (Vf); =0. Then 
the procedure is applied varying C;,, and so on. This 
will alsolead tothe minimum for f. In either case the 
subscript indexing and switching of non-homogeneous 
terms is done with a digital logic program. In the 
latter case, successive values of f are retained by a 
two or three word memory chain and comparisons made 
by comparators. 


With either method the step size must be determined 
by a gain factor which ensures proper convergence. 


CHAPTER 21 
MEDICAL APPLICATIONS 


Originally published by: Beckman Instruments, Inc. 
Author: Hiroshi Hara 


ANALOG SIMULATION OF VENTRICULAR PUMP 
ACTION 


Due to the complexity and non-linearity involved, the 
building of an analog model of blood circulation is one 
of the challenging problems in biomedical science. 
This example is by no means intended to coverthe 
subject, but illustrates how iterative programming is 
conveniently used in simulating ventricular pump 
action. 


The ventricle has an inlet valve (tricuspid valve) and 
an outlet valve (pulmonary valve). A cardiac cycle 
consists of two periods called 'the diastole andthe sys- 
tole, as shown in Figure 1. During the diastole, the 
heart muscle relaxes and blood flows into the ventri- 
cle as the inlet valve opens. At the completion of the 
diastole, the systole begins. During this period the 
heart muscle contracts andas the ventricular pressure 
exceeds the arterial pressure, the outlet valve opens, 
expelling bloodfrom the ventricle. An electric analog 
of the ventricle is shown in Figure 2. The analogous 
quantities are indicated below. 


Blood Pressure ~ Voltage 
Blood Flow ~ Current 
Blood Volume ~ Charge 

nth CYCLE (n+!) th CYCLE 











DIASTOLE , SYSTOLE | DIASTOLE SYSTOLE 


FIGURE 1 


The contraction and relaxation of the muscular wall of 
the ventricle is simulated by decrease and increase, 
respectively, of capacitance c, in some prescribed 
manner. If we let E = Cy 
the pump action is assumed to be described by the 
following equations: 


= elastance of the ventricle, 


“t 
Ealn =A- (A-Es_1) e +B [aa] n (1) 


_ ad _ 
leg], = Eat qt ¥ = fay) (2) 


where 
Egl., = Elastance during nth diastole 


fe in ~ Elastance during nth systole 


EG = Elastance at the end of nth diastole 
E14 = Elastance at the end of (n-1) th systole 


[aa] n = Blood volume during nth diastole 


Yn = Slope of [E | 7 to be determined as a 


function of qo , the blood volume at the 
end of nth diastole. 


A, B, and 7 are positive constants. 


Since the static work done by the heart ata particular” 
value of ad is known experimentally (Starling's Curve), 
ve must be foundas a function of ac such that experi- 


mental results are duplicated. The work done during 
the systole is 


Tsystole 
Ws =k J Voiyodt, 


where Vo = arterial pressure 
iyo = blood flow out of ventricle 
k = constant 


and consequently for a certain ao ’ a must be such 
that [w.| ne [w.| Starling. In other words, itera- 


tion with respect to w, is required. 


The passive circuit shown in Figure 2, equations (1) 
and (2), and an iteration circuit are implemented as 
shown in Figure 3. Note that two iterative cycles are 
required to simulate one cardiac cycle. The distinc- 
tion between the diastole and systole is made by afunc- 
tion relay (Relay KO) which in turn is controlled by 
one of the memory integrators in a chain of four. The 
action of this circuitry is shown by the timing diagram 
in Figure 4. Since the compute intervals for the dias- 
tole and systole differ, the comparator which controls 
the compute interval of the iterative control program 
must, necessarily, require two different reference 
voltages which are switched alternately by relay KO. 
The function of the M2- M1 memory chain is to make 
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ELECTRIC ANALOG OF A VENTRICLE 


Ri2 | Ve Res 


| — | = 
| iol C,(t) it2 


VENTRICLE 


Vo : VENOUS PRESSURE 


V, : VENTRICULAR PRESSURE 


Vo . ARTERIAL PRESSURE 


CHARGE, q 


be 


Ro; :INLET VALVE RESISTANCE 
Rig OUTLET VALVE RESISTANCE 
C, : VENTRICLE CAPACITANCE = +,WHERE E =ELASTANCE 


q © BLOOD VOLUME 


Ws | WORK DONE BY THE VENTRICLE DURING SYSTOLE = 
kf vo iy dt, WHERE k = CONSTANT 


FIGURE 2 


the final value of q (blood volume) at the end of the nth 
diastole available throughout the nth systole. (M de- 
notes complementary memory.) Note that relay KO 
prevents this memory chain from learning the end 
systolic blood volume. 


In order to determine the required functional relation-- 


ships between ao and ¥, , a simple iteration circuit is 


provided. With some ¥,, at the output of M3, systolic 


work is computed, and this is compared to the output 
of FG 2 which follows Starling's curve. The error is 
repeatedly used to up-date ¥,, until the desired ¥,, is 


obtained. The above procedure is repeatedfor various 
values of v,, until enough information is obtained to set 
up FG1. With FG1 inthecircuit, a ventricular model 
is ready to be used as a part of a circulation model 
which will not be discussed here. 
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CHAPTER 22 
A PRACTICAL APPROACH TO ADAPTIVE CONTROL 


Originally published in: CONTROL ENGINEERING, May 1966 
Authors: J. W. Banham Jr. and W. L. Smith 


Controllers are linear elements that are often re- 
quired to operate in nonlinear systems. Thus they 
cannot be expected to provide optimum performance 
over a wide range of system operating conditions. 
However, through a linear representation of the non- 
linear system a controller can be designed.with adap- 
tive features that does provide optimum compensa- 
tion for transient system requirements. The follow- 
ing text describes the design of a controller and its 
conversion to the "adaptive" type (one whose dynamic 
parameters are continuously established as a function 
of transients inthe system). The approach described 
has been breadboarded and tested for a shipboard sys- 
tem but the design techniques are applicable to other 
control systems. 


The maneuverability of anaval ship depends to a large 
degree upon the transient response characteristics of 
the steam generating and machinery systems. Most 
marine steam generators use forced draft; air for 
combustion is supplied to thefuel oil burners by vari- 
able speed blowers. The forced draft blowers in the 
machinery plant of a high speed combat ship are gen- 
erally powered by direct drive steam turbines and 
regulated by controlling the steam flow to the turbines. 
The rotational inertial characteristics of these ma- 
chines, over a wide range of operating speeds, contri- 
bute materially to the 'sluggishness" of process re- 
sponse to changes inappliedload. Designers interested 
in improving the system response must optimize the 
controller parameters of those system loops or sub- 
loops with the longer time constants. 


SYSTEM DESIGN STUDY 


In a conventional automatic combustion control sys- 
tem, controller settings are usually determined em- 
pirically at the load condition of maximum open loop 
sensitivity. Because controllers are linear elements 
and the forced draft blowersare nonlinear, the situa- 
tion is considerably less than optimum over a wide 
range of operating loads. 


a) Nonlinearities in the blower 


For any type of rotating machinery, the angular ac- 
celeration of the rotor assembly is proportional to the 
applied differential torque and inversely proportional 
to the polar moment of inertia. This linear relation- 
ship may be expressed by the simple differential equa- 
tion: 


(T, 7 T,.) 


Cont | be 


where 
“ 2 
w = angular acceleration, rad/sec 
J = polar moment of inertia, Ib-ft-sec” 


JT = driving torque developed by the prime 
mover, lb-ft 


T.. = load torque imposed by the blower, lb-ft 


Several factors influence the applied torques. The 
torque developed by a steam turbine is a function of 
the turbine steam rate and its rotational speed. A 
typical turbine performance graph illustrating the func- 
tional relationship is shown in Figure 1. The para- 
bolic curveonthis map represents the design load re- 
sistance for the turbine-fan unit. For this particular 
set of performance data, it was found practical to rep- 
resent the driving torque by the algebraic expression 


T. = os T -m Ss - 
. Cgmax re Gomax | \@max 
where 
G_/G = turbine steam flow, fraction of 
s’ ~smax : 
rating 

T = maximum locked rotor torque, lb-ft 
max 

m = aconstant, lb-ft 

w/w = angular velocity, fraction of rating 


max 


The steam flow to the turbine nozzles varies with the 
steam supply pressure and temperature and the steam 
valve flow coefficient and port area. The forced draft 
blowers considered in this study were equipped with 
a V-ported regulating throttle valve for which the 
characteristic flow-lift curve is as illustrated in Fig- 
ure 2. The supply steam state is reasonably constant 
throughout the load range. Theload torque varies di- 
rectly with the blower system resistance and with the 
square of the rotational speed. 


2 
T = = 
ro » (2nax) 








LINEAR REPRESENTATION OF THE 


NONLINEAR SYSTEM 





All of the elements of the combustion air flow control 
loop of amarine steam generating system are charac- 
terized by some nonlinearity. The following analysis 
illustrates one method of representing the forced draft 
blower nonlinearities by a set of linear functions. 


In the basic acceleration equation 


te 


1 
o = > (T, - T,) 


r 


the torque terms are functions of the turbine steam 
flow, Gg and the turbine-fan speed, w, as shown in 
the text of this article. The accelerating torque, Ty 


was found to be afunction of both steam flow and speed; 
the graphical relationship appears in Figure 1. In this 
study, the analytical relation between these variables 
was established by curve-fitting to manufacturer's 
performance curves, but can also be developed by a 
detailed derivation of momentum equations applied to 
the turbine wheel. For a fixed blower air system re- 
sistance, the load torque varies with the square of the 
rotationalspeed. The constant of proportionality, 6 , 
is established by this resistance. 


Steam flow to the turbine is regulated by a throttling 
valve characterized by the arbitrary function 


Gu. F(L) 


where L represents valve lift. 


Sufficient information is now available to permit the 
development of a set of transfer functions which de- 
scribes the open loop frequency response of the tur- 
bine-blower unit to small perturbations about a fixed 
operating level. The general perturbation response 
is given by 


ws oh 
Aw =F for, - AT, ,)dt 


If asmall disturbance is applied at a particular steady 


state initial condition (L Goo Wo)s the nonlinear 


torque terms may be reduced to linear functions ex- 
pressed by the following equations: 


AT, 


AG (T G 


- m 
max so G5 5%) 


AT. 2pw Aw 


Gs is a function of valve lift as shown in Figure 2; a 
tangent constructed at the point (Ly em) has aslope 
given by 


Substitution of these linearized functions into the ac- 
celeration equation and applying a Laplace operator 
yields 


~i}a | 
OUST a lat PO) p= n, MTmaxSso 


mG_w) - 


soo 240,00 


For any particular steady state operating level the 


terms d/dL F(L) IL=Lp, Trax Sso? MGso% and 26 w, 


are constants; the above equation may thus be expressed 
as a transfer function: 


d 
dL max~so _ mG,” | 


1 + [J/280 |S 


F(L)[T G 
Aw 
a2 (s) = 


This transfer function defines the response of shaft 
speed to valve lift for small perturbations given initial 


conditions defined by Lo Goo and Wo: 
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65 /Gs moy STEAM FLOW, FRACTION 
OF RATING 


T/T max Torque, fraction of rating 





2 0.4 06 08 1.0 
w/e mox Speed, fraction of rating 


Figure 1. Forced draft blower turbine performance 
curves. 





Gs/Gs mo, Steam flow, fraction of turbine rating 


Oo 0.2. 0.4 0.6 0.8 1.0 
L/L mox Valve lift, fraction of full stroke 


Figure 2. Turbine valve steam regulating charac- 
teristics. The arbitrary function G, = F (L) used in 


Figures 3 and 7 is illustrated. Construction of a tan- 
gent at OP Sau 0. 30, L/L, ax 0-53) illustrates 
the graphical approachtothe linear representation of 
a nonlinear system. The tangent slope represents a 
valve gain of 1.12 percent of steam flow per percent 
of valve lift at a point corresponding to 53 percent of 
lift. 


where g isa constantforaparticular system resist- 
ance, lb-ft. In the design study, the coefficient 8 
was permitted to vary with the effective fan discharge 
area (established by the number of burner air regis- 
ters in use in the steam generator). 


This informationled to the model of the system shown 
in Figure 3. The problem at this point was to devise 
a control law which would provide optimum compensa- 
tion for the nonlinear system over a wide range of 
operating speeds. For this design, the system is in- 
itially described in the frequency domain as a contin- 
uous series of linear systems, each similar in form 
to the others but having different parameters. 


b) Linear representation of the nonlinear system 


One approach to the linear representation of a non- 
linear system is to expand the nonlinear elements of 
the system by means of the Taylor series, then dis- 
regard all but the linear terms. This produces a 
linear system whichapproximates the real system re- 
sponse to small perturbations about a fixed operating 
point. Because the steam admission valve character- 
istic is an essentially arbitrary function, a graphical 
approach was used. In this method, all of the non- 
linear terms are represented as linear functions by 
constructing tangents at selected operating points along 
the curves which represent the functions of the alge- 
braic variables. The slope of each tangent curve is 
taken as the gain of the element under consideration, 
Figure 2. The resulting linear model is illustrated in 
Figure 4. The development of the values of the pa- 
rameters associated with this model is detailed in the 
box on page 22-2. N 


A linear model like the one in Figure 4 leads to a set 
of transfer functions instead of a single function. 
Each function in the set represents the open loop fre- 
quency response of the system to a small signal si- 
nusoidal disturbance superimposed on a particular 
steady state initial condition. The validity of the 
analytical model constructed in this fashion was con- 
firmed by comparison with actual frequency response 
data obtained from the operating units subjected to 
pulse test analysis at two load conditions, (Ref. 1). 
The family of transfer functions thus obtained was 
used as the basis for a series of control system de- 
signs ‘based on classical frequency response tech- 
niques applicable to linear systems. 


DESIGN OF THE CONTROL SYSTEM 


The following conditions were used as criteria for the 
dynamic performance of the closed-loop system: 


a. There must be no offset error in the steady 
state. 


b. No overshoot in the transient response to a 
"ramp" input disturbance is to be permitted. 


c. The time integral of the system absolute error 
is to be minimized. 


The conventional control loop employs aproportional- 
plus-integral controller which receives as its feed- 
back the output signal from a combustion air flow 
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DEVELOPMENT OF CONTROLLER DESIGN 


Conventional controllers used for regulating forced 
draft blowers are simple proportional-plus-integral 
devices characterized by fixed gain and reset rate 
over the entire operating range. The output signal 
of a linear controller of this type is given by 


1 
e=veat f eat 


error between input and feedback signals 
controller proportional gain 


controller reset rate 


An adaptive feature was incorporated in such a con- 
troller by servoing the proportional gain and reset 
rate (integral crossover frequency) settings as re- 
quired to maintain an optimum combination of closed 
loop system response and stability over the designed 
range of forced draft blower operation. These con- 
troller parameters are shown plotted against various 
loop operating levels in Figure 6. Curve fitting of 
these graphical functions produced the following alge- 
braic expressions: 


Y 


w/w 
/® max 


variable controller gain 


optimum controller gain at w = Oa 


steady state loop operating level 
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1 : 
where ie variable controller reset rate, rad/sec 


i = optimum controller reset rate 
oO 


Ww = 
at w/ ae 0 


= optimum controller reset rate 


at w/ sax = 1 


The output of the adaptive air flow controller may 
therefore be expressed by the equation 


1 1 
a= 7 ( Je +t feat 
w/ To 


“max 


1 1 


Ee Goa) 


T 
max max 


It is to be noted that the controller parameters have 
been expressed as functions of the control feedback 
signal (in this particular system, the blower rota- 
tional speed). To produce maximum loop gain during 
any transient disturbance, the controller gain and 
reset rate are automatically varied as a function of 
either the command signal ¢@ or the feedback signal 
w, Whichever of the two is least. An arrangement which 
produces this effect is illustrated inthe computer sche- 
matic, Figure 7. Using the notation indicated on this 
diagram, the controller equation becomes 


= y 


a) ey 


Pmax oO 


eg) fo 


Tmax max 








Figure 3. Forced draft blower block diagram. 





Figure 4. Block diagram of forced draft blower perturbation model at (Lo : wo) 


transmitter. For the "adaptive" design, it was elected 
to utilize the same basic controller concept but to vary 
the proportional gain and integral cross-over frequency 
to maintain optimum performance. 


The conventional design was initiated by selecting con- 
troller parameters which would give the entire control 
loop a gain margin of six db and a phase margin of 35 
deg; these margin goals were derived from previous 
experience. This conventional procedure was repeated 
for each of the forced draft blower transfer functions 
at the various operating load levels. The result, de- 
tailed in the box on page 22-4, was aset of empirically 

optimum controller settings which could be correlated 


with the loop input or output to develop analytic func- 
tions that would specify the adaptive nature of the 
design. 


The air flow controller thus obtained was simulated 
together with the forced draft blowers and air flow 
system onan analogcomputer. The transient response 
to input disturbances comparable to the effects of ship 
maneuvering were then studied. Typical results are 
shownin Figure 5. Note that a small amount of over- 
shoot occurred which violates one of the system de- 
sign criteria. This observation led to a repetition 
of the design procedure, this time using gain and phase 
margins of 9 db and 45 deg, respectively. Figure 6 
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TIME, SECONDS 


Figure 5. System transient response with initial 
optimization criteria. 


Controller reciprocal gain, |/y 


Controller reset rate/gain, 





Loop operating level, w/w max 


Figure 6. Controller parameters (servoed controller 
settings) as a function of loop operating level in the 
steady state. Open loop gain and phase marginsiare 
9 db and 45 deg. 
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illustrates the relationship between the controller 
parameters and the loop operating level for perform- 
ance with the revised open-loop margins. 


A revised analog computer program was then pre- 
pared; a simplified version of the basic program is 
shown in Figure 7. Transient response of this simu- 
lated system appears in the chart recordings of Fig- 
ure 8A; the corresponding response of the simulated 
adaptive control system appears for comparative pur- 
poses in Figure 8B. It is evident from these record- 
ings that the adaptive system is considerably superior 
to the conventional system. Inspection of the error 
recordings discloses a reduction of five to one in the 
time integral of absolute error. At the same time 
there is neither overshoot nor permanent offset. Thus 
the design criteria have been met. 


FROM SIMULATION OF ACTUAL TESTS 


Results obtained from the analog computer studies 
supported a decision to breadboard a pneumatic con- 
trol system from standard hardware. The basic com- 
ponents used in the adaptive controller were pneu- 
matic proportional and integral amplifiers with servo- 
driven gain settings. Two of these units provided the 
necessary proportional-plus- reset control mode. The 
gains of these two devices were automatically varied 
as a function of either the input to the closed air flow 
control loop or its output, whichever produced the 
greatest loop gain during any transient disturbance. 
In the steady state, these variables are equal and 
produce the same gain in the open loop. 


Tests of the loop were conducted by imposing "ramp" 
inputs which corresponded to those used in thecom- 
puter simulation, thus aiding in an evaluation of.con- 
trol system performance. The recorded results of 
one of these transients are shown in Figure 8C. The 
transient performance of the closed loop system 
closely paralleled that of the simulated system, see 
Figure 8B. 


The air flow control loop comprises one of the open- 
loop elements of the steam generator's automatic 
combustion control system. The principal benefit of 
the adaptive air flow control] system when used as a 
part of the overall steam generator control system 
is its ability to reduce the deviation of controlled 
steam pressure and to eliminate the possibility of 
incomplete combustion (and accompanying smoke) 
during maneuvering. Although hardware requirements 
for systems of this type are presently more com- 
plex, and less reliable than for conventional systems, 
adaptive design techniques become attractive when 
applied to systems constructed from encapsulated or 
modular design control circuitry. The U. S. Navy is 
currently investigating adaptive systems of the solid 
state pneumatic (pure fluid) and solid state electronic 
types for use in the control of shipboard main pro- 
pulsion and auxiliary machinery. 


PROPORTIONAL PLUS 
eg SUMMER RESET CONTROLLER 


SIGNAL 
LIMITER 





Figure 7. Simplified analog computer program used to study conventional 
and adaptive controller responses. 
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Figure 8. Transient response recordings. A-simu- 
lated conventional controller. B-simulated adaptive 
controller. C-breadboard adaptive controller. 
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Glossary of Abbreviations 


Capacitor 
Compute 
Clock 


Check 


Dor DN Down 


DCU 


DFG 


DVM 


ES 


Fb 


FF 


Decimal Counting Unit 


Diode Function Generator 


Digital Volt Meter 


Electronic Switch 
Feedback 
Flip-flop 


Forward-Reset Real Time 


Function (Operational) Relay 


Function Switch 
Hold 

Initial Condition 
Initial Junction 

Input 


Inverter 


NO 


O 


Iterative Operation 
Junction 

Relay 

Logical Variable 
Meter 

Normally Closed 
Normally Open 


Output 


Por POT Coef. Potentiometer 


R 


R 


R-H 


RO 


Reset 

Resistor 
Reverse-Hold 
Repetitive Operation 
Reset-T oggle-Set 
Set 

Summing Junction 
Toggle 

Trunk Line 


Up 
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The symbology used throughout this handbook is in conformity with the objectives 


of Simulation Councils, Inc., namely, to encourage the use of uniform analog com- 


puter graphics. This article was prepared by SCI, and appeared in the March 1966 


issue of SIMULATION. 


Permission has been granted to reprint this article. 





Presenting... 


Uniform graphics for SIMULATION 


The symbols and the methods of laying out analog 
and hybrid computer diagrams presented here are 
advocated to alleviate the confusion caused by the 
uncoordinated invention of new symbols and dia- 
gramming practices. The increasing use of hybrid 
techniques and equipment has aggravated an already 
bad situation to the point that it is often no longer 
possible for one worker in our field to read another’s 
diagram. Usually this is because symbols are devised 
and diagrams are drawn to include details peculiar 
to a particular kind of equipment. Such a wiring, or 
“patching,” diagram is of course necessary for set- 
ting up and checking out an actual simulation, but 
hardware-peculiar details are only confusing to those 
with other kinds of equipment. With few exceptions, 
the use of a simplified signal-flow diagram to illus- 
trate technical articles is much more effective. 
With the foregoing in mind an SCi committee 

composed of _ 

GEORGE BURGIN 

JOE HUSSEY 

HANS JORGENSEN 

GRANINO KORN 

JOHN McLEOD 


selected the symbols and offers the following sugges- 
tions for their use. Primary considerations were cur- 
rent usage, clarity, and simplicity. We devised no new 
symbols and, unless there were overriding indica- 
tions to the contrary, we adopted those already in 
widest use. Clarity and simplicity, we believe, will be 
enhanced by the choice of unique shapes to repre- 
sent different components, and the elimination of all 
unnecessary details in diagrams. 

There was no intent on the part of our committee 
to set up standards for the industry. However, all 
diagrams appearing in SIMULATION will be pre- 
pared according to the committee’s recommenda- 
tions (as they may be modified from time to time), 
and we hope that these recommendations will prove 
attractive to others. Suggestions for modifications 
and additions are solicited. 


General rules 

The following methods and symbols are recom- 
mended for the illustration of technical articles pre- 
pared for publication. Unless the purpose of the 
article is to describe. the use of a particular kind of 
equipment, and the hardware details are pertinent 
to the subject, such illustrations should not be “hard- 
ware-peculiar.” In other words, the objective should 
be to show signal flow, rather than “patching” 
details. 

The primary, or overall, system diagram should 
show only the essential signal flow. Where it is neces- 
sary to show details, separate diagrams should be 
made and referenced to the primary diagram by en- 
closing the detailed area of the primary diagram in 
dotted lines with suitable notation. 

The direction of signal flow should be indicated 
by arrowheads except where the shape of the sym- 
bols makes the direction of flow obvious. Primary 
signal flow (with the exception of feedback loops) 
should be from left to right, and, if practical, each 
“line” of symbols should be made to read like the 
mathematical relation it represents. 

The choice of whether to end a line and label it 
(preferably with the symbol of the variable that the 
signal represents) when it reaches the right-hand 
side of the diagram, and then indicate its continua- 
tion with the same label as it enters again at the 
left-hand side, instead of drawing in the connection, 
should be made on the basis of clarity; if a line re- 
turning the signal from right to left will cross many 
other lines and be hard to follow, it should not be 
drawn in. 

If a diagram involves a number of identical cir- 
cuits, only one should be shown in detail, while the 
others should be indicated by boxes with appropriate 
notation. 

Components should not be numbered unless they 
are referred to by number in the text. 

All amplifier gains should be shown just outside 
the amplifier at the point where the input enters. 
Unity gains should not be labeled. 


Always apply the test of clarity and simplicity. Ask 
yourself: “Is this the most understandable way to 
diagram this for those unfamiliar with the hardware, 
and less familiar with the subject, than 1?” 


SIMULATION 37 








SIMULATION 














© 
© 
a 
YS 
= 
UO 
< 
= 








